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 MATHEM 
_  INGENIEURWISSENSCHAFILICHE FORSCHUNGSARBEITEN 


.  „Volumentreue Verzerrungen bei endlichen Formänderungen. 
Von R. Moufang in Essen (jetzt Frankfurt a.M.). | 


74 BR eye ‚ Der aus dem Hookeschen Gesetz bekannte Übergang vom. Verzerrungsiensor zum Verzerrungsdeviator 


im Rahmen infinitesimaler Formänderungen wird auf endliche Formänderungen übertragen, indem von 
dem Verzerrungstensor ein geeigneter Kugeltensor abgezogen wird; er entsteht aus dem Einheitstensor durch 


ERDRLER mit einer Funktion der drei Tensorinvarianten, die als Lösung einer kubischen Gleichung 
esiimmt ist. ’ j s 


in case of infimitesimal deformations, is transferred to finite deformations by subtracting an appropriate 
spheric tensor from the distortion tensor. The spheric tensor is obtained by multiplying the unity tensor 
with a factor that is, as solution of a cubic, a function of the three tensorial invariants. 2 


Le passage du tenseur au deviateur de deformation dans le cas de deformations infinitesimales, connu 

par la.loi.de Hocke, est &tendu aux deformations d’ordre fini en separant dutenseur de deformation un tenseur 

‚ sph£rique approprie; ce dernier, ü son tours, est deduit du tenseur elementaire en le multipliant par une 
fonction des trois invariants tensoriells, qui doit repr ‘senter .une solution d’une equation cubique. 


llepexon oT Tensopa gehopmannmu K KeBHaTopy nedopmaunn, B CILyuae ÖecKOHeyHO MAIBIX 
- nebopmannf usBecTHpä 13 3akoHaT'yra, pacıpocrpauserca Ha cıyyal KOHeYHBIX TepopManai. 
- Ma Tensopa nehopMalun BhigenneTch NuHaTOHANBHBIH TEeH3OP, TONYYAMINHÜCH U3 ENUHUYHOTO 


IFTFÜRANGEWANDTE 
ATIK UND MECHANIK 


The transition from the distortion tensor to the distortion deviator, well knoun from H ooke’s law Ä 
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TeH30pa YMHO:KeHHM HA PYHEIMI TPpeX HHBAPNAHTOB TEeH30Pa MeoopMauun. ITa hyHruna 


YAOBNETBOPAeT HEKOTOPOMY KyOnyecKoMy YPaBHeHHI. 
Re Einleitung. 
Die Erforschung der Näturgesetze, nach denen die plastische Formänderung sich vollzieht, 
erfordert zunächst begriffliche Untersuchungen, die die mathematische Definition der durch 
die Verformungsgesetze miteinander zu verknüpfenden Größen betreffen. Die grundlegenden 


Begriffe sind in diesem Falle der Verzerrungs- und der Spannungstensor. Die Definition dieser 


Begriffe für endliche Formänderungen ist in der Literatur bereits erledigt?). 


Bei der Frage nach der möglichen Form der Plastizitätsgesetze?) liegt es nahe, in Analogie 
zum Hookeschen Gesetz zunächst gestaltändernde von volumenändernden Verformungen 
zu unterscheiden. Beiinfinitesimalen Formänderungen führt dies auf den Verzerrungsdeviator, 
der aus dem Verzerrungstensor dadurch entsteht, daß von den Elementen der Hauptdiagonale 
das arithmetische Mittel aus diesen Elementen abgezogen wird. Der so entstehende Tensor 
bringt keine Volumenänderung’mit sich, weilseine erste Invariante, durch die bei infinitesimalen 
Formänderungen die Volumendehnung gegeben ist, verschwindet. Beiendlichen Formänderungen 
ist die Volumendehnung nicht nur von der ersten Invarianten des Verzerrungstensors abhängig, 
sondern von allen drei Invarianten. Dies hat zur Folge, daß die Bildung des Deviators nach der 
obengenannten Vorschrift jetzt nicht mehr einen Verzerrungstensor mit ver- 
schwindender Volumendehnung liefert. Das Ziel der folgenden Betrachtungen ist die 
Konstruktion eines ‚‚Verzerrungsdeviators‘ bei endlichen Formänderungen, d.h. eines Tensors, 
der keine Volumendehnung bewirkt. 


$1. Die Definition des Verzerrungstensors. 
Der Vollständigkeit halber wird die Definition des Verzerrungstensors bei beliebigen 
Formänderungen hier noch einmal angeführt. 
Durch die Belastung erfahren die Punkte des elastischen Kontinuums, bezogen auf ein 
Koordinatensystem (x!, «2, x?) eine Lagenänderung, die beschrieben wird durch 


Mae Meiiii rar. 


1) Siehe z. B.E. Tre ff tz, Z. angew. Math. Mech. 13 (1933), S.160; F. D. Murnaghan, Amer. J. Math. 
59 (1927), 8.235; R. Kappus, Z. angew. Math. Mech. 19: (1939), S. 271, vgl. auch R.Moufang, Ing. Arch. 12 


1941),.8. 265. 
2) Eine Untersuchung hierüber im Rahmen unendlich kleiner Formänderungenfindetsich beiH. From m, Ing.-Arch. 4 


(1933), 8. 432. 
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‚Das Längenelement im Punkte P ist gegeben durch 


de = Egır(® (2) da'da* a ra r we, Se (2) 
,k= : 
GmeeprScheng das Längenelement im Punkte P’ durch 

der I gınle DE ee ee (2). 


Trägt man in Gl. (2) rechts für « die sich durch Auflösung von Gl. (1) ergebenden Ausdrücke. 


in den x’ ein, so wird: 
de — 5 Orle «') da" da" (ri = Gin). 


Die Koeffizienten S;; der quadratischen Differentialform 


3 2 \ 
ds’ — ds? = I (2(8) — Girle')) da da’ 
uk=1 DET ee (3) 


=28 8;:(® dad", a ) 
heißen die Komponenten des Verzerrungstensors in kovarianten Komponenten. Hierbei sind 


willkürlich die Koordinaten des Punktes in seiner Endlage als unabhängige Variable gewählt. 


Sollen die x als unabhängige Variable gewählt werden, so erfahren die folgenden Ausführungen 
eine ohne Schwierigkeit durchzuführende Abänderung, die gestützt ist auf die Definition des 
Verzerrungstensors (7;z(2)), wo jetzt 


ar dee N Tnleerae ren (3). 


E) 
ist. 


$2. Die Volumendehnung. 
Mit der Transformation Gl]. (1) ist eine Volumendehnung verbunden 
Al en en. (4) 
aV Ya“) : 


wobei //g:r// die Determinate der Matrix (g, Ist (gi*) die zu (g;r) inverse Matrix 
und setzt man 


Pr ER 


ı=1 
so wird nach dem Multiplikationssatz für Determinanten 


; Eülie xl! 
Ne Hr) elf 

IR au girl 
also wird 

dV : 

EL en (5) 
Wegen 
+ DIS Nee G;: 
ist 

= 2 gt gu— 2 NS = —28 an 

h 2 1 ı—k 


also hat man 
Na EI EI BT, 
J,—= 8 +84 = Summe der Hauptdiagonalelemente von (8!) 


= ae (8185 — 51 8?) = Summe der Hauptunterdeterminanten 6 
2. Ordnung von (S®) RE un (6) 
J, = //8}/|= Determinante von (8%) 


und nach Gl. (5) 


av k } k n 
Zr = Yld-28= Yı 29, 448% ae ee (7% 


| 
E 


Aus Gl. (8) erkennt man dagegen, daß bei e 

„daß (S;x) ein er ist, nicht tolet, daß 
on solcher r ist, denn. mit EL ELITE 
Er j EEE 2 SE nz 


ea * 


4 


Be X 2 RE RE e ad, 2a, ee FR IN TE $ 5 


Mr eree 


Pa g 3. Volumentreue ee hei endlichen Forint 


RS Gegeben sei ein Verzerrungstensor SER ‚gesucht ist Rn nr 1 Jr I AR A 
E so daß er Tensor ae ma . Y | 


gem a g A - a 1 PR N - . | Be i ; x h 

ao ‘ i 3 = 4 ; ae Sr ii a Te 100 a ae es Es Ber Di 
E ER eine EHEN Slumendehnung ee pist eine skalare Funktion der drei Invarianten 2 
E=,.- von ©. Der Einfachheit halber beziehen wir uns im folgenden auf ein grihopaiade Koordinaten- Ben’ 
IE a | 2 ee 
Be ge ie = 
8 Beer so Habe es möglich ist, von den kovarianten Tensorkomponenten zu den physikalischen Tensor- 
Br. Ft komponenten Sir überzugehen vermögen: 
£ = | Ss; . (10). 


| Hierbei geht der Fundamentaltensor & in den Einheitstensor € über. Man a leicht nach, 
daß die Invarianten von ($;;) den entsprechenden Ausdrücken, gebildet mit den S;; an Stelle 
E- E _ von Si, gleich werden: ; | 


In 


SEEN _ Damit wird nach GI.(7) . 


der Hauptunterdeterminanten 2. Ordnung, die Determinante des Tensors S=(8;,) bedeuten. 
F Die Bedingung, daß Weine verschwindende Volumendehnung besitzt, lautet nach Gl. (7) mit W 
g - an Stelle von © 


= Villen 2Bu/l; 0 für ;Fr- ET 
J | er Übergang zu physikalischen Tensorkomponenten erhält man also aus Gl. (9) die ‘Matricen- 
F gleichung gan er ö E 
® 17, AS ZylhS IN E nn), P: 
RER worin J,Jg Js der Reihe nach die Summe der Hauptdiagonalelemente, die Summe 


E:;-- BIN Se 
# oder ausgerechnet ‚ 
3 BR EN 8 I, = 13 EA 
Setzt man 


Ve a Re (18) 
F ATTERSEE N, 
so ergibt sich aus Gl. (12’) für A die kubische Gleichung 
a ee, 
14 
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Bei infinitesimalen Formänderungen, wo J; und J, gegen J, vernachlässigbar sind, wird in 
erster Näherung nach Gl. (14) ' 


J, 
also ; 
J% J! 
RI 


was nach Gl. (9) auf die bei infinitesimalen Formänderungen bekannte Deviatorformel führt 
BEE a ; 
YA=-65—-Z€E=E. 
2. Bei endlichen Formänderungen hat man die Gl. (14) aufzulösen. 
Mit 
124%, o=-2.+— 1 re ET 


geht sie über in die Normalform 
23 +3p2 +2q9=0 


12 

1 ji 

RE: ee (16). 
DR IL 

ee a > 


3p=— 


Die Discriminante ist 


AD 
[ER HEINE III) 
und es ist : | 
EN EBEN ee m 


In dieser Formel sind rechts nachträglich für die J’ wieder die Größen J einzusetzen. 


3. Die Diseriminante unterscheidet sich von zg um Glieder, die in den Verzerrungen min- 


destens von 3. Ordnung sind. Beinicht allzu großen Verzerrungen ist also 4>0, — dies ist z.B. 
sicher dann erfüllt, wenn dritte Potenzen der Verzerrungen gegen die ersten Potenzen vernach- 
lässigt werden können —, d.h. Gl. (16) besitzt genau 1 reelleLösung, die durch Gl. (18) dargestellt 
wird. Unterliegen die J/, J,, J| bzw. die Verzerrungen von vornherein hinsichtlich der Größe 
keiner Beschränkung, so kann A<O werden, so daß von vornherein die Möglichkeit nicht 
auszuschließen ist, daß Gl. (16) mehr als eine reelle Lösung hat, d.h. daß die ‚„‚Deviatorbildung‘“ 


Gl. (9) zu gegebenem © nicht eindeutig ist. Z.B. für 


NeBEan, Seo ne 
mit 
= 3 2 2 Rn, 
= oe Sa 0 = S;, Sr =0(ik), 
was durch 
2 R 
a: ads, ER 


realisiert wird, ist A <0. Es gibt also bei endlichen Formänderungen zu einem gegebenen Ver- 
zerrungstensor stets mindestens einen ‚‚Verzerrungsdeviator“ ; die Möglichkeit, daß die Deviator- 
bildung mehrdeutig (höchstens dreideutig) ist, kann von vornherein nicht ausgeschlossen werden. 
Dies bedarf bei der Aufstellung von Verformungsgesetzen noch einer besonderen Untersuchung, 
auf die hier nicht eingegangen wird. Die Mehrdeutigkeit entfällt aber von vornherein in einem 
gewissen Bereich der Variablen J,, J;, J,, der den Nullpunkt enthält, also für einen auf Haupt- 
axen bezogenen Verzerrungstensor den Fall verschwindender Hauptdehnungen umfaßt. 


lag GL = zu N 
unktionen erfordert, was im einzelnen Falle 


lism der Naturgesetze wenig handlich machen wird. Es lie Bahr De zu fragen, 
n kommt, wenn man sich in Gl. (9) auf ganzrationale Funktionen g beschränkt. Ganzrati 
erungsausdrücke für @ ergeben sich aus Gl. (18), wenn man die auftretenden Wurzel 
‚der: -Cardanoschen Formel nach dem Taylorschen Satz entwickelt. Diese Rechnung ist a 
wenig übersichtlich. Einfacher ist es; direkt mit einem Reihenansatz | in di, Kuh & 
EN für p hineinzugehen, wie das jetzt näher ausgeführt wird. Pr Ä ee 
2. Beschränkt ER sich i in pi auf Glieder, die von 1. Ordnung i in den Verzerrungen, 3 


so liefert. Gl. a2) 


a 


Er zu Mi Tone. a. ‚(9 )] Behönige Volumendehnung ist: dann « eine Funktion der Verzerrungen, Fi 


u; die mit den Gliedern 2. nr beginnt. nr BET 
Br. erh ae . == ee i 


3 & Een liefert in Gl. 112°) eingesetzt bei Vernachlässigung aller Glieder, die in En Verzerrungen. von 
Be 3. Ordnung sind, 


2 : ä EB: iz +7 5 Ber: 

B 215070 ee 

5 y ! er rn j REDE 1 ee . 2 
E. : = = (4 +21) N ee (19). 


Zu diesem Br gehört eine Yolumendebnung. die in den Verzerrungen von 3. Ordnung ist. 


3. Der Ansatz 
E e ’ 
a ges tghracdı | \ 


+4J3 +@3J1J2 +03J3 


liefert 
16 8 8 
es: 3 GR = +7: 
also ist 
| re] 2 VER ERE 4 4 er 
1 8-(SAtgN- tt gar). en) 


2 ein Tensor, dessen Volumendehnung in den Verzerrungen von 4. Ordnung ist. a 
4. Dieses Konstruktionsverfahren läßt sich in eindeutiger Weise beliebig fortsetzen. In e: 


der Tat, man setze 
Vene Far KR E- 


wo @% (J1,J2,J,) homogen von der Dimension » in. den a ist. Die 
Gl. (12°) ist zu erfüllen bis einschließlich der Glieder von der Dimension (n+ 1), und zwar iden- 
tisch. Sei für alle Glieder der Dimension n bewiesen, daß die Koeffizienten in den 9, bis @n 
eindeutig bestimmt sind nach der Methode der Koeffizientenvergleichung, so erhält man für 


die Glieder der Dimension (rn +1) aus Gl. (12°) 
St tP (II J3; 9 4 9) = 0: 
wo ® eine bekannte Funktion ihrer Argumente ist. In der letzten Gleichung stehen links mit 


de N Arche eingegangen, da sie für das vorliegen eP 
} en: die kleiner als 1 sind, kann man in der Praxis eine geeign 

rernachlässigen, so daß das eben geschilderte Verfahren eine ‚bra 
"Verzerrungsdeviator liefert. Bei beliebigen Verzerrungen steht. von v 
RE mit eulie von Gl. Be und 15) zur I 


"ber: en Formatierungen? 


EI gr der Spannungsdeviator, & der te) een das Sc 
ee Pr 


Ber NEUEN) nn 
‘wo f von der Beziehen Wahl der Bssaie nicht Me und (® Je bzw. (&)2 das Sara 


= Produkt von ®’ bzw. &’ mit sich selbst ist. Bezieht man den Verzerrungstensor auf seine Haupt- - 
. axen und bezeichnet, ‚seine Komponenten dann wie üblich mit &7,.&» Eu so lautet ur (22) bei 


pr 


taten, ; Ye RE, = 


infinitesimalen Formänderungen, wo & =6— 3 or, a re 


K o? Fe 4o3 — 0109-03 men 0,) + y2 (e? ei ra HN &2 — &2 men) 2a). 


. Die o; sind die wahren en ar &; hängen mit den Verschiebungen und hr ER 


A ingen in linearer Weise zusammen. Bei endlichen Formänderungen hat man A aus Gl. 9) 
an Stelle von © zu setzen, also 


2 MSN +TP & 
B> | E ws 3 de Bee 
Be: SH eh eye | A 
Er | =-iTrEB TE Tas A)(&ı +82 +8) np UssHR B= 

Bi 

Be wo A eine Wurzel von 

SE 2—2(, +, te)? +4 +. +5 Ay =l.... (12*) Er 

& ist. Das Verfestigungsgesetz lautet dann t 


VRW= Seh tet Hei FON le, Fo Fe) FELZAP) . . . . (2b). 


Hierin sind jetzt &; die Verzerrungen bei endlichen Formänderungen, die eingangs definiert 
worden sind und von den Verschiebungsableitungen quadratisch abhängen. Die unabhängigen 
Variablen in den e; und o; sind die Koordinaten des Punktes in der deformierten Lage, und 
die Spannungen sind auf den verformten iR bezogen. Für den Grenzfall ' 


i=1 In (eı #82 +83) 
geht die rechte Seite von Gl. (22b) in die rechte Seite von Gl. (22a) über. 


Hingegangen: Februar 1945. 
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"30 Non  Quade, Zur Theorie der ebenen stetigen Gaswelen 


_ Zur Theorie der ebenen stetigen Gaswellen von 
endlicher Schwingungsweite. 
Von MW. Quade in Karlsruhe. 


Die Theorie der ebenen Gaswellen von endlicher Schwingungsweite gewinnt an Einfachheit, wenn als. 7 


abhängige Veränderliche die longitudinale Verschiebung eingeführt wird. Die letztere genügt einer Monge- 
Ampere schen Differentialgleichnng, deren Integral mit Hilfe der Methode der C'harakteristiken sowohl für 
den Fall isothermer als auch für den adiabatischer Zustandsänderung in geschlossener Form angegeben 
‚werden kann ohne daß eine Quadratur erforderlich ist. Die Riemann-Hugoniotsche Erscheinung 
Jindet dabei ihren geometrischen Ausdruck in dem Auftreten einer Rückkehrkante in der Integralfläche 
‘der Verschiebung. 


. The theory of plane gas undulations of finite amplitude becomes more simple by setting the longitudinal 
displacement as dependant variable. It satisfies a Monge-Ampere differential equation, the integral. 


of which can be given in a closed form by the method of characteristics not only in Ihe case of isothermal 
but also in the case of adiabatic change of state-without any integration. The Riemann-Hugoniot 
phenomenon hereby finds its geometrical expression in the appearance of a return edge within the integral 
surface of the displacement. . “ e 


. La theorie des ondes planes de gaz, & amplitude finie prend de simplicite en etudiant le deöplacement 


longitudinal comme fonction des coordonnees de Lagrange. Le deplacement satisfait & une equation dif- 


ferentiellede Monge- Am p2re, dont l’integration se fait par la methode des caracteristiques sous forms finie 
et sans quadrature, non seulement pour le cas de transformation isotherme mais aussi pour celui de trans- 
formation adiabatique. Le phenomene de Riemann-Hugoniottrouve & cet egard sa representation 
geometrique par une arete de rebroussement dans la surface du deplacement. 

Teopua IIOCKEX TA30BEIX BOAIH C KOHeYHOK INMpoToi KoleÖaHuf SHAYHTEIBHO YIIPO- 
INAETCH, EeCJIU IIPUHATB B KAyecTBe HE3ABMUCHMOHÄ IIepeMeHHOH IMPONONBHOE MepeMelleHme. 
loesennee ynoBieTsopaer auhhepennnansHomy ypasHenmo .MOHSK-AMIIEPA, uererpas 
KOTOPOTO MO:KeT ÖBITb HAUMeH B 3AKPBITOM BUNE KAk MIA USOTEPMHYECKOTO, TAK U IA alua- 
ÖariyecKoro H3MeHeHHA CocToaHnus. Auıenne PUMAHH-XYTOHNOT’a uaxonur mpu 
3TOM CBOE TEOMETPHYECKOe BEIPA)KEHHEe B NOABJICHHN BOSBPATHOTO pe6pa Ha HOBEPXHOCTH 
HUHTerpana mepeMeıeHngn. 


Die Theorie der Gaswellen von endlicher Schwingungsweite verdankt ihre ersten grundlegen- 
den ErgebnisseB.RiemannundH.Hugoniot!). Zu dieser Theorie will die folgende Arbeit 
einen Beitrag liefern. Sie will zeigen, wie sich unter Anlehnung an klassische Methoden der Inte- 
grationstheorie der Monge-A mp ereschen Differentialgleichung und durch Einführung der 
Begriffe der ‚‚longitudinalen Verschiebung‘ und ‚Dehnung‘ eine besonders durchsichtige 
Darstellung dieser Theorie gewinnen läßt?). Nach Aufstellung der Monge-Ampere- 
schen Differentialgleichung für die Verschiebung wird zunächst eine Theorie der kleinen 
Schwingungen (Schallwellen) entwickelt, welche dadurch charakterisiert ist, daß sich die 
elastischen Eigenschaften des Gases durch das Hookesche Gesetz (Proportionalität zwischen 
Druck und Dehnung) beschreiben lassen. Anschließend wird die Integration der die Ausbreitung 
ebener Wellen von endlicher Schwingungsweite beschreibenden Monge-A mp £re schen Diife- 
rentialgleichung unter Verwendung der Charakteristikentheorie erledigt, zuerst für den Fall 
isothermisch-isentropischer Zustandsänderung, dann für den adiabatisch- 
isentropischer, wobeinur stetige Vorgänge betrachtet werden. In beiden Fällen 


“ ergibt sich die Lösung, zu deren Gewinnung nicht einmal eine Quadratur erforderlich ist, in ge- 


schlossener Form. Die Eigenschaften der Integralflächen der Verschiebung, besonders die Rie- 
mann-Hugoniotsche Erscheinung, die ihren geometrischen Ausdruck in dem Auf- 
treten von Rückkehrkanten der Integralfläche findet, werden beschrieben und an Hand 
von Beispielen erläutert. 


81. 
Aufstellung der Differentialgleichung. 


1. Verschiebung. Denkt man sich kartesische Keordinaten x, y, z eingeführt, die 
positive x-Achse in Richtung der Ausbreitung der Wellen, dann ist die Lage jeder der y, z-Ebene 
parallelen Ebene durch die Lagrangesche Koordinate a charakterisiert. Danach ist im an- 
fänglichen Gleichgewichtszustand des ruhenden Gases die Lage der in einer.Ebene senkrecht zür 
Fortpflanzungsrichtung liegenden Gasteilchen durch a bestimmt. Bei der Ausbreitung einer 
ebenen Störung erfährt die Ebene «eine longitudinale Verschiebung ss, die von der 


1)B. Riemann, Abh. d. Kgl. Ges. d. Wiss. Göttingen 8 (1860), Ges. Werke 2. Aufl. (1892), 8.43; H.Hugo- 
niot, Journ. €c. polyt. 33 (1887), S. 477. 
-2) Vgl. hierzu die Darstellungen von J. Ha damard, Lecons sur la propagation des ondes, Paris 1903, Kap. IV 


8.143; E. Goursat, Cours d’analyse, Bd. III, 5. Aufl., Paris 1942, S. 127, Fußnotel; K. Bechert, Ann. Physik 


(5) 37 (1940), S. 89, 38 (1940),.S.1; H. Pfriem, Forsch.-Ber.12 (1941), S. 51. 


sie ist Ihlensaihs wicder eine stetige Funktion von a und £ ebenso ie 
een Aarau a RE ee 

e # j = “ir \ De BR . 3 e: 
en 'e Frr 5 ; ng a een 


r z 2. erde Diiatation: 45 Deiden ee ee 


eine Scheibe von der Dicke h. Beider Ausbreitung einer Bahre, möge die a = aa Hr h> 
die Scheibe < x, x > übergehen, wobei Be : 


> ’ =a+s(a,t), EEE Ban 
Als De Ihn ungu wird dann der folgende Grenzwert eg: RT Er: 
un aeg ER a ir Blase 1) — 5 (a, t) Bi See % EN 
RE r g; 2 x ka sus - ERSRNT. u ee er RR da‘ ? KR: ’ 


ei Dehnung ist somit die Ableitung. der. Verschiebung s nach der Kooninare a. Er die 
“Dehnung % positiv, soliegt Verdünnungder Gasteilchen in der Ebene a vor, im ne x 
= Fall ‚Verdichtung. Ins 2 = 


...3.Bewegungsgleic hun 8. Am einfachsten gelangt man zu der ee. 

“ wenn man zunächst die Eulersche Koordinate x verwendet. Ist der Gasdruck 9 (z,t) eine 

differentiierbare Funktion und die Gasdichte 0 (x, t) stetig, so hat man, wenn keine äußeren Kräfte 
vorhanden sind, als Bewegungsgleichung i in Eulerschen Koordinaten: fee 


NER £ 
Li z 


x | d een. 
3 | | - o(®, i) A = 5. I : ee 
Führt man La grangesche Koordinaten ein, so folgt.wegen = a + s(a, t) ; Bo A 
: da 9x da ee ; | 

damit ergibt sich als Bewegungsgleichung i in Lagrangeschen Koordinaten 5 
- \ | Os 0.0 = SER 
ne ie E 


= Br 4.,Kontinuitätsgleichung. Die Kontinuisätsglerchung ergibt- sich aus der An- 
nahme, daß die Gasmasse bei der erhalten bleibt. Da die Volumina umgekehrt Propor- 


V € 
tional den Dichten sind, also Te 7 gilt, wobei der Index 0 die betreffenden Größen im Gleich- 
3 
gewichtszustand charakterisiert, und NS 0 2 = Fee 0 nn —=1-+ wist, folgt 
Ta 0 N Ba 
Er 
Mit (1.2) vereinfacht sich die Bewegungsgleichung (1.1) in | 
0°8 0» i 
A RIFF 73 ER N Tr in, © 2% (1,5). 
5. Zusammenhang zwischen Druck und Dehnung. Zur vollständigen / 


Beschreibung des Vorganges wird noch eine Beziehung zwischen Druck und Dehnung benötigt. 
Setzt,man voraus, daß die Temperatur jedes Gasteilchens zu Beginn des Vorganges dieselbe ist 


®) Bei der Betrachtung der Vorgänge in einem Gas, das mit der konstanten Geschwindigkeit ®,==0 in Richtung 
der. x-Achse strömt, hat man 
z=a+9,t + 8(a,t) | 
zu setzen. Da die Geschwindigkeit v, ohne Einfluß auf die im folgenden für s formulierten Randwertaufgaben ist, wird 
der Einfachheit halber immer 9», = 0 angenommen. 


ee 


‚das Verhältnis der spezifischen Wärmen bedeutet. 


6 Die Ditferentialgleichung der Verschiebung. reizt man 
\ u durch | a und differentüert ‚nach .Q, So Folgt, | 


Ca ODE 
Os\r ti da: 
Da 


FE 

man beider Mon ge-Ampere schien Die naneeene angelangt. welche die ee z 

| . vorgänge in Gasen bei isothermisch-isentropischer oder adiabatisch-isentropischer Zustandsände- = 
+ rung beschreibt. 


7.Kleine Schwingungen. Nimmt man an, daß die Dehnungen und damit auch die. 
Druckschwankungen klein sind, dann kann man die rechte Seite von (1.4) durch die beiden ersten 
Glieder der Eobwioklung nach u ersetzen und erhält so a 


ı . ’ € a 


FR - | = 

= Er ee 
ao . y Pe RE = 

RE Dies Kormet geometrisch Bau hinaus, die Kurve y= er durch ihre Tangente an der Stelle 


u 0 zu ersetzen. Da der Druck p nur - positiver Werte fähig ist, und p=f(uw) nach (1.4) nur ” 


für u Se erklärt ist, gilt (1.6) nur im Intervall —l<wu< = 


"Der durch a. 6) gegebene Zusammenhang zwischen Druck er Dehnung stellt die Zustands- 
gleichung eines Gases dar, dessen elastisches Verhalten durch das H o o ke sche Gesetz beschrieben 
wird, denn es ist nach (1.6) 

08 
Bam Traun u 
el E=xp % Elastizitätsmodul des Gases ist. Une: dieser Annahme erhält 
“man an Stelle der Differentialgleichung (1.5) eine Differentialgleichung mit kon-, 
stanten Koeffizienten vom hyperbolischen Typus* 
0°8 0% 8 % Do BE, 
„ — =0? — , a N: BE 
| © = (1.7) on 


Bu Sn Kal Lund na Aa BR N nn aa I Se 
x R I [ \ 
f 


Aa 
x 


PER TIATN 


i Zum Vergleich mögen noch die Differentialgleichungen angeschrieben werden, welche die 
g Ausbreitung einer Störung unter Benutzung der Eulerschen Koordinaten bei Annahme der 
ä Gültigkeit des Hookeschen Gesetzes beschreiben. Man gelangt dabei auf das folgende System 
4 partieller Differentialgleichungen für Druck und Geschwindigkeit: 


r 
\ 


(2,4 2) Er 23 
0% 9/x m + Mm —P 0% 
0V Be, 0Pp: 
x re a De hr Peer 
< Dieses System ist nicht linear und augenscheinlich von verwickelterer Form als die Differential- 


gleichung d. 7), die sich bei ans Lagrangescher Koordinaten ergibt. 


Fern en hier urch die nahme 
die Din a 0) führte, Die letztere ‚geht ge) Sr 


g a hie Schienen zu studieren, möge a die Lösung er Rand 

verden, welche die Schwingungen in einem durch die Ebenen «= 0 und a= 
ae ‚Im aa a) besitze. jedes Gasteilchen 0<a Bag eine Anfangs 
ne, | j s(a, 0) = Bl Be 
nt Ani (,0)=cy(a),esseiako ° N. gr St = 
= a en 2 (0) plain are 
ee) ara Y (a) seien im Intervall < 0, 1> erklärte Funktionen; o (a) muß, ‚das (a, %) stetige zweite & 
Ahleitungen besitzen soll, mit einer stetigen zweiten Ableitung versehen sein, y (a) aus demselben 


u 


a - Grunde mit einer stetigen ersten Ableitung. Da die Ebenen a=0unda=! dauernd i in Ruhe 5 


: ‚bleiben, muß. 


sein. uferdem sei 


: | ‚ er » ar e = h 1 m £ . r. k 2 | pe; \% 
= | 5 Aula ss 2 a 
> = vorausgesetzt. Im übrigen seien — Funktionen p (a) und (a) beliebig nz In den beiden. N 
NEN: Grenzebenen 4a= O unda=/seifürr>0 pe 


Ei f: 2(0, 2) 8l, a ER ER a8). 


Die eben formulierte Randwertaufgabe und deren Lösungist von der beiderseits. ten 


schwingenden Saite her bekannt. Während es sich aber bei der schwingenden Saite um 
transversale Schwiungngen handelt, liegen hier longitudinale Schwingungen vor, 
denn die Verschiebung s erfolgt hier in Richtung der @a-Achse . Das hat zur Folge, daß den Funk- 
tionen 9 (a) und y (a) gewisse Einschränkungen auferlegt werden müssen, welche ihren Ursprung 


in den geometrisch-kinematischen Eigenheiten der Tone iger Schwingungen haben. Im Zeit- 


“punkt ?= 0 ist die Eulersche Koordinate 


(lea Ht ala), Se RT AR 
- Durch (2.7) ist jedem a aus < 0,2 > ein x zugeordnet. Da u nach (1.6) auf das Intervall 
b<% er beschränkt ist und ı (a, 0) = Y'(a)ist, muß gelten / 


f 1 ' 
Wegen der Voraussetzung (2.5) ist diese Bedingung immer erfüllt. 2e= «+9 (a ) ist daher eine mit 


amonoton zunehmende Funktion und die durch (2. 7) erklärte Ankalduna des Intervalls 
<0,2> auf sich selbst ist umkehrbar eindeutig. 


Gibt man sich eine den eben ausgesprochenen Voraussetzungen genügende Anfangsverschie- 
bung 9 (a) vor, so hat man als anfängliche Druckverteilung nach (1.6) 
p(a,)=yla)=pmll—xp (a). —° 
Dabei hat man als Konstante p, diejenige positive Zahl zu nehmen, welche dem konstanten Gas- 


druck entspricht, "der vor der Verschiebung, also im Gleichgewichtszustand des Gases, vorhanden 
war. Gibt man sich umgekehrt eine anfängliche Druckverteilung x (a) vor, dann ist nach (1.6) _ 


Po —x(a) 
Po 


Dabeiist x (a) gewissen Einschränkungen zu unterwerfen. Zunächst muß x (a), dag (a) eine stetige 
zweite Ableitung besitzen soll, eine mit stetiger erster Ableitung versehene F auktion sein. Aus 
(2.8) ergibt sich durch ‚Integration und wegen pP Ve =V 


1 
En. 
(%) 


y' (a) = (2.8). 


ge A 


a EA DR aaa 


N VER DE 


. bung (a) folgt dann durch Integration aus (2.8) mit 9(0)=0. 


5 erhält man als Lösung der gestellten Randwertaufgabe): 


Er, rade und periodische Funktionen mit der Periode 21 fortgesetzt zu den- 
ken. Die Lösung zeigt, daß s (a, x) in bezug’auf 7 die Periode 2 besitzt, also in bezug auft die 


Periode 2. Naps I Be 


Auf Grund der Voraussetzung (2.5) folgt aus (2.11) al< . Infolgedessen ist die für jedes feste 
7 Vrducch er + s (a, 7) vermittelte Abbildung des Intervalls <O ‚2) auf sich selbst umkehrbar 


und die Dichte 


1 
= 
Ä 


Re 
De N „” 


* 


ügen. Zur Bestimmung einer geeigneten Anfangsverschiebung @ (a) kann man sich demn 
enügende anfängliche Druckverteilung x (a) vorgeben, die entsprechende Anfangsver 


Im Anschluß an die klassische Behandlung des Problems der s c hwi n ge nde n Sa 


art 


2a d)=zlpiatntpla—n] +2 [ve a 


E\ 


0 
2 . 


Dabei sind die Funktionen pund y über das Definitionsintervall <0, 1> hinaus als ung e % 


a-—r 


Aus dieser Lösung lassen sich alle interessierenden Größen ermitteln. Es ist die Eule r- 
sche Koordinate VEN, PERS ARE GEN EN; BERN. 
| 5 REN i x(a,r)=a+s(a,r). 

Die Dehnung ist gegeben durch _ 


ER 1 
u(a, 7) ESTER: 


ar “ i Bi 


Part + Warp]. . al. 


eindeutig, es tritt also bei der Verschiebung der Gasteilchen keine Vertauschung der ursprüng- 
lichen Reihenfolge der Ebenen ein, denn aus a, > a, folgt x (a, T)> x (a,,T). 
Die Geschwindigkeit ist 2% 


v(a, T) = cs. (a, N)= 5 [p’ («+ r) =:E (a—r)]+ = [(y(a-+ 7) 2: v(a—r)] 2128 h 


Ihr Betrag liegt wegen (2.5) immer unter == Die Geschwindigkeit des einzel- 


nen Gasteilchens bleibt demnach fürx)Jlimmer unter der Schall- 
geschwindigkeit .c. 
Der Druck ist nach (1.6) 
p(a,1)= m [1— rule, nm], 


EIIRTONL EN. 54 ee, 
1+u(a,r) Tg 


Alle diese Größen sind periodische Funktionen von z mit der Periode 21. 


o(a,T) = 


Wegen uk läßt sich aus = a + s (a, r) immer a als eindeutige Funktion von x und r = 


gewinnen, die stetige zweite Ableitungen nach x und r besitzt. Somit sind u, v, p und o eindeutige, 
mit stetigen ersten Ableitungen versehene Funktionen der Eulerschen Koordinate x und der 
Zeitt. Willman, um baispielsweise das Verhalten der Geschwindigkeit zu studieren, die Fläche 
darstellen, welche das Bild von v (z, r) ist, so bedient man sich zweckmäßig der folgenden Para- 


- meterdarstellung dieser Fläche: 


z=a+s(a,r), 
1 0sasl, EV. 
Be a Te 


Entsprechend geht man bei der Darstellung der Flächen von Dehnung, Druck usw. vor. 


“) vgl. z.B. L.Bieberbach, Differentialgleichungen, 3. Aufl., Berlin 1930, S. 363; J.Horn, Partielle Diffe- 
rentialgleichungen, 2. Aufl., Berlin u. Leipzig 1929, S.9; E.Kamke, Differentialgleichungen reeller Funktionen, Leipzig 


1930, S. 416, 


(0, &) ER 


N Es a dasjenige nieht! identisch verschwindende Integral : s (a, 7) von D. 1 ee 
‚Randbedingung (3.1) erfüllt und längs einer wegen s (0,0) — 0 durch den Anfangspunkt gehende 
und im ersten Quadranten der a, r-Ebene verlaufenden Kurve (C) mitsamt seinen ersten und 
zweiten Ableitungen verschwindet, welches also längs (©) mit Ber Integral s=0eine Be r ü h 
. rung zweiter Ordnungeingeht. RE TER 

a2 Das’letztere wird gefordert, um in "Übereinstimmung mit der Definition der Lösung einer 
EG: partiellen Differentialgleichung zweiter Ordnung zu bleiben, nach welcher nur solche Funktionen 

' als Lösung zugelassen sind, die stetige zweite Ableitungen besitzen 5). s(a, r) ist dann im ganzen 
" ‚ersten Quadranten der a, r-Ebene ein Integral von (2.1). Die unterhalb (0) gelegenen Punkte ent- 
sprechen den zur Zeit noch in Ruhe befindlichen Gasteilchen . a. Da die Kurve (O8 gleich- 
zeitig zwei verschiedenen‘ Integralflächen angehört, muß sie- 
eine Charakteristik der Differentialgleichung (2.l) sein. Als durch den 

- Anfangspunkt gehende im 1. Quadranten verlaufende Charakteristik von (2.1) ergibt sich a =r. 
Daraus folgt: Die ebene Welle breitet sich mit der Schallgeschwin- i 
dıokeitc-in Richtung der positiven a-Achse aus. 

Es bleibt noch übrig, ein nicht identisch verschwindendes Integral von (2.1) zu bestimmen, 
welches die eben ausgesprochene Randwertaufgabe löst. Hierzu sind die folgenden Voraussetzun- 
gen über die Funktion :o (r) zu machen. "Da die zweiten Ableitungen des Integrals stetig sein 
sollen, muß 0” (r) ebenfalls stetig sein und wegen s (0,0) = 0 muß o (0) = 0 sein. Außerdem sei 


ee ee ee 


und | 
1 % e 

° z<em<i ee A ee a 

vorausgesetzt. | "T 
Zur. Lösung der Randwertaufgabe. geht man von dem allgemeinen Integral von (2. 1) aus: 


s(a)=flrtatge—a). | h 
Da s auf‘der Charakteristik r— a verschwinden muß, folgt f(2a) +g(0)= 0, und es bleibt SER 


s(a,r) =g(t—a) —g (0). Mit der Anfangsbedingung (3.1) ergibt sich dann s (0,7) = o(r) = 
g(T) —g (0) und damit als Lösung der Randwertaufgabe 


s(a,rT)=o(T—.a) ee re er A SEE 
Aus der Voraussetzung o (0) = 0 folgt s (0,0)—= 0. Den aus (3.4) old Gleichungen 
ee Ar 5 = 0o'(T—.a) | 

oa BER, OR e N 

025 028 Br Be; 

Ne re) ; 
entnimmt man, daß längs der Charakteristik a = r wegen (3.2) die ersten und die zweiten Ab- 
leitungen verschwinden. Da a = —o’(r —a) ist, folgt auf Grund der Voraussetzung (3.3), daß , 


u die Bedingung —1<ü< = erfüllt. e 


Als Beispiel für die Auchreifung einer Welle kleiner Schwingungsweite möge die Lösung der 
eben behandelten Randwertaufgabe bestimmt werden, die sich für 


N Te Au SO Tee 
90 ® >, 
;) Vgl. E. Kamke, Differentialgleichungen reeller Funktionen, Leipzig 1930, 8. 377. 


A=0,2 
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Fa ee ln a 3 mn 
. =" € ER 
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a rs de al ie ca“ 
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es dsin (T—a)eos (ea). u 


ar 


A u re 
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In Bild 1 ist, die Fläche von. v (x, 7) skizziert*). Die Schnitte 7=konst. der Fläche v (« 


i gegenüber der Kurve — v(0,7)= —3 Asin? TcosT eine Verformung der Wellenlinie an, alı lie 


a Bi; 


Tr - 


= z Te 
Bild 1. 


der, die bei den unten behandelten Wellen mit endlicher Schwingungsweite zu beobachten ist | 


(Steilerwerden der Wellenfront). Insofern liefert die hier dargestellte Theorie der kleinen Schwin- 
gungen eine bessere Annäheıung an die in den beiden folgenden Paragraphen gegebene strenge 


Lösung als die mit Eulerschen Koordinaten operierende, durch welche die eben beschriebene 


Erscheinung nicht erfaßt wird. 
SA. j 
Ausbreitung einer ebenen Welle bei isothermisch-isentropischer Zustandsänderung. 
Das im vorigen Paragraphen unter der Annahme des Hookeschen Gesetzes behandelte 


-Problem der Ausbreitung einer ebenen Welle soll im folgenden unter der . 


Annahme isothermer Zustandsänderungen des Gases gelöst werden. Der Zusam- 
‚menhang zwischen Druck und Dehnung ist dann nach (1.4) durch 


p u 1 > 
Baisin le Erie: eintei an leugier vie. ol a ne Mlaucers ie A 
E e >. IHU J Be 
gegeben. Als Differentialgleichung für die Verschiebung s erhält man mit «= 1 aus (1.5) 
028 2 028 : 


Du = 2? 
ot (145%) 0a 


da 


*)In Übereinstimmung mit derim folgenden gewählten Bezeichnungsweise ist hier die r-Achse als „-Achse bezeichnet. 


. angew. Math. Mech. 
222 Quade, Zur Theorie der ebenen stetigen Gaswellen na.2s/27 Nr.B0 Nov./Dez. 1947 
oder unter Einführung von ct=7T 
| EIN RN VERS NA 
dr ( f) u da? Eis 
Ta 
a 9a 


Auf dieselbe Gleichung kommt man, wenn man eine adiabatisch e Z us tands- 
änderung mit nicht zu großen Dehnungen voraussetzt. Bei der letzteren gilt nach (1.4) 


BER 1 
Bor (tu 
Ersetzt man auf der rechten Seite (1 + w)* durch 1+xu, so hat man 
BU IB Er ee 
Po 1+ru r 
Die Hyperbel y = = besitzt an der Stelle v—= 0 gleiche Ordinate und gleiche Tangente 
wie die Kurve y= u. und stellt demnach bei nicht zu großen Werten von u eine gute An- 
U * 
näherung an die letztere dar. Mit (4.3) ergibt sich an Stelle von (4.2) die Differentialgleichung 
a 1 025 
(ge . dar” 
Ihre 
Setzt man noch x3 = 8, so folgt für $ die Differentialgleichung 
ST 1 028 
Omas EEE 
143) 


also wieder die Differentialgleichung (4.2). Die Behandlung des isothermen Falles vor dem adia- 
batischen ist dadurch gerechtfertigt, daß die Formeln im ersten Fall einfacher sind als im letzten, 
auch läßt sich schon amerstendie Riemann-HugoniotscheErscheinung studieren. 

Die Differentialgleichung (4.2) ist eine Monge-AmperescheDifferentialglei- 
chung. Um die in der Theorie dieser Gleichung üblichen abkürzenden Bezeichnungen verwenden 
zu können, sei im folgenden 


a=r, ER s=2, bzwW2S 2 
gesetzt. Mit den üblichen Abkürzungen 
02 02 022 022 022 
Berne im Na ar aeg ee 


schreibt sich dann (4.2) 
r— (IH P1=0. 2.10, ee 


und es ist die Dehnung « = p und die Geschwindigkeit v = cq. 

Um den Vorgang einer in der positiven x-Richtung sich ausbreitenden Welie zu beschreiben®), 
hat man ein im ersten Quadranten der x, y-Ebene nicht identisch verschwindendes Integral 
z=F(x, y) von (4.4) zu bestimmen, welches längs einer durch den Anfangspunkt gehenden und 
im ersten Quadranten der x, y-Ebene verlaufenden Kurve (C') mitsamt seinen ersten und zweiten 
Ableitungen verschwindet, welches also längs (©) mit dem Integral z = 0 eine Berührung zweiter 

„Ordnung eingeht. Weiter wird von der Integralfläche z—= F (x, y) verlangt, daß sie durch eine in 
der y, z-Ebene gelegene vorgegebene Kurve (I’) hindurchgehen soll. Die Kurve (7) ist dabei durch 
die Bewegung des Kolbens, also der linken Grenzebene x = 0 gegeben. (I') sei definiert durch 


RUN y—n z=@ln).. we Wlan 


wobei (n) die Auslenkung des Kolbens als Funktion von n= 7 bedeutet. Die Funktion & (n) 
erfülle die folgenden Voraussetzungen. Da die zweiten Ableitungen von F(x, y) stetig sein sollen, 
muß w’ (n) ebenfalls stetig sein, und wegen F (0,0) = 0 muß o (0) = 0 sein. Außerdem sei 


o' (0)= w’ (0) = 0. ee 


Zur Lösung der eben formulierten Randwertaufgabe wird im folgenden die Methode der 
Charakteristiken benutzt?). 


°) Vgl. die Formulierung der gleichen Randwertaufgabe in 3. 


’) Vgl. die Darstellung dieser Methode bei E. Kamke, Differentialgleichungen reeller Funktionen. Leipzig 1930, 


e: ei 


a e 
= & P>— Kr BI AR 
as s System v von n gewöhnlichen Ditterentielgleichungen Be 


Be Be ER = Fir +, A FL} 
EI u ER = on SR, 
ne & —g h 0 DS } = 


- Er AN 
ofinient, wenn er nur Punkte von @ enthält und x’ ad y nicht gleichzeitig a, 
die Kurve(Ö)gleichzeitig zwei verschiedenen Integralflächenan- 
gehört, muß sie eine Charakteristik sein. Da längs (O))a=p=g=0gil 

: ey es dem Er (4.7) (€) mit-einer der beiden durch den Aalaaepugl! gehenden, a 


een und y=—% zen 
nmenBllen Yon teen beiden Geraden kommt nur die erste in Frage, da nur sie im En 
- Quadranten verläuft. Die Kurve (OJist also die Winkelhalbierende des ersten Quadranten. Damit 
‚hat man als erstes Ergebnis: Dieebene Welle pflanzt sich mit der Schall- 
‘ges chwindigkeit e in Richtung.der positiven.z- Achseforul ou 
Be : Es bleibt noch ee ein. nicht identisch verschwindendes Integral z=F er Yy) für das 
E07.  Gebiet- A £ 
ar | Bi 
Ba: 2: ; ee vor 


aufzusuchen, denn unterhalb der Winkelhalbierenden y= xist z=0. Ein solches läßt sich hier 

angeben, da sich aus dem System (4.7) ein Vorintegral?°) von (4.4) gewinnen läßt. Nimmt 
man in (4.7) als Parameter die we Veränderliche x, so lauten die beiden letzten Glei- 

chungen von 1ie 7) 


UHR eg 
ER et man peide Seiten der er Gleichung und ersetzt y’? durch (1 + p)2, so hat man 
e-; | P®= (14 p%d*. 5 


Durch Aufspalten dieser oo... ergeben sich die beiden folgenden: 
—1+pd; PFlrDg, 
von denen jede ein a von (4.4) liefert, 
p=(,e!—l, p= 0, —1: 


Von diesen beiden Vorintegralen kommt hier nur das erste in Betracht. Denn von den beiden 
‘Scharen von Charakteristiken, welche man durch Aufspalten der mittleren Gleichung von (4.7) 


REN REN Re LE FOREN, arg N 
\ Ye, 2 Je 
s„ N ’ 


erhält: 
a 2 =pRtgy), 
y+{l+tpnz=0, ud 0 Eu a DE 
65 . 2 = pe —dy=d, ve" —gdy—=d, 


. muß zum Aufbau der Integralfläche diejenige Schar ausgewählt werden, welche die Charakte- 

“ristik y= x nicht enthält, und das ist gerade die aus dem ersten System sich ergebende Schar, 
zu welcher das erste der Vorintegrale gehört. 

- Esi u =yp>I 4 =M= dein auf der Winkelhalbierenden gelegenes Flächen- 

element. Dann ist auf dem dieses Flächenelement enthaltenden charakteristischen Streifen erster 

Ordnung wegen 9, = 9% = 0 die Konstante C,=1 zu setzen, und es besteht daher auf jedem 


solchen Streifen die Beziehung 


az 2 TE W u DIR, ET 


p=en—1 Se ee RE EN 


Damit ist die Aufgabe auf die Integration der partiellen D tere ichung 
erster Ordnung (4.9) zurückgeführt. Jede Lösung der letzteren ist eine Lösung von (4.4). 


°) Vgl. E. Kamke,a.a. O., S. 394. 
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® Auch diese Differentialgleichung erster Ordnung möge unter Verwendung der Charakteristiken- 
mzsthode integriert werden. Die Charakteristiken dieser Gleichung ergeben sich aus dem 
System ?): | Bi 
Y (a),=ent, 2 (2)=erı—1+get, «()=0.... . (410). 
Ein auf der Kurve (T') liegendes und der Integralfläche z = F (x, y) angehörendes Flächen- 
element ist nach (4.5) und (4.9) : 
=), yan, z=0o(n®» = ee mM—1, go um r1 
Aus der letzten der Charakteristikengleichungen (4.10) folgt FR 
q=konst.=w’(n), 
und damit lauten die beiden ersten dieser Gleichungen 


em "(= [1 Folie 1 Trier 
Die Integration dieser Gleichungen unter Berücksichtigung von (4.5) liefert schließlich 
yzxzeoM+n, :=[(1+o’m))e®w—1]z+to(n) . ..:. (4.13). 
Damit ist die gesuchte Integralfläche z— F (x, y) in der Parameterdarstellung 
y-y(un),.. - 2=y(lon), 220, nz=0 


gegeben. Wie (4.13) zu entnehmen ist, sind die Charakteristiken in der Integralfläche liegende 
Geraden, d.h. die Integralfläche ist eine Regelfläche mit der Leitkurve (JI') 1). 


Es bleibt zu zeigen, daß die gefundene Integralfläche längs der Winkelhalbierenden y = x 
eine Berührung zweiter Ordnungmitder Ebene z—= 0 eingeht. Nach der ersten der 
Gleichungen (4.13) kann nur diejenige der Charakteristiken (4.13) durch den Anfangspunkt gehen, 
bei welcher 7 = Oist. Da auf Grund der Voraussetzung (4.6) »’ (0) = 0 ist, fällt diese Charakteri- 
stik mit der Winkelhalbierenden zusammen. Dann verschwinden aber wegen qg= »’(0) und 
p= e-®’(W —1 die ersten Ableitungen p und g längs der Winkelhalbierenden; es findet also 
längs der Winkelhalbierenden eine Berührung erster Ordnung zwischen der Integralfläche und der 
Ebene z= 0 statt. Daß auch die Ableitungen zweiter Ordnung längs der Winkelhalbierenden 
übereinstimmen, sieht man so ein. Aus 


de=pda+gdy=|[e ""—1]ldax+o(n)dy 

folgt 

een year (4.14). 

Nach (4.13) ist 
| dy=e "Mdax-+[1 20" m)" Man. 

Unter der Voraussetzung A= 1 —xw’’ e7®’ +0 ergibt sich aus der letzten Gleichung 


1 


dn=(dy—e " de), 


und dies in (4.14) eingesetzt liefert 


rn 


[0] ’ o’ 


P2= (er dr —2e " daedy+dy) 
Daraus liest man ab 
a ... (4.15) 
A A ? 


Da nach der Voraussetzung (4.6) »’’ (0) = 0 ist, folgt aus (4.15), daß r, s und t für n=0(0,d.h. 
längs der Winkelhalbierenden, verschwinden. Die durch (4.13) dargestellte Inte- 
gralfläche geht daher mit der Ebene 2z=0 eine Berührung zweiter 
Ordnung ein. £ 

An Hand der eben angeschriebenen Werte von r, s und t bestätigt man übrigens unmittelbar 
daß die durch (4.13) dargestellte Integralfläche der Differentialgleichung (4.4) genügt, und daß 
überdies r? — s?—= 0 ist, die Integralfläche also eine abwickelbare Re gelfläche ist. 

Bei der Bestimmung der zweiten Ableitungen wurde A +0 vorausgesetzt. Ist A(&, n) in 
der Umgebung eines Punktes (x,,n,) von Null verschieden, so läßt sich die erste der Gleichungen 


®) Vgl. E. Kamke, a.a. O., 8. 348. 


1°) Vgl.z.B.W. Blaschke, Differentialgeometrie, Bd.1, S. 72 u. 109, oder V undK.K i 
der Raumkurven und krummen Flächen, Bd.1, Berlin u. Leipzig 1931, S. 50.  Reler 


| ler Integral { rve, | 
a ist. Für die Rückkehrkante “ ergibt si sich a aus  A= ar 
en ‚die folgende Dr oe, 


RE ET SH 2 wen 
ee : v-; er = In 


Sr = ar 


® win man den zeitlichen Naluk ar Bean einer Ebene : u an = 2, im Fe 
verfolgen, so hat man die seen 2 dieser Ebene als Funktion von yzu untersten; ebenso Fe 


Ber Parameterdarstellungen gegeben: 5 


met 3 a4) Ham. a 


ER 1— 2,0’ I 


"Schließlich sei noch die Darstellung der Geschwindigkeit i in Eu e sen Koodmales ge 
 goben Bedeutet X=x%+ 2die Eulersche Koordinate, so ist die GEERUnGES v als Funk- 
tion von = zul Y durch die folgende Parameterdarstellung gegeben: 


Z=-wnlite)to,  I=0. 


U REN ie ndsı, ir, ET 


IERTEN 


Bild 2. 


Ein Böispiel möge das geometrische Bild einer solchen Fläche vor Augen führen. In Bild 2 
Fe ist die durch (4.13) erklärte Fläche dargestellt, die sich für 


Asin®n für O<n<a, ! 
o(n)= { mo, ,. 4-02 
15 
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ergibt. In die Figur sind as Schnitte y= konst. = y, der Fläche sowie die Charakteristik „=x 
eingezeichnet. Für die Schnittkurven erhält man die Parameterdarstellung 


= (ya), 
; EN TORE RR Te es 4.17 
= (9 NN oO er SL, 
In Bild 3 ist die Projektion der Rückkehrkante in die x, y-Ebene dargestellt. Von dieser 
Kurve kommen hier nur diejenigen Zweige in Betracht, die im ersten Quadranten der x, y-Ebene 
liegen. Die Fläche besitzt hier erst für y > 3,12 eine Rückkehrkante. Die Fläche g=qg(X,y) 
ist in Bild 4 skizziert. Für verschiedene Werte von y= konst. — y, sind die Schnittkurven 


Bild 3. 


% ıge | 

. teilerwerden ee Welle“, tpu ® dab. 
Finder Integralfläche | z=F(x,y) die Erden der Rückkehrkante bemerkbar mac 
TER _ Punkte (X, y), in denen Köln eindeutiger Geschwindigkeitszustand vorhanden WERTE 

‚ Die eben beschriebene in dem Auftreten einer Rückkehrkante begründete Eigenart er 

 gralfläche und die damit verbundene Unstetigkeit i in der Beschleunigung von Ebenen = k 
 sowie.die Existenz von Punkten, in denen ein mehrdeutiger Geschwindigkeitszustand h 
wird als Riemann-HugoniotscheErscheinung bezeichnet. Offenbar’ hat di 
 scheinung ihren Grund in den Hypothesen, welche den Gleichungen des Vorgangs zugrunde liegen 
- und insofern der Natur des Vorgangs der Wellenausbreitung nicht voll gerecht: werden KUuRTEE 
‚sie die Einflüsse der Reibung, ‚Wärmeleitung u. a. ‚unberücksichtigt | lassen. a 
NA . Hier kann man fragen, ob es eine für > 0 erklärte, mit stetiger zweiter Ableitung versehene, SE 
den Bedingungen ® (0) = ®’ (0) =» (0) =0 genügende Funktion ® (n) gibt, für welche die 


Be; . zugehörige Integralfläche z=F (x, y) keine Rückkehrkante besitzt. Die Antwort auf die Frage 
KR: lautet: Jede mit einer nicht identisch verschwindenden, deneben 
Fa Br ausgesprochenen a genügenden Funktion o(n) kon- 

2 SSTLULEr te Integralfläche 2=F(&,y) besitzt eine Rückkehrkante. In. 


' der Tat, die zu einem Wert 7 20 ee x-Koordinate der Rückkehrkante sich aus A — 0, u 
also aus 7% 
p: | | zo, ER 2 
E _ Verschwindet o (n), also auch &”’ (n) nicht identisch, dann gibt es wegen der Stetigkeit von o’ein 
ein „-Intervall, in welchem &” von Null verschieden ist, und zu jedem » dieses Intervalls existiert 
auf Grund der Gleichungen (4.16) ein Punkt der Rückkehrkante. Ist dagegen o (n), also auch 
©" (n) identisch Null, so gibt es zu keinem n eine x-Koordinate der Rückkehrkante, also auch 
eine Rückkehrkante, und es liegt das triviale Integral z = 0 vor. 
Danach ist das Auftreten der Rückkehrkante unabhängig von der Ordnung der eh ET 
welche die im Gebiet g:x 20, y 2x erklärte Integralfläche 2=F (x,,y) längs der Winkelhal- 
‘ bierenden mit der.Ebene 2 = 0 eingeht. Nimmt man beispielsweise 


EN 


kai a 


KOREN NEE UROTETRIERENS a a 


[deren für, (<< 
o(n) = 


Ar jr 
DENE vr 


so hat man längs der Winkelhalbierenden eine Berührung beliebig hoher Ordnung, da jede Ableitung 
von @ (n) an den Stellen 0 und » verschwindet. Die mit dieser Funktion & (n). konstruierte Inte- 
gralfläche besitzt, da  (n) nicht identisch verschwindet, eine Rückkehrkante. z 
/ Von Interesse ist weiter die folgende Frage: Gibt es nicht identisch verschwindende Funk- 
tionen & (n), für welche entweder die »-Koordinate oder die y-Koordinate der Punkte der Rück- 
kehrkante nur negativer Werte fähigist. Wenn es solche Funktionen gibt, so tritt bei den diesen 
Funktionen entsprechenden Vorgängen der Wellenausbreitung de Riemann-Hugoniot- we 
sche Erscheinung nicht auf, da die gesuchte Integralfläche nur in dem Gebiet: x 20, y 0 in- en. 
teressiert. Die Antwort auf die Frage lautet: ‘Ist o(n) eine nichtidentisch ver- Be 
schwindende, denoben ausgesprochenen Bedingungen genügende a 
Funktion, welche für „20 der Ungleichung @’<s0. genügt, dann 
besitzt die mit einer solchen Funktion o(n) konstruierte Inte- 
gralfläche 2=F(a,y) eineRückkehrkante, die außerhalb des Ge- 
bietes z>20 y 20.1iegt. 
Dies folgt nach (4.16) unmittelbar aus den Gleichungen 


e@ 1 


Nimmt man beispielsweise die Funktion 


o(nN)= — (1-+n®)+n—arcign, 


deren zweite Ableitung — a ist, für welche also die y-Koordinate der Rückkehrkante 
1% 0 ist, so tritt bei dem zugehörigen Vorgang der Wellenausbreitung die Riemann- 


Hugoniotsche Erscheinung nicht auf. Die.eben angegebene Funktion & (n) beschreibt 
eine Kolbenbewegung, bei welcher o(n) <0 ist. 
- 15* 
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Im Anschluß an das eben erörterte Problem der Ausbreitung ebener Wellen bei vorgegebenem 
Bewegungsgesetz der linken Grenzebene x = 0 läßt sich sofort die Lösung des folgenden angeben. 
Es sei nun nicht mehr das Zeitgesetz für die Bewegung der linken Grenzebene gegeben, sondern 
das Zeitgesetz für den in der linken Grenzebene z=0 wirkenden 
Druck. Bedeutet P (x, y) den Gasdruck in der Ebene x zur „‚Zeit‘ yund P, den Gasdruck im 
Gleichgewichtszustand, so ist nach (4.1) die Dehnung 


Ph P>0, P=0- 


Ist P (0, n) bekannt, so ist es nach der letzten Gleichung auch 7 (0, n). Esseip (0), n)=P(n) > — 
eine für n > 0 mit stetiger Ableitung versehene Funktion; außerdem sei N 

0.= Eee el 
Aus dem Vorintegral (4.9) ergibt sich 


| O0, )=—-ıhil+em)=en. 
Mit & (0) = 0 folgt hieraus 


ni 
omM=—/inll+g(alda. ENTE PRO 


Man bestätigt leicht, daß auf Grund der über 9 (n) gemachten Voraussetzungen »’ (0) = w’’ (0)=0 
IREN=®- 
Damit ist die eben formulierte Aufgabe auf die eingangs behandelte zurückgeführt. Man hat 
nur für ein den Bedingungen (4.18) genügendes p (n) die Funktion ® (n) aus (4.19) durch eine 
Quadratur zu ermitteln. 


S5. 
Ausbreitung einer ebenen Welle bei adiabatischer Zustandsänderung. 


Der Fall der adiabatischen Zustandsänderung zeigt gegenüber dem im vorigen Paragraphen 
betrachteten nur geringfügige Abweichungen. Die Randwertaufgabe und das Integrationsver- 
fahren sind genau die gleichen, lediglich in den Formeln treten Unterschiede auf und über o (n) 
ist noch eine weitere Voraussetzung zu machen. 


An die Stelle von (4.4) tritt de Monge-Ampe&resche Differentialgleichung 
FF FAT HS 0 ee 
deren charakteristischer Streifen erster Ordnung den Gleichungen 


=putgy, 
y2%—(1+p)"tHe?=0, NN ee ED 
p' a’ —q’ y' =) 
genügt, wobei x eine beliebige reelle Zahl > 1 bedeutet. Das Vorintegral lautet hier 
ß 
= ”—1 


1, Mer eo Re 


»=(1+gq 5 


Da die rechte Seite von (5.3) bei beliebigem reellem positivem A nur für > _; erklärt ist, 
können hier nur solche Werte von q zugelassen werden, welche die Zahl — nicht unterschreiten. 
Charakteristiken dieser partiellen Differentialgleichung erster Ordnung folgen aus dem 

ystem 


y()=(1+Ag)”, (ii? Ita), gdia)=0, 
re rt 
2 2 ul 


Aus der Forderung, daß die Integralfläche z = F (x, y) durch die Kurve 


(2) =; y=n,2=o(n), n 20 
gehen soll, ergibt sich 


gq=w (n), pet 
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Wegen qg> ET muß hier die Funktion o (n) außer den im vorigen Paragraphen ausgesprochenen 


Voraussetzungen noch die zusätzliche 


o(M)>— A A (5.5) 


erfüllen. Aus (5.4) folgt dann 
=] Zm)=[1+AR ler (Ml—L, 


und daraus durch Integration 


y=&:[11+4o(n)] ”+n, 2=ef{l1+4o0(n)] ’I+uo(n)]—1}to(n) .. (6.6), 
womit die gesuchte Integralfläche in Parameterdarstellung gegeben ist. Auch hier sind die Cha- 
‚rakteristiken in der Integralfläche liegende Geraden, es liegt also wieder eine Regelfläche 
mit der Leitkurve (I) vor. Unter der Voraussetzung o (0)=o’(0)=0 geht 
diese Integralfläche längs der Winkelhalbierendeny=x eine Be- 
rührung zweiter Ordnung mit der Ebenez=0 ein. Inder Tat, die Winkel- 
halbierende ist die zu dem Werte n = 0 gehörige Charakteristik. .Da »’ (0) = 0 ist, verschwinden 
»p und g auf der Winkelhalbierenden. Der ersten der Gleichungen (5.6) entnimmt man 


dy=(1+Ao) "dc +[1—xuw”(1+%o) ” "]dn. 
Diese Gleichung ist nach dn auflösbar, wenn. 


A=1—uxo” (1 SON 
von Null verschieden ist. Auf demselben Wege wie in $4 ergibt sich 


’ 


7) o’’ 


a ı\ — 2V BER: DV MA sh 
lee) = ) ey 


woraus die Behauptung folgt. Die es, ist hier wieder eine ab w j ckelbare Re- 
gelfläche. 

Inden Punkten, indenenA verschwindet, existieren die zwei- 
ten Ableitungen nicht. Diese Punkte liegen auf der durch 


23 


era). 1440 
rare a +; 
a. FO ee (8.2) 
- 1 
riet, 1+u o —(1+40')”) + o(n) 


u@ 
gegebenen Rüc kkehrkante; es tritt also auch hier im allgemeinen die Riemann- 
Hugoniotsche Erscheinung auf. 

Die Bewegung einer Ebene x = konst. = x, wird durch die folgenden Gleichungen be- 
schrieben, welche die Verschiebung z, die Geschwindigkeit v» und die Beschleunigung b als Funk- 
tionen von y geben. 

y-=n(41+io)”’+n, 2=o,|(1 +40) "(14 uw‘) —1l+o(n), 
v 5 bet. (M) 
c & Also n) 


Bedeutet X=x- z wieder die Eulersche Koordinate, so ist die Geschwindigkeit v als 
Funktion von X und y durch die folgende Parameterdarstellung gegeben: 
® 


X=(y—n)(1+ue))tein) -=-d-6(n): 


C 


In Bild 5 ist die durch (5.6) erklärte Fläche dargestellt, die sich für dasselbe ® (n) wie bei 
dem in $4 behandelten Beispiel ergibt. Für die Schnittkurven dieser Fläche mit Ebenen 
y= konst. = y, erhält man die Parameterdarstellung . 

= (y—n)(1+Aw'), 
= (vn) +Roe— (1+Ao)’])+e(n) 

Bild 6 gibt die Projektion der Rückkehrkante in die x, y-Ebene. Die Fläche v = v (X, y) 

ist in Bild 7 dargestellt. ; 


ew. Mein! Mech 
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Bild 5. 


Was das Auftreten einer Rückkehrkante anbetrifft,.so gilt hier ein entsprechendes Ergebnis 
wie im isothermen Fall: Ist o(n) eine nicht identisch verschwindende, 
den in $4 ausgesprochenen Voraussetzungen und (55) genügende 
Bunktion, dann besitzt die mitieiner soleben Funktion aiy) kon 
Struierte. Integraliläche = Fl(z,y)einerARuckkehrkante, Pas ao 
sich unmittelbar aus A = 0, also aus 

ko’ c=(l1-FAorr! 


und (5.6). Auch hier besitzt nur das aus ® (n)==0 folgende triviale Integral z= 0 keine Rück- 
kehrkante. 


IE 
X 


EIN: 


Bild 7. 


+ Me Antwort auf die u Frage, ob es nicht identisch verschwindende Funktionen © (n) 
gibt, bei denen die Rückkehrkante von z= F (x, y) außerhalb des Gebietes x >0, y > 0 liegt, 


lautet hier: Ist ®(n) eine nicht identisch verschwindende, denin SA 


ausgesprochenen Bedingungen und (55) genügende Funktion, 
welche für 720 der Ungleichung @"s0 - genügt, dann’besitzt die 
nit .esner ‘solchen. Hunktion om konstruierte "Lntegralfiäche 


2=F (x, y) eine Rückkehrkante, die außerhalb des Gebietes 20, yzd 
legt. Dies folgt unmittelbar aus den Gleichungen (5.7) der Rückkehrkante. 


Nimmt man beispielsweise die Funktion 


- für welche die y-Koordinate der Rückkehrkante — —<< 0 ist, so tritt bei dem zugehörigen 


N : 
Vorgang der Wellenausbreitung die Riemann-Hugoniotsche Erscheinung nicht auf. 
Dieser Funktion & (n) entspricht eine Kolbenbewegung, bei welcher &’ <O ist. 
- Es sei noch eine Erscheinung erwähnt, die im isothermen Fall nicht auftreten kann. 
Nimmt man also(n)eine Funktion, deren Ableitung an der Stelle 


Hensden Wert — unterschreitet, bei welcher also die Voraus- 


‚setzung (dö)nicht erfüllt ist, dann zeigt die Integralflächez= F(x,y) 


im Punkte (0, n,) ein singuläres Verhalten, sie ist dort nicht diffe- 
rentiierbar. In der Tat, da die Gleichungen (5.6) die Integralfläche z=F (x, y) nur für 
n<n, erklären, ist die Fläche au gen Geradenstück = 0, y > n, nicht erklärt. Es möge 


daher für y> n, der Grenzwert en of, y) untersucht werden. Der ersten der Gleichungen 


(5.6) entnimmt man 
| e=(y—-n)(lHAo), 0O<n<y, y>ne- 


Danach geht &—-+ 0, wenn „> 9 —0 strebt. Aus der zweiten der Gleichungen (5.6) ergibt 


sich nach Einsetzen des eben angeschriebenen Wertes von & 


=(y—n)(1 ERNT y—n)i+io)P ton). 


ee TATRA | 
, - 5 . « 
. ü Y 
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. Läßt man „> 79 — 0 gehen, so strebt 
2> (y—no)Il+ won) lt® m). 


| 1 
daher gilt wegen &’ (no) = Er 
m Fley=(y—n)(l—r) tom),  y>n- 


Die Grenzkurve ist somit eine Gerade mit dem Steigungsmaß 1 —» durch den Punkt [»,, om}: 
Da die durch 


®(n) : für 0<n<m 
("nl —)to(no) > n>no 


definierte Grenzkurve an der Stelle 7, keine Ableitung besitzt, ist = F (x, y) im Punkte (0, No) 


nicht differentiierbar. Den Gleichungen (5.7) entnimmt man, daß [7,,@(n,)] ein Punkt der Rück- 
kehrkante ist. 


Bild &. 
In Bild 8 ist ein Beispiel für eine solche Fläche dargestellt 11). Dabei ist © (n) = —10 sin? n 
angenommen; es ist dann n, = 0,4440.. 


Für die Lösung der ändwertäufgnhe one hier nur der in dem Gebiet: 2 20,2 <y<n, 
liegende Teil der Integralfläche z= F (x, y) in Betracht, da nur in diesem Gebiet die Stetigkeit 
der Ableitungen von z=F (x, y) gesichert ist. In der linken Grenzebene 2=0 geht für 


Y> no — 0 die Dehnung p— + , der Druck am Kolben somit gegen Null. 
Bei dem am Schluß des $4 erörterten Problem ist hier 


on) = [114 PM) 11a 
| 


zu setzen, wobei wieder @ (0) =’ (0) = 0 und @ (n) > —1 vorauszusetzen ist. 


Abschließend sei noch bemerkt, daß sich für den Fall beliebigen Zusammen- 
hanges zwischen Druck und Dichte 


MR | 0) 
Po ! 00 
wieder eineMonge-Amperesche Differentialgleichung ergibt, deren Integration sich in der- 


selben Weise durchführen läßt, wie dies oben für den isothermen und adiabatischen Fall im ein- 
zelnen dargelegt wurde. 


ıı) Für die Berechnung und Zeichnung der Kurven bin ich Herrn stud. el. M. Würz zu Dank verpflichtet. 
Eingegangen im Oktober 1944, 


_ Einfluß einer schwachen W 

RE der laminaren Grenzschichten. 
| 5 BAER ee ee. 
e= | Von H.Görtler in Freiburg i.B.  — BSR 


Für das im Titel genannte Problem wird eine Theorie der ersten Ordnung für schwache Welligkeiten 
er unter Beschränkung auf kurze Wellenlängen A (relativ zur überströmten Wandlänge L) entwickelt und nm 
en numerischen Beispielen erläutert. Der Grenzfall sehr kleiner Werte von A/L ist mathematisch besonders leicht 
ER zugänglich und führt auf einfache Gesetzmäßigkeiten. | ; h Be 

For ihe problem mentioned. in the above title a first order theory for small undulations with restrichion 
Ei to short wave lengths A (in comparison with the length L of the wall overfloum) is given and numerical 
5. : examples are evaluated. The limiting case of very small values of A/L is of especially easy access 
By ‚analysis and leads to simple relations. | N % FR 

En ce qui concerne le probleme cite dans l’en-tete une theorie du premier ordre est developpe concernant 
des ondulations petites avec restriction aux longeurs courtes d’onde X (relatives ü la longueur du mur inonde 

5  L) et expliquee par des examples numeriques. La limite des valeurs trds petites de AL est ires facilement 
4 abordable au point de vue mathematique et resulte dans des lois simples. - ; 2 


OrpaHnunpasch CIyyaeM HeSHAYUHTEIBHOH BOIHHCTOCTH CTEHOK C KOPOTKOH AUIUHOH BOIHH 
(10 oTHOmeHn® K AIHHe CaMOÄ -CTEHKH) BEIBONNTCH TeOPHA IIePBOTO HOPANKa YRASaHHOH 


| andwelligkeit auf den Verlauf 


Be B 3T0JIOBKe [3ayaym M MINDCTPHPYeTcH UNCHeHHBIM ImpuMmepoMm. B umpenensHoMm caysae ° 
OYeHb MauIEIX 3HayeHuf 2/L MATEMATHYECKHE BEIBOMBI OCOÖEHHO IPOCTH H IIPHBONAT K IPOocToA 
3&KOHOMEPHOCTH. 32 2 
5 ey $ 1. Einleitung 


Eine Untersuchung von J. Pretsch [5]über die Stabilität der laminaren Grenzschichten - 
an welligen Wänden lenkte die Aufmerksamkeit darauf, daß selbst sehr schwache periodische 
Schwankungen der Wandkontur erheblichen Einfluß auf den Verlauf der Grenzschichtströmung 
haben und insbesondere zur Ablösung der Grenzschicht führen können. Dieses für die Frage nach 
dem laminar-turbulenten Umschlag und dem Widerstand der Strömungskörper wesentliche Er- _ 
gebnis wurde angezweifelt, da die Grenzschichtbereehnung nach dem nicht immer zu brauchbaren 
Näherungsaussagen führenden Pohlhausen -Verfahren durchgeführt worden war. Auch wenn 
die Grenzschicht unter den von einer schwachen Wandwelligkeit bewirkten kleinen Druckschwan- 
kungen in einer dünnen Untergrenzschicht in nächster Wandnähe Rückströmung aufweist, kann 
das Pohlhausen -Verfahren nur die aus der einparametrigen Grenzschichtprofilklasse zu ent- 

. nehmenden voll ausgebildeten normalen Ablösungsprofile zur Darstellung dieser Strömung liefern, = 
was selbstverständlich ein im Stromlinienverlauf und im ‚Hinblick auf Stabilität falsches Bild er-. = 
geben würde. 

Neuerdings haben A. W. Quick und K. Schröder [9] für eine Reihe von Beispielen 
den Verlauf der Grenzschichtströmung an welligen Wänden numerisch ermittelt. Sie benutzten 
dazu ein zuvor vonK. Schröder [6]angegebenes numerisches Berechnungsverfahren auf Grund. 
eines Differenzenansatzes für die strengen Differentialgleichungen. Die Zahlenbeispiele von Quick 
und Schröder bestätigen die starke Empfindlichkeit der laminaren Grenzschichten gegen- 
über geringen periodischen Druckschwankungen. Sie zeigen ferner, daß die Wandwelligkeit eine 

; nach außen rasch abklingende Welligkeit der Stromlinien zur Folge hat, d.h. die für die äußere 
Potentialströmung wirksame Stromlinienwelligkeit am äußeren Rande der Grenzschicht tritt gegen- 
über der Wandwelligkeit stark abgeschwächt in Erscheinung. Dieses Abklingen erfolgt im wesent- 
lichen bereits im wandnächsten Teil der Grenzschicht. Die Störungen der Grenzschicht durch die 
Druckschwankungen sind entsprechend im wandnahen Teil der Grenzschicht besonders stark. 
Quick und Schröder folgern aus ihren Beispielen, daß jede noch so schwache periodische 
Druckschwankung schließlich, wenn man die Strömung rechnerisch nur genügend weit stromab- 
wärts verfolgt, Ablösung der Grenzschicht liefert. 2 
Quickund Schröder haben ihre Beispiele numerisch über die Ablösungsstelle hinaus 
verfolgt. Es muß als besonderer Vorzug des Sehröderschen Fortsetzungsverfahrens angesehen 
werden, daß die Rechnung ohne Genauigkeitsverlust über diese Stelle hinwegkommt. Es ist jedoch 
fragwürdig, ob es überhaupt einen Sinn hat, die Rechnungen nach einem Fortsetzungsverfahren 
in Hauptströmungsrichtung bis ins Rückströmungsgebiet hinein fortzuführen. Selbst wenn die 
Grenzschichtvernachlässigungen hier noch zu Recht bestehen, ist, wie schon L. Prandt] gele- 
gentlich zum Ausdruck gebracht hat, zu bemerken, daß die hier angetroffenen Flüssigkeitsteile aus 
stromabwärts gelegenen Gebieten kommen, in denen der Strömungsverlauf noch nicht erfaßt 
“ worden ist, jedoch bestimmend den Strömungsverlauf hinter der Ablösungsstelle beeinflussen kann. 
Es ist also zu erwarten, daß ein genähertes numerisches Fortsetzungsverfahren im Rückströmungs- 
gebiet nicht zu einer eindeutigen Lösung führen kann, sondern bestenfalls eine der noch möglichen 
Lösungen liefert. Aus den von Quick und Schröder berechneten, nach längerem nahezu 
periodischem Verlauf der Grenzschicht plötzlich einsetzenden, zum Teil recht heftigen Rückströ- 


a A 1 | A a ae Li an ai 
\ . s x 


en & 


mungserscheinungen dürften dahe, N ne zuverlässigen Schlüs 
mungsverlauf im Ablösungsgebiet zu ziehen sein IE LES 


kung darin, daß man numerisch nur immer eine Schar von lehrreichen Beispielen gewinnen kan 
Der rechnerische Aufwand bleibt dabei auch heute [7], [8], [9] noch angesichts der großen Zahl le 
 „ erforderlichen Rechenschritte in Fortsetzungsrichtung erheblich. Um allgemeine Gesetzmäßi 
2 keiten und einen tieferen Einblick in die analytischen Zusammenhänge zu erhalten, muß man das 


Die sehr verdienstvollen Rechnungen von Quickund Schrö re B 


> Problem für variable Störungsamplituden, Wellenlängen und Re y noldssche Zahlen zu Be 5 
‘suchen. Einen ersten Versuch in dieser Richtung unternimmt die vorliegende Arbeit. Einsolcker —_ 
Versuch gelingt angesichts der mathematischen Schwierigkeiten des Problems natürlich ‚nur auf Be; 
Kosten gewisser Vernachlässigungen. Diese beruhen im folgenden einmal in einer vorläufigen Be- >> Kae 
schränkung auf so schwache Störungsamplituden, daß eine linearisierte Störungstheorie durch- 


geführt werden kann. (Ein Theorie der zweiten Ordnung ist inzwischen in Angriff genommen 


Er worden.) Ferner werden Schwankungen kleiner Wellenlänge vorausgesetzt, so daß — ähnlich wie 
in der Tollmienschen Stabilitätstheorie der Grenzschichten — das ungestörte Geschwindig- 
_ keitsprofil als gegenüber den Störungen mit x langsam veränderlich angesehen werden kann. Be- 
sonders einfach werden die Verhältnisse, wenn bei sehr schwacher Wandwelligkeit die überströmte 


Wandbogenlänge sehr groß ist gegenüber der Störungswellenlänge. Dann erweist sich die Dicke 
der „„Störungsgrenzschicht‘““ als klein gegenüber der Dicke der gesamten Grenzschicht, so 
daß sich die Störungsvorgänge im wesentlichen im wandnächsten Teil der Grenzschicht, wo das 


‚ungestörte Geschwindigkeitsprofil praktisch noch linear verläuft, abspielen. Dieser theoretische 


Grenzfall wurde zur Gewinnung allgemeiner Gesetze besonders untersucht. Die extrapolatorische 


‘ Anwendung dieser Gesetze auf Gebiete weiter stromaufwärts dürfte ein größenordnungsmäßig 
. zutreffendes Bild der Verhältnisse liefern. Jedenfalls ist insbesondere die Nachprüfung des 


gewonnenen allgemeinen Ablösungsgesetzes an den Quick-Sc hröder schen Beispielen 
qualitativ und auch größenordnungsmäßig quantitativ befriedigend. Die oben geschilderten Eigen- 
arten der Strömung, wie sie durch die Quick-Schröderschen Beispiele aufgewiesen wurden, 


finden durch die vorliegende Theorie ihre Bestätigung, wenn man von den von Quick und 4 


Schröder berechneten Erscheinungen im Rückströmungsgebiet absieht. 


$ 2. Voraussetzungen der Theorie. 


Es seien x und % die üblichen, in der Nachbarschaft der umströmten Wandkontur definierten 
Grenzschichtkoordinaten: Ein Rechtssystem, in welchem x die Wandbogenlänge und % den senk- 
rechten Wandabstand bedeuten. Am ‚äußeren Rande‘‘ der Grenzschicht — dort also, wo die Strö- 
mung praktisch reibungsfrei verläuft — besitze die als stationär angenommene Strömung in z- 
Richtung die Geschwindigkeitskomponente 7, (x), somit in der Grenzschicht das Druckgefälle 
dp,/d& = —0Ü,Ü';g. Die danach auf Grund einer Grenzschichtrechnung sich ergebende, im fol- 
genden als ‚‚ungestört‘‘ bezeichnete Grenzschichtströmung habe in x- bzw. y-Richtung die Kompo- 
ienten U, (2,9), Vo(x, y). 

Wir denken uns nunmehr diesen Strömungsverlauf gestört. Die Geschwindigkeits- 
komponenten der gestörten Grenzschichtströmung seien U (x, y) und V (x, y). Wir schrei- 
ben als neue gestörte äußere Geschwindigkeit am Rande der Grenzschicht vor: 


U(2)=.U,4 8605 Ba 2 (2.1) 


und setzen dabei voraus: a) die Störungsamplitude e und die resultierende Störung der Strömungs- 
geschwindigkeit innerhalb der Grenzschicht sollen vom Betrage so klein sein, daß die nachfolgende 
Linearisierung der Theorie bezüglich dieser Größen gerechtfertigt erscheint; b) über mehrere 
Wellenlängen A — 2/ß der Störung möge sich die ungestörte Grenzschicht nicht oder jedenfalls 


' so wenig ändern, daß man dort überall in der nachfolgenden Theorie T, durch einen konstanten 


Mittelwert D,„ in guter Näherung ersetzen, A D’,/U, als vom Betrage klein gegen 1 ansehen und 
überhaupt in einem solchen x-Intervall einiger Wellenlängen A das Grenzschichtprofil U,= U, (Y) 
setzen kann. Wir beschränken uns, kurz gesagt, auf hinreichend kleine Wellenlängen!). 


Auf Grund der gemachten Voraussetzungen kann man das Druckgefälle längs der Wand- 
kontur durch 


d Ep ID E gi 
= OU = —0{0, 0, —U,ßesnßel .....2..009 
approximieren. Setzt man 
U, y)=Ud,+Eeun&Y, Va y)=Vo ten (ay, . nu Tr (2.3) 


!) An Stelle der durch diese Voraussetzung gesicherten Vernachlässigung der Veränderlichkeit der ungestörten 
Grenzschicht mit & könnte man nach einer Bemerkung von L.Prandtlein geeignetes zusätzliches Kraftfeld annehmen, 
welches gerade einen von x unabhängigen Verlauf von U, bewirkt. Man enthebt sich bei Zugrundelegung eines solchen 
Gedankenmodells der Notwendigkeit, Vernachlässigungen durchzuführen. 


” = Aussagen allgemeiner Art. EX | 

00... Die mit den obigen Voraussetzungen erzielbaren Vereinfachungen erlauben es, A: Störu 
„mit einem angemessenen Aufwand analytisch zu erfassen. Durch einen in x periodischen Ansatz . 

für diese Störung in Form einer Fourierreihe, von der für kleine & je nach Grad der Näherung £ 


Differentialgleichungen geführt. Wir werden uns in der vorliegenden Untersuchung auf eine :2 3 h 


 Druckgefälle, nämlich 


A'zu berücksichtigen. (Wenn A || <e angenommen wird, statt wie bisher nur A DAR << 
und |e/7,|<&1 anzunehmen, kann der 3. Term für eine Theorie der 2. Ordnung wieder weggelassen 


'($5) annehmen, daß die Störüng durch eine schwache Aufwellung der bisher betrachteten Wand- 


‚gewählten Grenzschichtkoordinaten x und % der ungestörten (,‚glatten‘‘) Wand beibehalten werden. 


Bar c“ \ 


1 oraussetzung b) einen in x T eriödischen ‚Ver ieser St 
wur als Approximation des Verlaufs für das betrachtete klei vall \ 
at der Rechnung erlaubt jedoch durch Variation der auftretenden. Parameter en ül 
ein beliebiges solches Intervall der Grenzschichtströmung und Auhre gamıl zu pn 


bezüglich e nur wenige Beiwerte ermittelt zu werden brauchen, wird man auf lineare gewöhnliche “ 


Theorie der 1. Ordnung in e beschränken. Auf diese ist ja auch die Approximation (2.2) bereits 
zugeschnitten. Für‘eine Störungstheorie der 2. Ordnung wäre der vollständige Ausdruck für GE Be, 


- 


‚dp 


ee U) — Ouße sin Bc + Diecos ßx  Spsin2 pe] BE. a 


werden. In den Störungsgeschwindigkeiten bleiben dann die aus der linearisierten Theorie berech- 
neten Fourier-Terme unverändert, es kommen nur weitere bisher un Terme additiv 
hinzu.) $ 

Es ist noch nichts che gesagt worden, wie nun eigentlich die angenommene Störung 
verursacht zu denken ist. Diese Frage soll auch noch zurückgestellt werden. Wir werden später 


kontur bedingt wird. Unbeschadet dieser nachträglich zu vollziehenden Deutung sollen aber die 


Der Leser wolle sich vorläufig vorstellen, es handele sich in (2.1) um einen von außen her vorge- 
schriebenen Verlauf der Geschwindigkeit am äußeren Rande der Grenzschicht an der glatten Wand 
(etwa durch entsprechende Lenkung des Strömungsverlaufs außerhalb der Grenzschicht). Das 
Ergebnis der Rechnung wird uns zwingen, diese Vorstellung aufzugeben. Eine bekannte Eigenschaft 
der Lösungen der Grenzschichtgleichungen gibt uns aber die Möglichkeit, leicht den Übergang zum 
eigentlichen Problem der Grenzschicht längs einer schwachwelligen Wand zu vollziehen und die 
Ergebnisse der Rechnung zutreffend auf diese anzuwenden. 


$ 3. Einführung von Dimensionslosen. 


Es bezeichnet im folgenden ö, = const ein geeignetes Maß für die Dicke der ungestörten B 
Grenzschicht im betrachteten x-Intervall. Wir führen folgende Dimensionslosen ein: S 
Er 07 = RE), y* — yJöo = 

UN UIO,; W= RE Pa Dr TR 3.1 E 

U*= U/D, + = V/ßöUm  &*=e/Un $: 

Stellt y(x, y) die vermöge U = da y/ay, V—.— 9y/dx definierte Stromfunktion der Grenzschicht- . 


strömung dar (entsprechend y, (x, y) diejenige der ungestörten Grenzschichtströmung), so defi- 
nieren wir ferner die dimensionslose Stromfunktion 


vr = Yyld, Um ee ee (3.2). 


Schließlieh führen wir noch I Zahl o vermöge 
DO, 
une - LE N ee En (3.3) 
[ 


ein. 

Neben den dimensionsbehafteten Konstanten ß, A, Öy, DT, o und ».und den dimensionslosen 
Störungsgrößen «, und v, werden im folgenden nur die in (3.1) bis (3.3) definierten dimensions- 
losen Größen Verwendung finden. Wir lassen daher weiterhin die Sterne 
wieder weg. 


$4. Herleitung der Grundgleiehungen. Allgemeine Lösungsansätze. 


Indem die Stromfunktion den Rechnungen zugrunde gelegt wird, ist die Kontinuitätsglei- 
chung von vornherein durch den Ansatz befriedigt, Die allein noch zu berücksichtigende Grenz- 
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schichtdifferentialgleichung lautet in den in $ 3 eingeführten Dimensionslosen unter Berücksichti- 
gung von (2.2) 


YıYan — Va Yyy Yu = — D, Der Bee BER: (4.1). 
Die Stromfunktion y, (z, y) erfüllt (4.1) für e= 0. Macht man den Ansatz Ei 
y=yvlsy)telig)esctn)sine, ...... er 


so erhält man daher unter Weglassen aller Glieder von höherer als erster Ordnung in & 
U,{—f/sinz + gicos a} + U,{fsin x —g, cos x} —o {fi cos «+ gi’ sin 2} = —sinw.. | 


- Dabei ist nachträglich U, = U,(y) gesetzt worden und Striche bedeuten Ableitungen nach y. 
Es ergibt sich somit zur Bestimmung von f(y) und g, (y) das System von Grundgleichungen: 


UWi-Uy)hteg —1 RE 

UV) U, y)n—eh”—=0 ’ ; 

Die Randbedingungen lauten: 

AN)=FLM)=0, Ft) hr 
1 )=g,0)=g(w)=0. | 
Ohne die Gestalt des ungestörten Geschwindigkeitsprofils U, (y) in spezieller Weise festzu- 
legen, sollen zunächst Lösungsansätze allgemeiner Art mitgeteilt werden, welche die numerische 
Berechnung von f, und g, für die interessierenden Parameterwerte in jedem Falle ermöglichen. 
Der Allgemeinheit dieser Aufgabenstellung entsprechend ist man auf die bei Grenzschichtrech- 
nungen und auch sonst gebräuchliche Methode angewiesen, Reihenentwicklungen von der Wand 
= O0 her nach Potenzen von y mit asymptotischen Lösungen für große Wandabstände numerisch 
abzugleichen. 
Die Reihenansätze 


_ yhy ou W 
E Ve FEN 91 MT RE ya cn (4.4), . 
welche bereits alle Randbedingungen in y= 0 erfüllen, geben zusammen mit der vorgegebenen 


mindestens für y < |y,| konvergenten Entwicklung 


bei Einführung in (4.3) und nachfolgendem Koeffizientenvergleich 
AY)= RAY) + Y2As(y) + As(y) 


EN ER 4.6 
1(y) = —RaA;(Y) + Yadı(y) + B;(y) J = 
‚mit den für o +0 und |y| < |ya| konvergenten Entwicklungen 
| _Y. Ay ma, y 
Als ee : 
AED WI 20 Ye DIR IE 
; 0.512... 0-71 v0.8108020} ’ 
4,(y) = REEL SE. a 
N a Aa rel “ 
1:y° 104? y? 
B,(y) = re aa ern 


Die Parameter Ra und Y. bleiben frei. Es sind ja erst vier Randbedingungen erfüllt worden. Wir. 
beschränken uns hier auf diese Angaben. Für den numerischen Anschluß an asymptotische Lö- 


Be wird man im allgemeinen weitere Glieder berechnen müssen, was keine Schwierigkeiten 
ereitet. 


Für hinreichend große Wandabstände y verläuft die Strömung praktisch reibungslos. Hier 
draußen ist also in erster. Näherung 


uf, — U, = 1, 
U, g, Zr: U, =0 
mit den den Randbedingungen in y= © genügenden „reibungslosen Lösungen‘ 


d 
Awm=v, | +Av, ee 
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(A, B verfügbare Konstanten). Iterative Berücksichtigung der Reibungsglieder führt auf Potenz- 
entwicklungen nach der mit der Zähigkeit' proportionalen Zahl o. Unter Einführung von 


Y 
5 d 
RER (4.7) 
Y os 
und des Operators : 
B Y 
F LE 
an! = Us v3 I Re E22 (4.8) 
Ya 


(Ya > 0) erhält man als erste Näherung für den Reibungseinfluß, d. h. bis zu den linearen Potenzen 
in o: 


A) = fly) + AU,(y) —c BJ (U,) | a 
9(y) = BU,(y) +cAJ(U,)+0J(f) I (J (g))) 


In Klammern sind die Fehlerglieder angegeben, mit denen die Iteration leicht fortgesetzt werden 
kann?). Wie es sein muß, ist- wegen Au U,(y)=1:lim’f, (y) = Zahl y, 9 (y) = Zahl. 
Statt von den „reibungslosen ‘‘ Dolungen aussugehen: hätte man die Iteration® auch auf Grund 
dieser asymptotischen Beziehungen, die ja nur die asymptotisch zu verstehenden Randbedingungen 
in y= © darstellen, in Gang setzen können. Man hätte daraus dann umgekehrt die Belang- 
losigkeit der Reibungsglieder, für eine erste Näherung weitab von der Wand folgern können. 

Die eingeführte feste Integrationsgrenze y= y, > 0 bezeichne nunmehr die Anschlußstelle 
der asymptotischen Näherungen an die von y = 0 aus entwickelten und bei einer genügend hohen 
Potenz abgebrochenen Reihen. Aus (4.9) für y= y, folgt: 


A = A (y)/Uo( (Ya) Begt YO al) ha (4.10). 


Bildet man fi (y,) und gi (y,) nach der ersten asymptotischen Näherung (4.9) und jetzt. man die 
eben gewonnenen Ausdrücke für A und B ein, so wird 


r ; U } 
Gh uh-irer =0 
0 


in ya Ya: 

u Ir: 

Un —- on -o0-- Ho 
os1 1 U, Jı U3% 


Setzt man hier die aus den abgebrochenen und für y= y, berechneten Reihen (4.6a) folgenden 
Werte für f,, fi, 9, und gi ein, so erhält man zwei lineare inhomogene Gleichungen für die noch 
zu ermittelnden Koeffizienten 8, und y;. Sie haben die Gestalt 


aß rt ap=dı) 
N TER 4.11 
— 0, ßa + ya dh] 


und liefern (e2 +.c3 +0) die gesuchten Zahlen ß, und y,, womit die Reihenentwicklungen und 
damit inabesöndere die daraus folgenden Näherungswerte für f, (y,) und g, (y.) festgelegt sind. 
Damit sind nach (4.10) auch A und B vollkommen bestimmt. 

In den a enngen müssen soviele Glieder berücksichtigt werden, daß man den 
Anschlußpunkt y, in einem endlichen Intervall, wo die abgebrochenen Reihen noch eine gute 
Näherung darstellen, variieren kann, ohne die Zahlenergebnisse für die Integrationskonstanten 
ß, und yzim Rahmen der gestellten Genauigkeitsansprüche zu ändern. Dies muß ja der Fall sein, 
wenn in diesem Interyall die gewählte asymptotische Näherung bereits ausreicht. 


—U 


$5. Grenzschiehtverlauf an einer schwach gewellten Platte (numerisches Beispiel). 

Die bisher entwickelte Theorie soll hier zunächst an einem konkreten Beispiel bis zum nume- 
rischen Endergebnis durchgerechnet werden. Als ungestörte Strömung wird die von H. Bla- 
sius [1] berechnete Grenzschicht an einer längsangeströmten Platte gewählt. Hier ist U,— 
const. (= 587 

Es sei L die überströmte Länge von der Plattenvorderkante bis zu einer mittleren Stelle 
des im folgenden zu betrachtenden x-Intervalls der periodischen Störung. Mit 

a (5.1) 
als für die Rechnung zweckmäßig ee Maß für die mittlere Dicke der ungestörten Grenz- 


») Man erkennt, daß man so für f, (und analog für g,) Glied für Glied die aufsteigenden Potenzen des auf die aus- 
geschriebene Näherung 2, + 4AU,— 0 BJ(U,) angewandten Operators — 0? J? erhält, 
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schicht in der betrachteten Umgebung der Stelle Z, erhält man wegen o!= ßög' U, 6o/? 
EBEN 2m. LjA oe ee (3.2). 

Bis auf einen konstanten Zahlenfaktor ergibt sich somit als charakteristischer Parameter o der 
Störung das Verhältnis von Störungswellenlänge zur überströmten Länge, also eine rein geome- 
metrische Kennzahl. Die folgenden Ergebnisse sind daher in weitem Umfang von der Reyr 
noldsschen Zahl unabhängig (im Rahmen natürlich der für diese Theorie allein zulässigen hohen 
Rey.noldsschen Zahlen). ; EEE 

Bezeichnet ö, die tatsächliche Dicke der ungestörten Grenzschicht nach Blasius an der 
Stelle L und ö, diejenige in L +4, so ist wegen 


{ Ös ur Fu HE = yı + D.TE 02% . je Pe N ET ( ) 
vorauszusetzen, daß i 
a WE 3 ae RELEEE u 
27 


(d.h. A<L) ist. Dann ist nämlich der relative Dickenzuwachs der ungestörten Grenzschicht pro 
Wellenlänge den Voraussetzungen unserer Theorie entsprechend klein und man kann über einige 
Wellenlängen U, = U,(y) setzen. Für das im folgenden zu untersuchende Zahlenbeispiel wählen 


wır 
a er (5.5). 


Wegen Ö,/ö, = 1 +ro ist der Zuwachs pro Wellenlänge hier etwa 3%, also erträglich klein. Die 
Wahl o = 10? bedeutet: Unsere Untersuchung bezieht sich auf eine Wandstelle, die in L= 100 4/2 z, 
d.h. ca. 16 Wellenlängen stromabwärts von der Plattenvorderkante liegt. Die Wellenlänge A 
bleibt frei. Te 

Amplitude e und Reynoldssche Zahl brauchen für die numerische Durchrechnung der 
Störungstheorie nicht zahlenmäßig festgelegt zu werden. Erst für das Bildmaterial werden be- 
stimmte Zahlenwerte angenommen, insbesondere e = 10°? und.-5 - 103, 


Nach Blasius a.a.0. (siehe auch etwa L. Prandtl1]2)) ist 


ER HN Mur 
U — = N Re ET 
Il. 6.) 
mit 
= 0,33206 2,0 ne Re a (5.6. a) 
(dies nach L. Howarth [4]) und 
Gl, Gl; G=1l, Or 3% 0 = A189 22 ro 
Im Ausdruck (4.5) der allgemeinen Entwicklung für U, (y) ist also nun speziell 
1 11 375 I 
= &, yu=—zh, %=.Ta, Au=— 3 Re Ds (9.7), 
während alle anderen a,(v+3r+1; n=0,1,2...) verschwinden. Das vereinfacht die 


FR weitere Rechnung erheblich, soweit sie sich auf die Ermittlung der 
Reihen (4.4) bezieht. Insbesondere verschwinden alle ß3, + ı und 
Y3y+1.für v= 1: 2 

Die numerische Rechnung im einzelnen, die zur Ermittlung 
der Parameter ß,,y,, Aund B führt, mag hier nur kurz gestreift 
werden. Sie ist trotz aller Vereinfachungen des vorliegenden 
Spezialfalls noch langwierig. Zwar ist die rekursive Berechnung 
der Koeffizienten in den Reihen (4.6a) hier besonders einfach, 
Er man muß aber für eine genügende Approximation bis zuhohen 
Plattengrenzechicht nach H. Blasins Fotenzen gehen. Wir haben bis zu den Gliedern mit y!® ge- 
und Ableitung. Verdrängungsdicke rechnet und konnten dann den Anschluß an 'die asympto- 


1 2 3 4 y 


6*= 1,72 6, tischen Lösungen in y,=1, d.h. im Wandabstand ö, Vvoll- 
TR ziehen (siehe hierzu Bild1). Der für die Verwendung der 


asymptotischen Lösungen (4.9) benötigte Verlauf von er 
’ d e ee | 
J(U,), J(U,) = ir J w). J (jı) und J’ (j,) mit y, = 1 ist in Bild 2a bis 2e dargestellt. Die 
Rechnung lieferte für unser o = 10°? bei einem zugelassenen Rechenfehler von 1% 


Aa= 21,08 »=—37,3 und A-2283, B--1915,...08), 


ee an [Dez. 1947 Gö rtler, Einfluß PLıor schwachen u Wandwelligkeit ) 


Eine Kontrollrechnung wurde in folgender Weise durchgeführt. Einmalige Integration der 
Differentialgleichungen (4.3) liefert 


22 1 Y ’ 
Hr =, Un —2] Duady), 


ER 


also für y— oo wegen U,(%) = 


fie er. (5.9) 
& Are Us ayı 


9) +2 Uıg ayı, 


. 
= [tim ( DU au). 


Bild2a. Verlauf von f, und df,/dy Bild2b. Verlauf von J(f,) und Bild 2c. Verlauf von J(U,) und 


mit Y, =1 als Funktion von y. aJ(fı)/day mit Ya =1 als Funk- d4J(U,)/dy mit v„ =1 als Funk- 
tion von Y. tion von y. 


Aus den asymptotischen Lösungen (4.9) ergibt sich wegen 
lim (fi — 4) = 3.825, lim J (U,) = —0,2093, lim J(7)=—1,662 . . . (5.11) 
y>» y>» Yy>& 
nun 
lim (A —y)= 443,825 + 0 B:-0,2093 
Ve SE (5.12), 
9 (o)= B—o : 1,662 — 0A : 0,2093 

also nach (3.10) und mit o= 10? 


—0,002093 A + B= 0,01662 — 0,01 BE 27 U! g,dy a a 


Nimmt man die Forderungen des ste- 
tigen Anschlusses von f, und g, in 
Y= Y„ = 1 hinzu, die A und B linear 
in ß, und y, ausdrücken, so bekommt 
man bei Einsetzen derselben in (5.13) 
zwei Gleichungen für ß, und y,. Setzt 
man nun in den Koeffizienten, insbe- 
sondere in den Integralen rechts in 
(5.13) den bei der bisherigen Rechnung 
numerisch ermittelten Verlauf von /f, 
und 9,in 0 <y < wein, und berechnet 
erneut ß, und y,, so erhält man wieder 
Ba = 21,08 "und, = — 37,33. Diese 
Zahlen dürften also mit der angege- 
benen Stellenzahl richtig sein, was 
mehr ist, als oben beansprucht wurde. 
Angesichts des Umfangs der früheren 
numerischen Rechnung war diese Probe 
jedenfalls eine erwünschte Bestäti- 
gung. In Bild3 ist der füro—= 10° Bila3. Verlauf von f, //, 9, und g/ in Abhängigkeit von y für 
berechnete Gesamtverlauf von f,, !/,9ı «= 10. 
und g; dargestellt. 

Allgemein sehen wir, daß £,; und y, nur von o, nicht aber von Re abhängen. Da in unseren 
Dimensionslosen 


U (x, Y) = un (Y) = € U, (%, Y), 
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also 


a Rt RN ee 


5% (x, 0) = 0,33206+ e{ ß, cos «+ Y,sin « } 


ist, hängt der Ort x der Ablösung, wenn eine solche eintritt, nur von dem Verhältnis e der Stö- 
rungsamplitude der äußeren Geschwindigkeit zur ungestörten äußeren Geschwindigkeit U, und 
von dem geometrischen Verhältnis o!= 2r.L]A ab. Dieses Ergebnis trügt jedoch. Bei der Deu- 
tung unserer Ergebnisse als Grenzschichtströmung an einer geeigneten welligen Wand =: Wir 
haben bisher nur eine ebene Platte bei periodischem Druck betrachtet und werden diese Hilfs- 
vorstellung bald aufgeben müssen — gewinnt auch die Reynoldssche Zahl Bedeutung, ins- 
besondere hängt & bei vorgegebener Wandwelle vonder Reynoldsschen Zahlab. ar, 

Für o=10?insbesonderetritt Ablösungnach den obigen Zahlenergebnissenerstmalig beiwach- 


du 
sendem e dann ein, wenn e den Wert 0;00775 überschreitet. Für e= 0;00775 Nissen (x, 0) gerade 


Null in 180°&/2r=119,5°. Für e=0,01 beispielsweise tritt Ablösung bei 180° x/2r 80° 
ein und Wiederanlegen erfolgt bei = 159°. Für beliebiges e und o = 10”? liest man entsprechende 


5 0 
Aussagen leicht nach (5.14) unmittelbar aus dem in Bild 4 dargestellten Verlauf von er (2,0) ab. 


Um nunmehr ein Bild von derjenigen 
Strömung an einer welligen Wand zu gewinnen, 
welche durch unsere Ergebnisse dargestellt wird, 


() -330-u)|,,- , 


Bild 4. Verlauf von 
öu,/@y längs der 
Wand -(y = 0) für 
0, — 100°. Zaun 
Phasenvergleich, 
darunter die 
Druckschwankung. 


machen wir uns die folgende bekannte Eigen- 
schaft der Lösungen der Grenzschichtdifferen- 
tialgleichung zunutze (siehe hierzu etwa [3]): 
Stellen U (x, y), V (x, y) die Geschwindigkeits- 
komponenten einer Grenzschichtströmung dar, 
so sind auch 


U*(2,y)=U(®,y+f(z)) . . (5.15) 


und der Kontinuität gemäß 


— (a) U y+f@)f 


Geschwindigkeitskomponenten einer Lösung. 
Dabei ist vorauszusetzen, daß die an sich will- 
kürliche Funktion f(x) derart verläuft, daß die . 
grenzschichttheoretischen Vernachlässigungen 
bezüglich der x-Richtung als ausgezeichneter 
Hauptströmungsrichtung auch für die neue 
Lösung zulässig sind und ferner, daß die an 
diese Grenzschichtströmung _anschließbare 
äußere Potentialströmung den vorgegebenen, 
der Grenzschicht eingeprägten Druckverlauf 


Bild 5. Verlaufdes .. 
liefert. 


Störungsanteils 
ger re Auf Grund dieses Sachverhalts ist der Über- 
eo 1ürk. gang von der bisherigen vorläufigen Vorstellung 


einer geradlinigen Wandkontur (z-Achse) zu der 
wirklichen schwach welligen Kontur vorgezeich- 
net: Dieser Übergang hatso zu erfolgen, daß die 
sich ergebende Stromlinie am äußeren Rande der 
berechneten Grenzschicht mit jener überein- 
| stimmt, welche dort eine Potentialströmung mit 
dem von uns angenommenen Druckverlauf besitzt. Nennen wir diese letztere Stromliniey = 49 (8). 
Schwankt die Verdrängungsdicke ö* der Grenzschicht vernachlässigbar wenig gegenüber der 
Schwankung von %, (x), so wird man als Wand die um den Mittelwert von ö* gegen y= y,(%) 
nach rückwärts versetzte Parallelkurve zu wählen haben. (Bei den meisten Problemen, mit:denen 
sich die Grenzschichttheorie üblicherweise zu befassen hat, bleibt 6* überhaupt klein gegenüber 
den maßgebenden Abmessungen des Körpers, so daß man auch von diesem Übergang zu einer 
Parallelkurve absehen und y= %,(x) als Wandkurve annehmen kann.) Ist jedoch die Schwankung 
von ö*(x) mit jener von y,(x) vergleichbar oder gar groß im Vergleich zu dieser, so muß man von 
y=y,(x) als Stromlinie am äußeren Rande der Grenzschicht um ö*(x) nach rückwärts gehen 


Ce ES 


„Für de Stromlinie 0 erhält man ‚daher in erster Näherung für Kae & 
ar .y= a ne ER IE 
3 Wegen. Pi= Up Reit d): nit Re(L)=U, oLiv. U ; 
Es  yo(®) =c0 Rel2(L)cos ®. 
3 = ‚Hieran möge Be eine Zwischenbetrachtung angeschlossen werden. +2, 
Br Das Geschwindigkeitsfeld der Potentialströmung ist in y-Richtung auf den ten Teil 


 geklungen, wenn man in Richtung wachsender y den Weg 1/ß = A/2n AEBeRBe hat. Nur 
beträgt die Verdrängnuigsdicke ö, der ungestörten ae) u 


= 5 | | 8, = 1,7208 d, z= re ER ol 


_r .der Grenzschicht in der ne ee noch über, einer , Breite. in : 
Be: y-Richtung praktisch annähernd volle äußere. Potentialströmungsgeschwindigkeit herrscht 3). 
Be; "Anders gesagt: Man muß fordern, daß die Dicke der Zone, in welcher die äußeren Geschwindig- 
ee - en abklingen, groß ist gegenüber der Dicke der Grenzschicht selbst. Nun ist 


= S 2 BO la! ReslD) ann 2, ee . 5.19). 


Pe Für a= 107 rd Re (L) = = 10° wird ßö,= 0,172, also schon etwas groß, aber immerhin noch 
‚größenordnungsmäßig im Rahmen unserer auf eine erste Näherung zielenden Fiechnung einigenz 
 _maßen befriedigend. 

e Man kann nun (siehe oben) in (5. 15) nur dann f(x) = y, (x) + const. eine wenn gegen- 
u ‘ über der Schwankung der für die Potentialströmung angenommenen Wand (y, (x) nach (5. 17)) 
E35: die Schwankung der Verdrängungsdicke ö* der gestörten Grenzschicht vernachlässigbar klein ist. 


2 — Die Verdrängungsdicke der gestörten Grenzschicht berechnet sich zu: 
+ — % (0-0. | 
f U Y 
) = 
Franz u 7208 2 nn (y—Fı(y)) cos —g, (@) sin a]) 


also, wegen (5.12) und unter Berücksichtigung nur des linearen Ausdruckes in ein der Entwicklung . 


von 1/(1-+ ecos «) 
6*—= d,+ 8 6,{ — (5,546 + A + 0,2093 0 B) cos x 


—(B—1,662 0 —0,2093 0 A) sin «}) ' ee 
Insbesondere ergibt sich für unser 'numerisches Beispiel o = 10 (hierzu siehe Bild 5) 
& ee F- 
| ee ed Ött+ ef 28,33 cos ©+19,21sin x} | . (5.208). 


— 6, {1,7208 —34,23 & cos (e+34,7°)Y J 


(Die Extremwerte von ö*(x) werden in 180° x/2 = 145° bzw. 325° angenommen.) Die Gesamt- 
schwankung der Verdrängungsdicke ö* (x) beträgt demnach 


a ad aan a N NN ale 13,3 357 2 1," Alan ar nl + ln nel Da a a NR 
un x \ ER, u Si 
’ f 
A 


| ee en, 

Nun beträgt aber nach (5.17) die Schwankung der Linie %, (x) (‚‚potentialtheoretische Wandwelle‘“) 

I a en mn......2...0, Se 
somit für o= 10? 5 

er lan Rer dm ....2...2.:062) 


») Hierauf machte mich Herr Prandt]l aufmerksam. 
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Berae N a ee. in Furbydeite ee Re on We erteilen würde 
über den kritischen Wert. hinausgehen. Selbst noch für Re= = 10° betrüge das a 
i meta etwa ein Drittel. 


Die wirkliche Ba weicht somit gewiß erheblich v von y= Yo (@) = ‚7208 ab. 


ur online y=yYy ( &). Diese muß als Stromlinie am äußeten Ra der ren iehleh ‚gewi ıl 
werden, damit die von außen der Grenzschicht aufgeprägte Druckvertei ung der Potentialströmung 
- mit der von uns angesetzten übereinstimmt. Von diesem Rande nach innen hat man also di 
- berechnete Grenzschicht anzuschließen. Als Wandstromlinie resultiert dann die gegen die Strom- 
linie I =% em um die Verdrängungsdicke ö* (x) nach rückwärts RER Kurve. Somit wird R 


Be: | Ce or I0EE es 62 


Die Wandwelle wird insbesondere Mo 1072 nach (5. 20a) = re 
_Yw(x) = 0,01 Reu2; cos &— 1,7208 + e (28, > cos c—19,21 sin.«) 


—1,7208 + 55” 00 (+9) a TS 


‚mit Ay=e 3 [(0, 01 Reil2 1 28,33)2 + 19 ‚2122 | 
"und tg o—= 19,21/(0,01 Re’2- 28,33) ee ar Kö 


| Beispielsweise wird für Re= 10°: Aw = 42,87 ed, und 9= 26,6° und für Re= 105: Aw=36, 89 € Ö, 
"und 9=31,4°. Da die y-Abmessungen mit ö, dimensionslos gemacht wurden, ist die Länge Pr 
die Amplitude der Wandwelle. Diese ist Sean, ß als vorgegeben anzusehen. Aus ihr resultiert als 


15 


1 re 
Bild 6. Verlauf von u, (x, y) im Falle o = 10-?in Abhängigkeit von Bu 
y für verschiedene im Winkelmaß angegebene z-Werte, 


*) In Übereinstimm R h 
Re ung mit den Ergebnissen der von Quick und Schröder [9] durchgerechneten Zahlen- 
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Amplitude A,=e/ß der ‚„wirksamen‘‘ Welle x) am äußeren Rande der G i 
ee yo) renzschicht wegen 


A,— Ay 0,01 Re\2] [(0,01 ReV24 28,332 + 19212? . . ...... (9.26), 
insbesondere für Re= 10° wird A,= 10 Aw/42,87 z0,2Aw, für Re= 105 ergibt sich A, = 0,1 Ay. 


0% 02 et 10 12 
Bild 7. Geschwindigkeitsprofile U(x, y) für co = 10-2, sr 102% 

Die vorliegende Theorie hat somit schon in Übereinstimmung mit den Quick-Schrö- 
derschen numerischen Beispielen das sehr beruhigende Ergebnis geliefert, daß eine Wandwellig- 
keit durch die Grenzschicht hindurch stark geglättet (und übrigens auch phasenverschoben) als 
„wirksame“ Welligkeit der äußeren Strömung in Erscheinung tritt. 

Zum Verlauf der Strömung im einzelnen geben einige Bilder Aufschluß. Bild 6 gibt den 
Verlauf von u, (2, y) (das Geschwindigkeitsprofil ist, siehe oben, U= U,(y) + eu, (%, y)) für ver- 
schiedene in Graden angegebene x- 8 B ß 
Werte. Eswurden ferner die Beispiele E 
€ —= !/ıo, (mit Ablösung, siehe oben) 
und e=!/,,. zahlenmäßig ausgewer- 7 
tet. Die Bilder 7 und 8 geben die y 
zugehörigen Geschwindigkeitspro- 
file, die Bilder 9 und 10 den Strom- 
linienverlauf. Für den Stromlinien- 
verlauf ist Re (L) = 10% angenom- 5 
men. Man erkennt, daß im Falle 
e= 1/00 Ablösung zwar eintritt, je- 
doch nur in einer sehr schmalen 
Unterschicht der Grenzschicht (vgl. 
y-Maßstab mit Bild 7 bzw. 81); die 
Geschwindigkeit im geschlossenen 
Wirbel des Ablösungsbereichs ist 
praktisch Null. 2 

Die vorstehenden Rechenergeb- 
nisse werden durch die Betrachtun- 
gen des folgenden $ 6in vollem Um- 
fang bestätigt werden. Im ganzen 
kann gesagt werden: Die Ergebnisse 


geben zwar schon manchen Auf- 02 0,4 05 08 10 
schluß, daaber des numerischen Auf- 0 0,2 0,4 06°, 08 1,0 
wandes wegen nur-der Fall o= 10? Bild 8. Geschwindigkeitsprofile U(z,y) für a = 10-:, 8 = 0,5 - 10-2. 


5) Quick und Schröder prägten dieses Bezeichnung [9]. 
16* 
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ausgewertet wurde (eund Re (L) 


—., 


beliebig), bekommt man noch nicht den gewünschten analytisch 


vermittelten Gesamtüberblick. Es muß also eine Vereinfachung der Theorie angestrebt werden, 
die die Bindung an eine Auswertung für einen oder wenige einzelne o-Werte behebt und allgemeinere 


N ee 


nn 


_ 


Aussagen ermöglicht. Ein gangbarer Weg wird im folgenden Paragraphen beschritten. Ferner 
muß gefragt werden, inwieweit die bezüglich e linearisierte Theorie bei Störungen von der Stärke 
der in obigen Bildern für o = 10”? dargestellten noch zuverlässige Werte liefert. Zweifellos werden 


et 


nn 
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Bild 9. Stromlinienverlauf für o = 10-2, e=10-2, Re(L) = 10°. Überhöhung 5 r/3-( = 5,24) fach. 
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Min. 


a 
0° 30° 60° 99° 120° 150° 180° 210° 240° 270° „ 200° 930° 360° 


4,4 
40 
3,6 
32 
2,8 


Bild 10. Stromlinienverlauf für o = 10-?, e = 0,5 -10-?, Re(L) =10°. Überhöhung 5 n/3-(= 5,24) fach. 


die Glieder zweiter Ordnung in e hier schon merklich ins Gewicht fallen. Die Ergebnisse sind jedoch 
sicherlich qualitativ und größenordnungsmäßig quantitativ zutreffend. Mehr wird man von seiten 


der Praxis von diesen Betrachtungen gegenwärtig kaum fordern. 


Eingegangen im April 1947. 


Fortsetzung im nächsten Heft. 
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Bemerkungen zu einer Entwicklungsformel aus der Theorie 
| der Zylinderfunktionen. 
Von Herbert Buchholz in Heidelberg. 


In der vorliegenden Arbeit werden für die Differenz der Produkte zweier linear unabhängiger Zylinder- _ 
Junktionen ungleichen Argumentes, die noch mit den Quotienten zweier Bes sel schen Funktionen behaftet 
sind, und für einige damit verwandte Formen mittels des klassischen Verfahrens der Residuenrechnung die 
ihnen zugehörigen Partialbruchzerlegungen hergeleitet. Dabei wird auf die fehlerhafte Angabe einer solchen 
Entwicklung in dem Buche von G. N. Watson und in der Formelsammlung von Magnus undOber- 
hettinger hingewiesen. 


In the present treatise, partial fraction developments are deduced for the differences of the products of 
two linearly independent Bessel funktions with unequal arguments multiplied. with the quotient of two Bessel 
Junctions and for some other similar forms, by means of the classical method of the residue calculus. Besides, 
an incorrect development of that kind given in the book by G.N. Watson and in the formulary by 
Magnus and Dberhettinger is indicated. 


La these sous les yeux a pour but de deduire pour la difference des produits de deux fonctions de eylindre 
bineavres et absolues de plus encore chargees des quotients de deux fonctions de Bessel et de quelgues especes 
approximatives, au moyen de la methode classique du calcul des residus la decomposition des fractions par- 
tielles en considerations An faisant cela lattentions est fixee sur le rapport faux d’une telle demonstration 
en dans son livre par G.N. Watson et möme dans le recueil de formules d Magnus et d’Ober- 

ettiinger. : 


B noneskamefi yCMOTpeHnm padore, MIA PASHOCTH IPOMsBeNeHHÄ NBYX JIHHeAPHO HesaBHu- 
CHMEIX PYHKIHH UMAIHHAPa PasıImYHBIX APTYMEHTOB, KOTOPEIE KPOME TOTO EINE ABIAWTCH OTBET- 
CTBEHHEIMM 34, KOHPHNHSHT ABYX dyHrumü Beccena, u ANA HeKOTOPHX PONCTBeHHBIX HM 
POPM IIPOBONUTCH KIIACCHYECKHM CHOCOOOM, Ho Merony PesunyeH, ucyncseHune u pas1o:keHue 
HPUHANTESKAIMHX K HHM YACTHYHBIX APodeii. 

IIpu s9Tom cıenyer 0O6paTuTb BHUMAHNe Ha besommmÖbOoUHOe YRABaHNe IIONOÖHOTO PA8BHTUA B 
xHure T.H. Barcon a B cöopnuze dopmyı Markyca u Odeprertriunrep. 


1. Einleitung. 


In dem bekannten Buch über Zylinderfunktionen von G.N.Watson [1]!) wird auf S. 499 
für reelle Werte von x und X mit 0=Zx<X<<1 und beliebige komplexe Werte von 2 und » eine 
Entwicklung für den Ausdruck (1.1a) angegeben, die als eine Partialbruchentwicklung nach den 


Pu) = N (Ku) Al) FL ee SE 
By) = HE (Ku) I, u) N (u) IutK u) en, 


Nullstellen des Nenners aufgefaßt werden kann. Es scheint bisher nicht bemerkt worden zu sein, 
daß diese Entwicklungsformel an der angeführten Stelle nicht bloß falsch wiedergegeben worden 
ist, sondern daß der fehlerhaften Angabe eine irrtümliche Behauptung in der Beweisführung zu- 
grunde liegt. In den beiden an dieser Stelle von W at.s o n zitierten älteren Arbeiten von A.Som- 
merfeld[2jundH.S.Carlslaw [3], in denen die diskriminierte Formel auf anderem Wege 
gewonnen wird, ist sie in ihrer richtigen Form zu finden. Da das Watson sche Buch eine große 
Verbreitung gefunden hat und die besagte Formel in ihrer fehlerhaften Gestalt auch in die erst 
neuerdings herausgekommene Formelsammlung von Magnus und Oberhettinger [4] 
aufgenommen worden ist, so erscheint es dem Verfasser, der diesen Fehler schon vor einigen 
Jahren im Laufe eigener Arbeiten auffand, nunmehr doch an der Zeit, auf ihn hinzuweisen. 


Aber ein weiterer Umstand, der inzwischen hinzugekommen ist, läßt es geraten erscheinen, 
die Frage noch einmal von Grund auf anzupacken. Es läßt sich nämlich ohne viel weitere 
Mühe eine ganz ähnliche Formel für die unter (1.1b) angeführte Funktion herleiten sowie noch 
für zwei andere Varianten dieser beiden Formen, die wir erst später angeben werden, und da 
diese vier Funktionen bei der rechnerischen Behandlung des Ausbreitungsvorganges elektro- 
magnetischer Wellen durch hohle Leiter [5] eine große Rolle spielen und die Kenntnis ihrer Par- 
tialbruchentwicklung hier das beinahe einzige Mittel darstellt, um Integrale mit einer dieser 
vier Funktionen im Integranden aufzulösen, so dürfte es auch einem weiteren Leserkreis willkom- 
men sein, die analogen Entwicklungen in den anderen drei Fällen kennen zu lernen. Bei der 
Herleitung dieser Entwicklungen werden wir uns wie Watson der funktionentheoretischen Me- 
“#hode bedienen, jedoch nicht ohne dabei inzwischen neu gewonnene Erkenntnisse heranzu- 
ziehen. 

Wir bemerken noch generell, daß im folgenden z, w und » beliebige komplexe Zahlen sein 


1) Die Zahlen in eckigen Klammern verweisen auf das Schrifttumsverzeichnis am Schluß der Arbeit. 
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können, von denen wir nur über » die einschränkende Voraussetzung machen wollen, daß zum 
mindesten stets v» # —1,—2, ... ist. Die Parameter x und X sollen hingegen mit x = X rein 
reelle Zahlen sein, die auf den Bereich-0 ...1 beschränkt sind. 


2. Die wesentlichsten Eigenschaften der Funktionen F,(w) und F,(w) 
in bezug auf die Veränderliche w. 


Ein wichtiges neueres Hilfsmittel, um den Charakter der Funktionen F,und F, zu erkennen, 


bilden die beiden Integraldarstellungen (2.1) und (2.2) für die in den obigen beiden Definitions- 
gleichungen auftretenden Klammerausdrücke, von denen wir bei der Gelegenheit gleichzeitig noch 
je vier andere nützliche und gleichwertige Formen angeben. Diese beiden Integraldarstellungen, 


Fa) HK) | WEL 1 KEY) 
EVA Xu) rw) IJulKu) | 2 WEN (Kw)| Sin(mr) |I,(w) J,(X%) 

+lnX (2.1): 

Jo) IJuXw)| df&r ee Er 
ey) ref: I, (w- Y1 + X? —2X :Cofp) dp 
: —InX i 

I,(w) J,(Xw) |_ EP (w)HR(Ku) IE WED (Ku) I, Jn(Ku) 
HEN Ku) wm) Alu) | ZEN WEL N, IK) 

+nX 2.2) 


J,(w) J,(Xw) 
Y,(w) Y,(Xw) 


J(w-Y1+X2—2XCojp) 
1+X?—2X6C0jo 


a 


.(1- X’Cojp)-do 


Zn BN S [XP+ XP] be Ka PB? 
TW IT 


ner 


die für beliebige reelle oder komplexe Werte von 2 und » gelten, wurden vom Verfasser bereits 
an anderer Stelle [6] hergeleitet. Sie lassen zusammen mit den beiden Umlaufsrelationen (2.3a, b) 
- sofort erkennen, daß es sich bei den beiden Funktionen F‘, (w) und F,, (w) trotz der. Mehrdeutigkeit 


J,(zw-etmri)=etmair. (zw) ./. (2.32), Julawermr)—= —etrmeir.J, (zw). . (2.3b) 


der darin vorkommenden Neumannschen oder Hankelschen Funktionen um zwei in der 
ganzen w-Ebene eindeutige Funktionen handelt, die zunächst einmal sicherlich in den Punkten 
w= 4 jy, oder j,, mit p=1,2,3... als den Nullstellen von J,(z) oder J, (z) lauter einfache 
Pole haben mit dem unendlich fernen Punkt als Häufungsstelle. Wegen der mithin wesentlich 
singulären Natur dieses Punktes sind demnach beide Funktionen ihrer Art nach meromorphe 
Funktionen. Auf ihr Verhalten in der Umgebung des unendlich fernen Punktes werden wir sogleich 
zu sprechen kommen. Überdies sind sie in bezug auf wgerade Funktionen. In der gegentei- 
ligen Behauptung bei W atso.n liegt die eigentliche Ursache für die Herleitung der fehlerhaften 
Entwicklung. 


Aus den beiden Integralstellungen lassen sich auch sehr leicht die Werte der beiden Funk- 
tionen F,, (w) und F, (w) im Nullpunkt der w-Ebene entnehmen. Es ist nämlich offenbar 


mM Se 


x : Ä 
20) 7, [PX] (v bel.).. . (2.4a), P,(0)=— [X + x] +0)... (2,4b), 
Für v= 0 wirdF, (0) = 2i/r In (X). Hingegen hat F/ (0) im Nullpunkt keinen endlichen Grenz- 
wert mehr. Tatsächlich zeigt die Gl. (1.1b) im Verein mit (2.2), daß die Funktion F/ (w) im Punkte 
w= 0 einen Pol zweiter Ordnung besitzt, in dessen Umgebung die Entwicklung (2.4b,) besteht. 
Ant 3 \ 
na Air re ET or 


Fo (w) = 
DONE 
Die Residuen der beiden betrachteten meromorphen Funktionen in ihren bezüglichen Polen be- 


stimmen sich am bequemsten unmittelbar aus ihren Definitionsgleichungen. 
Beachtet man dabei, daß z.B. 


ne ee 
Jy(jv») T)vp 
ist und auf Grund der Differentialgleichung für J,(z) die Gleichung E 
RD Er R 
Bin=—-(i-Z) nun. a 
Jvp 


j 
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. besteht, so berechnen sich für die Residuen von F,(w) und F,, (u) in. den Polen w= +j,, und 
w= + j,, die Beziehungen: 


land Fk 

Ref (Lin) + . en RT (2.7), 
Sr en uinsge 

Nei {F,(+35n)} — en en EM RR (2; 7b). 


Die noch ausstehenden Angaben über das Verhalten der beiden Funktionen F', und F, für 
. sehr große Werte von w lassen sich am bequemsten aus der dritten der fünf in den Gl]. (2.1) und 
(2.2) zusammengestellten Darstellungsmöglichkeiten für die geschweiften Klammern in den Gl. 
(1.1a,b)herleiten. Setzt man nämlich hierin für diebeidenHankelschenFunktionen und ihre 
Ableitungen ihre bekannten asymptotischen Entwicklungen für großes Argument und kleinen 
Zeigerwert ein, so entstehen nach dem Ersatz von wdurch R- eirfür F,(w) und F, (w) die beiden 
folgenden asymptotischen Ausdrücke: 


sgno Be. ar 2 i Ne: 
F,(w) m —T _,o—R-sin|p|- (X—2) +irsgn(p)- Reosp- (XK—Ü)—-iPy , , .., (2.82), 
en) aR-yeX 7 
= =) 
F, (w) en u Br .e-Rsinjpl (X a) ti-sen(p)-Recosp (X m)—ip , , .—,, (2.86). 


nR-YaX 
Sie gelten auf zwei unendlich fernen Kreisbögen in dem gegen die reelle Achse zu nicht abgeschlos- 
senen Winkelraum des ersten und vierten Quadranten der w-Ebene, und da die besagten beiden 
Funktionen gerade Funktionen sind, so haben wir damit zugleich auch Kenntnis von dem Ver- 
halten von Fund F,, längs den zu den ersten beiden spiegelbildlich zur imaginären Achse gelegenen 
Kreisbögen im zweiten und dritten Quadranten. Für ein X >x verschwinden demnach diese 
beiden Funktionen unter Ausschluß der nächsten Umgebung der reellen Achse nach einem Ex- 
ponentialgesetz, aber selbst für ein X = x findet immer noch ein Verschwinden statt, das wie 
bei 1/R vor sich geht. 

Das Verhalten, das die Funktionen F,, und F, für große Werte von win der Nachbarschaft 
der reellen Achse zeigen, untersuchen wir am zweckmäßigsten längs einer weit entfernten, zur 
imaginären Achse parallelen Geraden, die in der Mitte zwischen zwei benachbarten’ Polen dieser 
Funktionen hindurchgeht. Bezeichnen wir die Schnittpunkte dieser beiden Geraden mit u, und Ug, 
so ist demnach zu setzen: 

nv v0 


1 
w=nr(dr D+ +7 w=alat a: a 


mit v—= + iv,und ganzzahliges > -+ ©, und es wird dann nach bekannten Formeln für 
wv=u-+ti'v 


Sinn X+2)- Yeti iX-a) 
ne X: Cof ie): (u, + iv) 


Sind =XH 2 &)]-i-Coile RE A. 


-i-&oj|e(1-X x).+ = 


De ua] 


F,(w) > (2.8%), 


1 
5 5 iv ti u (X + 2)| 


Fy,(w) = (2.8). 


n:Yx X: Cof 22): (u, + iv) 


Für ol= R-sin |p)— © entstehen daraus wieder die obigen Gl. (2.8a,b), denn wegen X > x 
überwiegt dann im Zähler der letzten beiden Gleichungen stets der erste, Summand. Darüber 
hinaus unterrichten die beiden Gl. (2.8«,ß) darüber, daß auch in Richtung der reellen Achse oder 
einer Parallelen zu ihr, auf der die weit entfernten Pole von F', oder F,, liegen, diese beiden Funk- 
‚tionen zum mindesten wie 1/R abnehmen, sofern veraussetzungsgemäß dafür gesorgt. wird, daß 
bei dem Hinausrücken der betrachtete Punkt stets zwischen zwei Pole zu liegen kommt. 


3. Die Partialbruchentwicklungen für F',(z) und F, (2). 
Wir: betrachten an erster Stelle das komplexe Integral (3.1), dessen symmetrisch zur 


1 w-dw , j 
are > Mr 202 ne 
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u- und v-Achse der w-Ebene verlaufender Weg etwa oben und unten aus zwei Kreisbögen besteht, 


die links und rechts durch die im vorigen Abschnitt eingeführten, zurimaginären Achse parallelen 
Geraden + u, zu einem geschlossenen Weg ergänzt werden. Wir denken dabei von vornherein 
an den Grenzfalleines — 0, in dem der Radius der beiden Kreisbögen unendlich groß wird. Dann 
ändert der unter dem. Integralzeichen außer F (w) stehende Faktor w/(@? —w) nichts an der Tat- 
sache, daß für g, R— der in der ganzen w-Ebene eindeutige Integrand von (3.1) in allen Rich- 
tungen dieser Ebene mindestens wie 1/R verschwindet. Somit ist auch der Wert des Integrales 
selbst gleich Null. 2 

Der Integrand von (3.1) hat nun außer den schon in Abschnitt 2 erwähnten einfachen Polen 
in den Punkten w= +5, mitp=1,2,3... auch noch zwei einfache Pole in den beiden Stellen 
w= + z mit den wegen der Geradwertigkeit von F,(w) einander gleichen Residuen (3.2a). 


1 1 ? 

Re eg a 
£ AR een 

i Re ri Re Re, 


Auch in den beiden Polen w= + f, haben jetzt gemäß (3.2b) und (2.7a) die beiden zugehörigen Resi- 
duen einerlei Vorzeichen. Bei der Wahl des Integranden F,, (w)/(2? — w?), wie er bei Watson 
a.a. ©. vorkommt, würden die Residuen in’ den Polen w—= + z und w= + j,, entgegengesetzte 
Vorzeichen erhalten, und ein von Null verschiedenes Resultat könnte sich dann bei der Inte- 
gration nur einstellen, wenn F,(w) eine ungerade Funktion wäre. 

Die Anwendung des Residuensatzes führt auf Grund dieser Feststellungen unmittelbar zu 
der für alle reellen und komplexen Werte von » und z gültigen Gl. (3.3), in der die rechtsstehende 


nr ,Jy(®2), = © I, (2, )J,(X, "= 1 
4 J,(@) Wr) Wa Ze N 


Pi (jv») 2 —Jvp 


Reihe für alle Werte von x und X mit O2 x <X <1 absolut und gleichmäßig-konvergiert. 
Die daher zulässige Differentiation von (3.3) nach X liefert nach dem Übergang zu dem speziellen 
Wert X = 1 die Beziehung (3.3a), in der jetzt die rechtsstehende Reihe nur noch bedingt kon- 


J,(&2) s VITIIBELACHT) 

Ze ae a a 
J,(2) 2 J, (jun) 
vergiert. Sie ist aus diesem Grunde auch nicht mehr in dem ganzen Intervall0O= x <1gleich- 
mäßig konvergent, und damit hängt es zusammen, daß sie für x = 1 versagt, denn links entsteht 
für diesen Wert 1, während in der Reihe für = 1 jedes Glied für sich verschwindet. Die Ent- 
wicklung (3.3) ist identisch mit den a.a.O. von Sommerfeldund Carlslaw mitgeteilten 


Entwicklungen, und sie ist es, die an die Stelle der fehlerhaften Entwicklung bei Watson zu’ 


treten hat. Die aus ihr hergeleitete Partialbruchentwicklung (3.3a) ist auch bei N. Nielsen 
zu finden. 


An zweiter Stelle befassen wir uns mit dem Integral (3.4), über dessen Integrationsweg. 


die gleichen Bemerkungen gelten wie oben. Auch in diesem Falle muß das Integral wegen. des 


1 ’ w+ dw 
a Hu mit 0, (0 Te ar er (3.4) 
q 

Verhaltens von F, (w) in den unendlich fernen Teilen der w-Ebene einen verschwindenden Wert 
haben, und auch die Residuen des Integranden, für die mit F, und j,, an Stelle von F, und j, 

ebenfalls die Gl]. (3.2a, b) gelten, dürfen in den beiden Polpaaren w= +z und w= Fr wie- 
derum keinen Einfluß des Vorzeichenwechsels erkennen lassen. Die Anwendung des Residuen- 
satzes führt damit im vorliegenden Falle zu der Entwicklung (3.5). Sie ist mit Ausnahme der 
Polstellen gleichfalls für alle Werte von zund» # —1, —2, —3 ... gültigund konvergiert für alle & 


J SET EINBRUNE.S, 
s ET IT RSS EHEN a en -.- az 


2 . 2 2 ./ 
2 DL JPD EN Jy(jv2) 2? —jyp 


N NET 


und X des Bereichs 0<x<X <I1 absolut und gleichmäßig. Das erste, einzelne Glied auf der 
rechten Seite von Gl. (3.5), das nur im Falle » = 0 einen von Null verschiedenen Wert hat, ver- 
dankt sein Auftreten dem Umstande, daß dann für den Integranden von (3.4) in Rücksicht auf 
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die Gl. (2.4b,) zu den bisherigen Polen noch der Pol w = 0 als einfacher Pol mit dem Residuum 
4 ill 22) hinzutritt. 

Setzt man in (3.5) X = 1, so entsteht die einfachere Gl. (3. 5a), die der früheren Beziehung 
(3.3a) entspricht, jedoch mit dem Unterschied, daß sie ebenso wie die Ausgangsgleichung (3.5) 


J,(%2) jv „(&j»») T 
Der 2 De ı a le 
2 J,(2) dort ee J,(jv2) FURB PA 2 m 


absolut und gleichmäßig ee, In (3.5a) darf demnach «= 1 gesetzt werden. Da nun Ir die 
kleinste unter allen Wurzeln j,,ist, so darf man für ein (z) <j,, beide Seiten der Gl. (3.5a) in eine 
nach steigenden Potenzen von 2? fortschreitende Reihe entwickeln, und man gelangt so zu der 
für allev-# —1, —2, ... gültigen Darstellung (3.6). Der Vergleich gleich hoher Potenzen von zführt 


> el; 
SEN E- 1 1 Sn 2 
zu 27 2 |=5°" 9, ee Dar): 22 Een (3.6) 
2=0 - Ip=1Jvp (v9) zy +24) 9, 2 A=0 
a At £ 


dann für die ersten fünf Summen in der Een eckigen Klammer zu den Formeln (3.6a...e). 
Sie gelten natürlich, wie man an ihrem Aufbau ohne weiteres erkennt, nur für Werte von » = 0, 
—l, —2, —3 usw. Für » = 0 beginnt die im Nenner von (3.6) stehende Reihe nicht mit der Potenz 


Dee EREAEERERUN) Sn ae (3.60), 


PP 2y — jvr Be 4»2(y-+1) 
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(3.6e) 
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2°, sondern mit dem Glied —?/2. Nach der Division der beiden Reihen für J, (2) und 2: J, (z) 
heben sich dann aber in (3.6) die beiden Glieder mit der Potenz 2”? wieder Be so daß sich also 
tatsächlich, wie es die Gl. (3.6) angibt, auch im Falle » = 0 rechts eine polfreie Potenzreihen- 
entwicklung nach 2? einstellt. Es werden nur eben die Berechnungsformeln für a; (0) andere wie 
für a, (v). Es erübrigt sich jedoch, auf diesen Grenzfall noch ausführlicher einzugehen, da für 
v=0 J,() = —Jı (2) und 5,» = ji, wird. Für die Wurzeln j,, sind aber Beziehungen der Art, 
wie sie die Gl. (3.6a... e) darstellen, längst bekannt und sowohl beiWatson[1]als auch bei 
Nielsen [7] zu finden, Dagegen scheinen die obigen, für j, aufgestellten Reihen bisher nicht 
bekannt geworden zu sein. - 


4. Definition und Eigenschaften der Funktionen 6,(w) und & (w). 


Man gelangt zu einer recht nützlichen Abwandlung der soeben aufgestellten Entwick- 
lungen, wenn man an Stelle der Funktionen F,(w) und F, (w) von Gl. (1.1a, b) die durch die 
Gl. (4.1a, b) erklärten Funktionen @, (w) und @, (w) verwendet. Ihr auf den beiden äußersten 


= AA Ku) Aut) ER (0). Wa a a, 
u. {HP (Xu) -J,(w)— HD (w)-J,(Xw)}= ar Fy(w) . . . (4.16) 


Gleichungsseiten angegebener Zusammenhang mit den zuerst eingeführten Funktionen. F', (w) 
und F,, (w) läßt es als das Ratsamste erscheinen, die Eigenschaften der beiden neuen Funktionen 
aus denen der beiden früheren herzuleiten. Im Hinblick auf die beiden Gl. (2.3a,b) können 
wir auf Grund dessen sofort sagen, daß es sich auch bei den Funktionen @G,(w) und @G, (w) um 
zweiin der ganzen w-Ebene eindeutige, meromorphe Funktionen handelt, die jedoch zum Unter- 
schied gegenüber den beiden anderen Funktionen nicht gerade, sondern ungerade Funktionen 
sind. 

Ferner erkennen wir, daß man, um die Werte von @,(w) und @,(w) in der unmittel- 
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baren Nachbarschaft des Nullpunktes der w-Ebene zu erhalten, nur nötig hat, die früheren 
Gl. (2.4a, b, b,) mit dem Faktor 


v/(zw) für v»=0 oder —aw/2 fürv—=0 


zu multiplizieren. Die Funktionen @,(w) und @,(w) besitzen also nicht bloß eine sich von 
— co nach + oo hinziehende Kette einfacher Pole in den Punkten w= + jyp oder w= +» 


er 23 I, (&jvn) Io (Av )) j (4.2) 
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mit den Residuen (4.2a, b), sondern sie haben auch noch im Nullpunkt der w-Ebene einen 
einfachen Pol mit dem Residuum (4.2) oder (4.2). Für »—= 0 hat demnach die Funktion 


Ref, (0) = a: ra (MID) Re (4.20), 
RO)! —i/ ne! (PHX)...(42P), Re, (0) = ——- RAN 


G,(w) im Nullpunkt ein verschwindendes Residuum, die Funktion @, (w) hat aber dort das 
von Null verschiedene Residuum (4.2 ß,). e 

Hinsichtlich des asymptotischen Verhaltens der beiden Funktionen G, (w) ‚und G, (w) 
kann aus den’ Definitionsgleichungen sofort geschlossen werden, daß auch diese beiden 
Funktionen mit wachsenden Beträgen von w in jeder Richtung der w-Ebene verschwinden, 
sei es nun, daß dies exponentiell geschieht, wenn inw= R-exp(ip) O<|p|< x ist, oder 
daß es wie in den beiden Richtungen der reellen Achse in der Ordnung 1/R vor sich geht. 


5. Die Partialbruehentwicklungen von G,(w) und 6; (w). 
Wir betrachten in Verbindung mit den Funktionen G,(w) und @,(w) zunächst das Inte- 
gral (5.1), das über denselben Integrationsweg zu erstrecken ist wie die Integrale (3.1) und 


1 dw & 
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(3.4). Dann wird auch dieses Integral bei einem immer weiter ins Unendliche hinausrückenden 
Integrationsweg wegen des Verhaltens von @G, (w) (2? — w?) schließlich gegen den Grenzwert 0 
konvergieren. Die Notwendigkeit des Übergangs von dem früheren Faktor w/(z2? —w?) im 
Integranden, der eine ungerade Funktion von w ist, zu dem jetzigen, in w geraden Faktor 
1/(2? — w?) ist dadurch begründet, daß auch mit dem Übergang von der Funktion F,(w) zur 
Funktion G,(w) ein Wechsel von der Geradwertigkeit zur Ungeradwertigkeit verbunden ist, 
und nur dadurch kann es erreicht werden, daß auch im vorliegenden Falle gemäß den Gl. (5.2a,b) 
. die Residuen in den beiden symmetrisch zum Nullpunkt gelegenen Polpaaren w= +z und 
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w= +2 2 u) 22 3 22 R man) 

ugs AN LIE - 
Re N 
(w = + j»» wiederum das gleiche Vorzeichen bekommen. Das Residuum des Integranden von 


3.1) in dem einfachen Pol w = 0 ist aber offenbar lediglich das 22-fache des durch Gl. (4.2.6) 

angegebenen Wertes. Daraus ergibt sich dann nach dem Residuensatz sofort die folgende im 

Bereich 0Z« <X <1 gültige, in 2 mit Ausnahme der Polstellen absolut und gleichmäßig 
an 2d, (2.2) 
4 J,(z) 


2 1 
D {J,(2) Y „(X 2) —J,(X2) Y,@)}= gr ee XP) 
. (5.3) 


. 


od Jvn In? (Ivo) 2 —jyp 
konvergente Partialbruchentwicklung (5.3). Differenziert man diese Gleichung nach X, was 
wegen der gleichmäßigen Konvergenz erlaubt ist, und setzt hinterher X = 1, so entsteht die 
See a x Aula), = 5 

J,(z) SEE! Jy (vn) 2 —Iyp 
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einfachere Entwicklung (9.38), aus der ihrerseits für & = 1 die bekannte Summenformel (5.3b) 
entfließt. Behandelt man diese Gleichung in entsprechender Weise wie die Gl. (3.5a), so gelangt 


J,(2) “2 J,(z) u rj, 


man zu den schon vorhin erwähnten Summenformeln für die geradzahligen reziproken Potenzen 

von Jvp, ‚die bei Watson [1] und N. Nielsen [7] angegeben sind. 

Bei dem vierten und letzten der hier zu betrachtenden Integrale (5.4), dessen Verschwin- 

den aus denselben Gründen sicher steht wie bei den drei übrigen, gelten für die Residuen des 
Br 

ai 
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Integranden in den Polen w= +z und w= + j,, nach dem Ersatz von G,(w) und j,, durch 
G,(w) und j,, dieselben Gl. (5.2a, b) wie oben, und auch das Residuum im Pol w = ergibt 
sich wiederum durch Multiplikation der Gl. (4.28) mit 22. Demgemäß kommt im vorliegenden 
Falle die folgende, in z mit Ausnahme der Pole absolut und gleichmäßig konvergente Partial- 
bruchentwicklung zustande: Es ist für ale 0=<zx= <X <1 


ENDE {I (2) Y,(X2)—J,(X2) Y,(2)} = ec X]. (1—6gy+ don: ©?) 
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Für X = 1 nimmt sie die einfachere Gestalt der Gl. (5.5a) an. Differenziert man diese Gleichung 
(a2) Ser  _ Inlaion) 
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nach &, setzt sodann x —=1 und führt rechts den aus der Differentialgleichung für J,(2) fol- 
genden Wert für J,(j,,) ein; so gelangt man zu der Gleichung 


diInJ,;@ »—1 | 
oz Fe 22. I aa 
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deren Integration nach z zwischen den Grenzen O0 und 2 unmittelbar zu.der Produktdarstellung 
(5.6) für die Funktion J,(z) führt. Sie besitzt einen ganz ähnlichen Aufbau wie die bekannte 
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Produktdarstellung für J, (2). 
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6. Schlußbemerkungen. 

Von den vier im vorstehenden hergeleiteten Partialbruchentwicklungen (3.3), (3.5), (9.3) 
und (5.5) beziehen sich die erste und dritte auf Entwicklungen nach den Wurzeln von J,(z) 
und die zweite und vierte auf solche nach den Wurzeln von J,(2). Stellt man durch Differen- 
tiation der Gl. (3.3) und (3.5) nach dem Parameter x deren Gleichheit mit den linken Gleichungs- 
seiten von (5.3) und (5.5) her, so entstehen dadurch bei der differenzierten Gleichungsgruppe 
auf der rechten Seite in beiden Fällen zwei bedingt konvergente Reihen, während die Reihen 
der anderen Gleichungsgruppe absolut konvergieren. Da die Subtraktion der beiden Gleichungen 
innerhalb jedes Paares ohne Änderung der Reihenfolge der Glieder erfolgen kann, so ist sie er- 
laubt, und man gelangt dann zu den beiden von z unabhängigen Reihenentwicklungen (6.1) 
und (6.2), die für jedes <<‘ X innerhalb des Bereichs 0...1 gültig sind. 


3 en ori 2) RR ee in. 
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Durch andere Kombination der vier Partialbruchentwicklungen kann man auch zu Dar- 
stellungen gelangen, bei denen auf der linken Seite von Gl. (5.3) z. B. an Stelle des Quotienten 
J,(x2)/J,(z) die Quotienten J,rı (©2)/J,(2) stehen. Es zeigt sich aber, daß diese Entwick- 
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lungen nicht mehr einen so einfachen Aufbau haben wie die oben behandelten. Deswegen soll 


auch davon abgesehen werden, sie hier besonders anzuschreiben. > 
Praktisch bedeutungsvoller ist die Frage, welche Funktion von den in den Gl. (3.3), (3.3), 
(5.3) und (5.5) rechter Hand stehenden Partialbruchentwicklungen beschrieben werden, falls 


die reelle Veränderliche x, die ja in diesen Entwicklungen stets kleiner oder höchstens gleich 


X sein sollte, größer als X wird. Für die Gl. (3.3) und (3.5), deren Reihenglieder bei gleichzeitig 
absoluter Konvergenz symmetrische Funktionen von x und X sind, beantwortet sich diese 
Frage von selbst. Es braucht nämlich in diesen Gleichungen im Falle eines & > X auch linker 
Hand nur & mit X vertauscht zu werden. | E, Ren 

Um dieselbe Frage für die Entwicklungen (5.3) und (5.5) zu klären, liegt der sicherste 
und am weitesten bringende Weg darin, von den beiden in w ungeraden Funktionen (6.3a) 


Au) = HER (X) du) EN) IL} Re Se 6) 
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und (6.4a) auszugehen und sie auf die nämliche Art wie in den vorangegangenen Fällen in Par- 
tialbrüche zu zerlegen. Bei der dafür im vorstehenden bereits geleisteten Vorarbeit erfordert 
dieses Vorgehen nur geringe Mühe, und es entstehen auf diese Weise nach kurzer Rechnung 
die folgenden beiden Darstellungen: Für alle 0O=x= <X <1 ist 


reed) ,(lTz) r } 1» - 
eo ir) + gl + Rt | 


(irn I 


An Iv-(Iyp) 2 —jyp 


y% EB J, I Br | ah (6.3), 


el 
IT 2:J,(%2) 74 7 ” Y£ 1 En 
WOHIN IK} +, ZI 
© N > BEE ke 
a Jvp ‚JIlzjvp) I (Kjvo). 22 
zum nu 


Wie der Vergleich mit den in (5.3) und (5.5) auftretenden Reihen zeigt, geben sie in der Tat 
auf die oben aufgeworfene Frage unmittelbar Auskunft. Die in (6.3) und (6.4) links stehenden 
Ausdrücke kann man auch aus den Gl. (3.3) und (3.5) durch Differentiation nach X herstellen. 
Die dabei rechts entstehenden Reihen sind aber verschieden von denen der beiden letzten Glei- 
chungen und im Gegensatz zu diesen nur bedingt konvergent. Ihre Subtraktion führt zu dem 
Summenwert der für alle x <X des Bereichs 0...1 gültigen beiden Reihen (6.5) und (6.6). 


TH IE) 1 
N aa 7 ST RI NN = 
en = [ An ee cu: 
ee) Lane 
: 3 ; a I ET RR NE er 
u Rn el x ] (6.6). 


Sie ergänzen die beiden Formeln (6.1) und (6.2) dahin, daß sie die Summen der darin vorkom- 
menden Reihen auch für den Fall angeben, daß die Argumente der abgeleiteten Besselschen 
Funktionen nunmehr den größeren Parameter enthalten. 
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Mutungsgrenzen und Mutungswahrscheinlichkeit. 
(Aus der Statistischen Abteilung des W. G. Kerckhoff-Instituts.) 


Von M. P. @eppert in Bad-Nauheim. 


Eingangs werden die von Gini gegen die statistische Testmethodik erhobenen Einwände widerlegt und 
logischer Inhalt und Grenzen der statischen Zahlenprüfungen analysiert. Durch schrittweise Übertragung 
des Gedankenganges auf das Problem .der Schätzung eines Parameters in einem Kollektiv aus einer Stich- 
probe wird versucht, die Begriffe „Mutungsgrenzen‘‘ und „Mutungsverteilung‘‘ in Einklang zu bringen 
und ihren logischen Sinn und Zusammenhang zu klären. 


In the beginning, Gini’s objections against the statistical test criteria are refuted, and logical contents 
and limits of the tests of significance are amalyzed.- By transferring by steps the reasoning to ihe problem 
of estimation of a parameter in a population from a sample, the A. tries to bring the concepts ”confidence 
lIimits‘‘ and ”fidueial distribution“ to harmonize, and to explain their logical significance and coherence. 


D’abord, les objections soulevees par Gini contre les methodes de la statistigue critique sont r&futees, 
et leur logique contenu et leurs confins sont analyzes. En transferant par degres le raisonnement au problöme 
de l’estimation d’un parametre dans un collectif par un Echantillon, on essaie de meltre en accord les notions 
de »limites de confidence« et de »distribution fiduciale«, et d’Eclaircir leur logigue signification et coherence, 

B Hayare ompoBepraertca Bo3pa:keHune l:kuHu IPOTHB CTATUCTHYeCKOH MeTONUKH opo6 u 
AHAJIMSUPYIOTCH NOTHYECKO® Copep:KaHue MH TPAHHIIBI IPUMeHeHHA CTATUCTHYECKUX YUCHEHHBIX 
HCHBITAHHH. 


Ilyrem tocremeHHoro mepeHoca 3TUX coo6öpakeHnä Ha 3anayy OIECHKU Hapamerpa - 
KONIeEKTHBA TIPU TIOMOINM OTAEJIBHBIX. npoÖ AelaeTcH TIONBITKA COTNACOBATb -U YTOYHHTB 
AHOTUYeCKHÜ CMEbICH HOHATAH: „BePOATHEIE NPenellgl‘‘ U \,,BEPOATHO® PacıpepeleHne‘‘. 


Je zahlreichere und mannigfaltigere Methoden einerseits die moderne mathematische 
Statistik — speziell die englische Schule — entwickelt, je stärker andererseits diese moderne 
Methodik die praktische Statistik durchdringt und je breitere Anwendungsgebiete sie sich in 
der Wissenschaft — ganz besonders in den Naturwissenschaften —, in der Technik, Wirtschaft 
usw. erobert, desto mehr darf man vom mathematischen Statistiker fordern, daß er sich über 
Inhalt, Bedeutung, Berechtigung und Grenzen der von ihm geschaffenen oder mindestens von 
ihm empfohlenen und weiter verbreiteten Methoden ohne Vorurteile und Rücksichten restlos 
eindeutig klar werde. 

Die Grundlage vieler moderner Methoden der Statistik — wenn man von der rein be- 
schreibenden Statistik absieht, die ihre Betrachtung auf die Beobachtungsergebnisse selbst 
beschränkt und auf jede Verallgemeinerung derselben bewußt verzichtet — bildet das Prinzip 
der statistischen Zufallsprüfung oder Bedeutungsprüfung (engl. test of significance). Mit der 
Möglichkeit, auf wahrscheinlichkeitstheoretischer Grundlage empirisch gewonnene Zahlen mit 
einem bestimmten Zufallsschema zu vergleichen und hieraus Schlüsse auf das Zutreffen des 
angenommenen Zufallsschemas zu ziehen, steht und fällt die fundamentale Bedeutung der 
statistischen Methode für die Naturwissenschaft und andere Wissenschaften. 


Ehe wir zum eigentlichen Thema des Aufsatzes kommen, ist es daher erforderlich, ganz 
kurz das Wesen einer jeden statistischen Zufallsprüfung, aus jeglicher speziellen Fragestellung 
herausgeschält, darzulegen, um so mehr, als die Berechtigung derselben von Gini [3a, b], dem 
führenden Statistiker Italiens, in geistvoller Weise in Frage gestellt worden ist!). 

Wir nehmen an, das Zufallsereignis Z könne sich unter zwei sich ausschließenden und 
komplementären Voraussetzungen H und H = Nicht H verwirklichen, und zwar unter der 
Annahme H mit der Wahrscheinlichkeit 97 (E), unter der Annahme H mit der Wahrscheinlich- 
keit oz (E); analog seien 


GH) l—gn(E) bzw. rlE=1— all)... ...0.. 208 


die Wahrscheinlichkeiten für das Eintreffen des komplementären Ereignisses # = Nicht E unter 


der Voraussetzung H bzw. H. Ist ferner » die apriorische Wahrscheinlichkeit für die Gültigkeit 
der Annahme H, mithin 


ee a nr Ss 


diejenige für die Gültigkeit von H, so ist nach der Bayesschen Regel die aposteriorische Wahr- 
scheinlichkeit für die Annahme H, wenn das Ereignis E eingetroffen ist, 


EN: p*Pn(E) Be p"on(E) a (3) 
we (mM + gu) B-galE) FI—p) ga) 


‘ A) Eine eingehende Stellungnahme zu Ginis Einwänden findet sich in [2b]. 
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und diejenige für HZ, wenn E nicht eingetroffen ist, 


War ode 

p-pn(E) +P pa(E) 
Beide Werte lassen sich nur dann berechnen, wenn @y(E),yz(E) und p bekannt sind, und 
können, sofern nicht einer der Werte 9z(E), ya (E), pa(E), va (E) gleich Null ist, mit vari- 
ierendem p alle Werte zwischen 0 und 1 annehmen. Da aber im allgemeinen nur 9 (#) bekannt 
ist, während p und zumeist auch gz(E) nicht gegeben sind, zieht Gini hieraus den pessimi- 
stischen Schluß, daß ein Rückschluß vom Prüfereignis E auf die zu prüfende Hypothese H 


p- Yu (B) 


2 U-gaMm) 


wg(H) — pP Yn(#) 
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unmöglich sei und die üblichen statistischen Kriterien auf Trugschlüssen beruhen müssen. 


In Wahrheit beruht jedoch die Zufallsprüfung nicht auf der (im allgemeinen praktisch 
unmöglichen) Berechnung der aposteriorischen Wahrseheinlichkeit Gl. (3) oder (4) für H, 
sondern sie beruht lediglich auf einer zweckmäßigen und sinnvollen 
kleine positive Zahl e, 


als „Zufallsgrenze‘ festgelegt. Dann wird folgende Verabredung getroffen: Tritt das Ereignis E 


ein, obwohl 


1 
0<e<z 


Ya(E) se 


ist, so wird die Annahme H abgelehnt. Tritt E ein und ist 
Ya(E)>e. 


so ist H nicht abzulehnen, aber auch nicht erwiesen; es wird dann auf eine Aussage verzichtet. 


Entsprechend wird, wenn das Ereignis E = Nicht E eintritt, im Falle 


die Hypothese H, abgelehnt, im ‚Falle 


über HZ nichts ausgesagt. 


Selbstverständlich stellt diese Schlußweise nicht einen zwingenden logischen Schluß dar: 
denn auch bei noch so kleinem e kann E trotz Gl. (6a) auf Grund von H eintreffen. Sie ist, 
wie sowohl Neyman und Pearson [5, 6], welchen das Verdienst zukommt, das Wesen der 
Zufallsprüfungen in einwandfreier Weise begründet und formuliert zu haben, alsauch Koller ([4a], 
S. 3—4, 6-7; [4b], S. 131, 140, 206) aufs schärfste betonen, nichts weiter als eine willkürliche 
Verabredung. Jedoch ist diese Verabredung außerordentlich vernünftig, da sie es gestattet, die 


Ya(E) ze 


Yp(E)>e 


Verabredung. Es sei eine 


RR) 


2:08 


MR ee 


(72) 


>78 


Wahrscheinlichkeit, bei Anwendung der Regel eine Fehlentscheidung zu treffen, wie folgt, 


abzuschätzen. 


Tafel 1 gibt eine Übersicht über die möglichen Fälle und die entsprechenden Entschei- 
dungen gemäß der getroffenen Verabredung. Je nach der in Wahrheit zutreffenden Voraus- 
setzung H oder H und dem tatsächlich eingetroffenen Prüfereignis E oder E sind insgesamt 
vier Fälle (I—IV) möglich, deren Wahrscheinlichkeiten in der 4. Spalte von Tafel 1 verzeichnet 
sind. Von den zur Entscheidung über H herangezogenen Ungleichungen (6a) bis (7b) ist entweder 
Gl. (6a) oder (7a) oder (6b) und (7b) gleichzeitig erfüllt; diein den entsprechenden zwölf Unter- 
fällen laut Verabredung zu fällenden Entscheidungen über H sind in den Spalten 5, 7, 9 der 


Tafel 1 angegeben. 


Tafel 1. Mögliche Fälle und Entscheidungen, wenn die Voraussetzung H nach dem Eintreffen 
des Prüfereignisses Z oder E beurteilt wird. 


an tatsächlich 
wahre eine 
Voraus- 


setzung Ereignis 


Wahr- 


getroffenes | Scheinlich- 


eit 


pH(E)<Se 


PyHE)<e 


9y(E)>s Py(E)>e 


Entscheidung 
falsch 
oder 
richtig 


Verworfene 
Hypothese 


Entscheidung 
falsch 
oder 
richtig 


Verworfene 
Hypothese 


Entscheidung 
Verworfere falsch 
Hypothese oder 
richtig 


I H 

II H 
III H 
IV H 


Sr ES 


HE) H falsch 
:om(B) H richtig 
:pH\E) -_ leer 
on(#) leer 


— leer 
leer 
falsch 
richtig 


NS | 


— leer 
= leer 
= leer 
= leer 
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Aus dieser Übersicht erkennt man, daß die Wahrscheinlichkeit, eine Fehlentscheidung 
zu treffen, die sog. Irrtumswahrscheinlichkeit J(E, H; e), wenn Gl. (6a) gilt, gleich derjenigen 
des Falles I, wenn Gl. (7a) gilt, gleich derjenigen des Falles III, und wenn Gl. (6b) und (7b) 


zusammen gelten, gleich Null ist. Das Resultat läßt sich in Form eines Vierfelderschemas 
angeben: 


Tafel 2. Irrtumswa hrscheinlichkeit i 


(6a) ". leib] 
ya(E)<e | eHa(E)>e 
(Ta) | all) <e — JB, H;e)= p-pn(B)<e 
(7b) ga(B)>e J(E,H;e)=p'pu(B)<e J(E, H;e) = 


oder in Gleichungsform: 


P'Ya(E) Se für Yu(#) se 
J(B,H;e)= ip pd)se „ Yu(E) <e | ae 
0 »» Yu(E)>e, Yu(E)>e, 
und genügt stets; unabhängig vom Werte », der Ungleichung 


ES Be Se RE, en. - (9). 

Die Wahrscheinlichkeit, bei Anwendung der getroffenen Verabredung auf 
dieHypothese Hund das Prüfereignis Z oder E eine Fehlentscheidung zu treffen, 
ist demnach unter Berücksichtigung der unbekannten apriorischen Wahrscheinlichkeit » der 
Voraussetzung H stets höchstens gleich der der Verabredung zugrundegelegten 
Zufallsziffer e. 

Diese Aussage ist nicht zu verwechseln mit der folgenden leider nicht selten anzutreffenden, 
absolut falschen Formulierung der statistischen Kriterien; da trotz Gl. (6a) Z beobachtet worden 
sei, sei H mit einer höchstens e betragenden Wahrscheinlichkeit richtig, also mit einer wenigstens 
1—.e betragenden Wahrscheinlichkeit unzutreffend. Eine solche Aussage dürfte sich nur an 
die aposteriorischen Wahrscheinlichkeiten Gl. (3) oder (4) knüpfen, die, wie oben dargelegt, 
im allgemeinen nicht berechnet werden können und keinerlei angebbaren Schranken unterliegen. 
Deutlich unterscheiden sich die aposteriorischen Wahrscheinlichkeiten wa(H) bzw. wz(H) 
einerseits und die Irrtumswahrscheinlichkeit J(E, H;e) andererseits durch das Kollektiv, auf 
welches sie sich beziehen. wz (H) bzw. wz (H) bezieht sich auf die Gesamtheit der mit E bzw. E 
verträglichen Hypothesen H, H, dagegen J(E, H; e) auf die Gesamtheit aller aus dem Eintreffen 
von Z oder von E bezüglich der Gültigkeit der Voraussetzung H getroffenen Entscheidungen; 
bei wz(H) bzw. wg (H) werden in Tafell nur die Fälle I, II bzw. III, IV in Betracht gezogen, 
bei J(E, H;e) dagegen alle Fälle I—IV. 

Ginis Einwand trifft demnach nicht das Grundprinzip der statistischen Zufallsprüfungen, 
sondern lediglich die — allerdings in den Anwendungen der Statistik recht häufigen und in der 
statistischen Literatur weit verbreiteten — leichtsinnigen und gedankenlosen Formulierungen 
der statistischen Kriterien sowie die übertriebenen Schlüsse, die aus ihnen gefolgert werden. 

Dieses allen statistischen Prüfungen zugrundeliegende Prinzip wird nun folgendermaßen 
auf den Rückschluß von einer in einer Stichprobe beobachteten oder gemessenen Größe auf,die 
ihr in der Gesamtheit zugrundeliegende Größe angewandt. Eine statistische Gesamtheit @ 
besitze einen Parameter © (z.B. O=Mittelwert bzw. mittlere quadratische Abweichung einer 
statistischen Zahlenreihe G). Die Variable 7’ sei eine Schätzung desselben in Stichproben be- 
stimmten Umfanges, d.h. zur Schätzung des unbekannten Wertes werde bei zufallsmäßiger 
Auswahl einer Stichprobe © (bestimmten Umfanges) aus @ der Wert 7 benutzt, den ein als 
Funktion der Stichprobenglieder sinnvoll definierter Parameter in der betreffenden Stichprobe 
annimmt (bei den obigen Beispielen: 7’=Mittelwert bzw. mittlere quadratische Abweichung 
in der Stichprobe ©). Die Verteilung des Stichprobenparameters 7 in Stichproben gleichen 
Umfanges sei bekannt und kontinuierlich mit der Wahrscheinlichkeitsdichte 


MRION SV en ans ne 
d.h. die Wahrscheinlichkeit, bei Entnahme einer Stichprobe aus einer Gesamtheit mit dem 
Gesamtheitsparameter © einen Stichprobenparameter mit einem im Intervall 7, T+dT 
gelegenen Wert anzutreffen, sei f(T|O)dT. Die Summenwahrscheinlichkeit dieser bedingten 
Stichprobenverteilung sei 


P=F(T|O) = (fit!) ne 0 Re Be u. 
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wobei die Integration der Einfachheit halber über alle reellen 7’ erstreckt ist und f (T|O) außer- 
halb seines Definitionsbereiches gleich Null gesetzt ist. Wegen der Monotonie von F(T|O) als 
Funktion von T läßt sich GI. (11) bei festem © nach 7 auflösen.: 
- TED(0|P) 2 ee 
wo ® bei festem © und wachsendem P monoton zunimmt. 
Wir setzen ferner voraus, F(T|®) sei als Funktion von © bei festem 7 monoton, und. 
zwar — ohne Einschränkung, da man den umgekehrten Fall durch einfache Bezeichnungs- 
änderung leicht auf diesen zurückführt — monoton fallend mit wachsendem ©. Infolgedessen 


fällt auch die Umkehrfunktion 
= PTR 2a SE ee 


beifestem T und wachsendem P monoton, und mithin wächst 7 (T|P) beifestem P und wachsen- 
dem T und ®(©|P) bei festem P und steigendem © monoton. a 
Die zusätzlich vorausgesetzte Monotonie von F(T|®) als'Funktion von © ist also gleichbedeutend mit 
derjenigen von ®(®|P) als Funktion von © und besagt demnach, daß hei fester Gesamtwahrscheinlichkeit P 
zunehmenden ®-Werten steigende T-Werte entsprechen — eine Forderung, die man an eine ‚„‚Schätzung‘‘ 7’ des 
Parameters © billigerweise stellen kann. ee 

Ziehen wir nun die Möglichkeit in Betracht, daß die Gesamtheit @ ihrerseits einer Ober- 
gesamtheit derartiger Gesamtheiten mit verschiedenen @-Werten zufallsmäßig entnommen sei, 
so ist auch der. Gesamtheitsparameter © in dieser Obergesamtheit nach einem apriorischen 
Wahrscheinlichkeitsgesetz mit der Dichte (©) verteilt, so daß 

o(9) >20, [oe a0=1, 2.02. 20 . (14) 

1St, 
Wir nehmen nun an, in einer Stichprobe S aus der Gesamtheit @ sei der spezielle Para- 
meterwert T’—=t beobachtet worden, und-wir wollen untersuchen, ob der unbekannte Parameter 
@= ® von @ ober- oder unterhalb einer bestimmteu Grenze ®, liege. 

Dann können wir zunächst als zu prüfende Hypothese die Ungleichung 


& A ee 
als dazu komplementäre 

d<d, : Veen (H»,) 
auffassen. Als Prüfereignis wählen wir die Gültigkeit der Ungleichungen: 

tn. a SS 

Te 
wobei i, vorläufig eine beliebige, aber feste Zahl bedeute. 
Tafel 3. 


4’ [65 
ey(E<e Yy(E)<e en(E)>e, gy(E)>e 
Es gibt mindestens ein 
9,’ und @,’in 
-Fit|0) Ze Ftl@)21—s | 9728, 0,72%, 
füralle 02%, füralle >98, für welches 


F(t|9,’)>e, 
F(t,|93’) <1—e 


Ver- 


Ver- Ent- Ver- Ent- Entscheidung 


Wahre Tatsächlich 


ectmoftenen Ä ER N worfene scheidung | worfene| scheidung |worfene falsch 
Voraussetzung etenie Wahrscheinlichkeit Hypo- | falsch oder | Hypo- | falsch oder | Hypo- oder 
these richtig | these richtig these richtig 


oo 


IH, (929) EZ, (t<t,) [ot9-F(|O)do falsch 


TH, (#<d,)E, (St) 


richtig 


ut) „ (929%) 8, (>t,) | fo(0)-1—Fit|O)]d6 


falsch — 


® 


Je(9-L—Ftt|9)1d6 


IVA, (#<9,) E, (t>t,) richtig | — 


Zu a De 


up # , N 2 e 1 “ ’ + 
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Bei Aufstellung der der Tafel 1 entsprechenden Tafel 3’ ist zu beachten, daß das Ereignis 
(H,,) sich auf verschiedene Weisen (verschiedene @-Werte) verwirklichen kann, denen jeweils 
andere bedingte Wahrscheinlichkeiten für das Ereignis Z entsprechen; es ist mithin bei der 
Berechnung der Wahrscheinlichkeiten der vier verschiedenen Fälle I-IV außer dem Multi- 
plikationssatz auch der Additionssatz der Wahrscheinlichkeitsrechnung: zu verwenden. 


\ 


Der Fall S | : 
Fit ZSex- für. ale 2m ra, Be) 
liegt wegen der Monotonie von F(t,|O) dann und nur dann vor, wenn 
Pd) SE. N re er 
ist. Der Fall 
1—F(,)0) ze Tür alle; Oo 20: 2, ER ee 


kann, wenn die Funktion. F(7|©) für jedes feste 7’ mit wachsendem © bis Null absinkt: 


Min. F(T|0) = 0, 
[c) 


überhaupt nicht eintreten; liegt für alle T’-Werte der entsprechende (in dem dem betreffenden 
T-Wert entsprechenden größten O-Wert erreichte) Minimalwert von F(T7|®) unterhalb einer 
von 7 unabhängigen, unterhalb 1 gelegenen Schranke 


Min F(T|9)> 0'<1, 
(2) 


so ist der Fall (B’) bei genügend kleinem e(e <1—o’) ebenfalls unmöglich; nur, wenn es zu 
jedem noch so nahe an 1 gelegenen 0’<1 stets einen 7',-Wert gibt, für welchen der minimale 
F-Wert o’ übertrifft: 
Min F(T,|O) >e', 
©) 


kann der Fall (B’) bei genügend kleinem i, überhaupt eintreffen. 
Man kann jedoch ebenso umgekehrt die Ungleichung 


® = E 2 . . . . . . . er re De, (Hs) 
Tarelrs“. 
> 2 2 
a) se RE re) >&,95(M>e 


Es gibt mindestens ein 
0,’’und ®,’’in 


F(t,|0)zZ1-—e Fit |OQ)se 9. 0 9'<S%, 
füralle O< 9, für alle  @< 9, für welches 
F(,|0,’) <1—e, 
PA] 95) >8 
en; Ver- Ent- Ver- Ent.- Ver- | Entscheidung 
Wahre Tatsächlich inli ä worfene| scheidung |worfene| scheidung | worfene falsch 
Voraussetzung | eingetroffenes Wahrscheinlichkeit Hypo- |falsch oder | Hypo- |falsch oder | Hypo- oder 
Ereignis these richtig. | these richtig these richtig 
oo 
IH, (929) Z, (tt) 4 ©(9)-F(t,|O)d® — = H,, | richtig | — a: 
0 
® 
II, (d<3) 2, (sth) | [o(9)-Fltlo)dO _ = H, | falsch | — 3 


|, (# >88, (1>t) | [o(9)-TI—FitlO)]d0] H, richtig | — 
d, 


Va, #<3,)l#, (>t) Foto -raond H,, | falsch | 


| 1 Bl 
als zu prüfende, (4,,) als die komplementäre Hypothese auffassen und die Gültigkeit der 


Ungleichung x 
ee ee 


Er, 


55 MS Zah " 

az er rg f 

r eh * E 
Et Lan 3 2 
f Re 

e 


als Prüfereignis benutzen. Der Tafel 1 Bene BE Tafel & 3”, wobei v 
1— F(t,|0) se | für ale © sh 
ns ER nur dann eintritt, wenn | 


f 


ist, und Fall "1 höchstens dann — und zwar bei hinreichend großem t,— Sintrokfen PUR | 


En ' 


ro) ehr ER a 


wenn es zu jedem noch so kleinen o””>0 stets einen Wert e gibt, für den 


- 


ist.- 


Max F(T|O)< 0” 
o 


Da demnach bei der ersten Festsetzung der Fall B’, bei der zweiten der Fall B’’ praktisch E 5 
bedeutungslos und schwer nachzuprüfen ist, werden die beiden Verabredungen in eine verr 


einfachte zusammengefaßt: Ist <t, und gilt die Ungleichung 


so wird die Hypothese 


"abgelehnt, also 


gefolgert; ist t 
so wird 


abgelehnt und 


>t, und gilt 


gefolgert. Gelten, dagegen gleichzeitig die beiden Ungleichungen 


Bil ee RE EN. (15a’), 3 
ET - -(H,) | = 
ER Rn... (16°) | 
Tri) (15a”), | 
DEI a : RE A RE RE 
Sr EN ee ENT 
Erd) Ir re IR (15b), 


so wird auf jegliche Aussage verzichtet. 
Die verschiedenen Möglichkeiten sind in Tafel 3 zusammengestellt. Die Irrtumswahr- 


» Wahre 
Voraussetzung 


1H,, 8%) 
II H,, (d<%,) 


IA, (#2%,) 


IV H,, (<H,) 


scheinlichkeit J(t,]9,]e), 


J(6,j.|e) = 


Tafel 3. 


(153°) (15a’’) 
F (to] SE 
ER Ver Ent Ver- Ent- Ver 
Tatsächlich i i 
ne 2 worfene| scheidung |worfene| scheidung |worfene 
ne Wahrscheinlichkeit Hypo- | falsch oder | Hypo- | falsch oder | Hypo- 
these richtig these richtig these 


(st,) 


fo(9-F@lo)ae |'H,, | falsch 
% 


F(t,|0) d@ H,, | richtig 


d.h. die Wahrscheinlichkeit, bei Befol 
Fehlentscheidung zu treffen, ergibt sich hieraus zu 


richtig 


falsch 


For) dO<e für Fitd)<e 


d, 


— 


0 


So(0) 1 —F 


AI AOSEe „ 1-—-Flid)<se 


” Ee< F (&o)90) < 1—;, 


Entscheidung 
falsch 
oder 
richtig 


gung dieser Verabredung eine 


. (17) 


Dun Ze U, si 
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® 
oder in Tabellenform geschrieben 
Tafel 4. 
| 1 — F(t,|®,) se | e 1 — F(t,|9)) > € 
Fit) Se — / Ktoldole) = 1 (9) -Fit,)O)dO <Se 
Pl 


d ; 
Fiod)>e | Foldole=fo(9)-U—FM)]A9<e | Ittuldule) = 0 


und genügt offenbar stets der Ungleichung : 
Jod) SER (18). 


Während zwar die Irrtumswahrscheinlichkeit J für jede beliebige apriorische Verteilung (©) 
der Ungleichung (18) genügt, hängt natürlich der genaue Wert von J nach Gl. (17) wesentlich 
von der Gestalt der Funktion &(©) ab. 


Zunächst darf die dem Kriterium (£,-9,-Kriterium) zugrundeliegende Grenze ?, nicht 
von der Beobachtung T=t abhängig gewählt werden, sondern müßte eigentlich eine feste 
Größe sein. Andererseits steht es einem frei, auf eine Beobachtung 7’=t nacheinander dasselbe 
Kriterium mit verschiedenen i,-Werten anzuwenden, vorausgesetzt, daß sich die Ergebnisse 
nicht widersprechen. Eine eingehende Analyse aller möglichen Fälle zeigt, daß, wenn die Grenze ?, 
von oben oder unten her dem Beobachtungswert 7’=i genähert wird, die von den einzelnen 
t,.Kriterien gelieferten Aussagen sich niemals widersprechen, aber immer mehr verschärft 
werden, d.h. daß bei.der Verschiebung von ?, auf it zu zwar entschiedene Aussagen über die 
Lage des unbekannten # zur Grenze 9, neu gewonnen werden, aber nie verloren gehen können. 
Da man andererseits, wenn mehrere verschiedene Kriterien zur Verfügung stehen, berechtigt 
ist, das wirksamste derselben anzuwenden, d.h. dasjenige, welches möglichst selten zu leere 
und möglichst häufig zu entschiedenen Antworten führt, setzt man demnach, wie es allgemein 
in der Literatur üblich ist, 2,—=? und trifft also die vereinfachte Verabredung: 


Gilt 
PtlD or Ser ee u ae ae rer (19a’), 
so wird 
<< Do Pe (16°) 
gefolgert; gilt 
DEE IO e. e ml (19a) 
so wird 
Be Be (16°) 
gefolgert; gelten gleichzeitig die beiden Ungleichungen 
ES Er. n , (19b), 


so wird auf eine Aussage verzichtet. 

Die Bestimmung der Irrtumswahrscheinlichkeit J(£|d,|e) dieses d,-Kriteriums stützt sich 
folgendermaßen auf die für das t,-9,-Kriterium geltende Gl. (17) oder Tafel 4. Zur Ver- 
anschaulichung seien in Bild 1 7 und © als rechtwinklige Koordinaten in einer Ebene ein- 
gezeichnet, über der wir uns die Funktion »(©)-f(T|©) ausgespannt denken. In ihr seien 
ferner die beiden gemäß der Monotonievoraussetzung ansteigenden Kurven 


G=.2(71:) bzw.r T =D (Glen 2 are (&,) 
d 
— = PTH—s)-bzu ss raB(Oled ne oe ee 


eingezeichnet. Die Horizontale =, schneidet €, und €, allgemein in den beiden Punkten 
P,, P, mit den Abszissen D(d,|e) und D(d9,|1—e). Liegt nun der Punkt P,(io; %,) links ober- 
halb der Kurve €,, so ist die Irrtumswahrscheinlichkeit J(t,|d,|e) des tu-Öu-Kriteriums 
Gl. (17) gleich dem Integral der Funktion &(©) -f(7T|©), erstreckt über das von den durch ER 
hindurchgehenden Achsenparallelen rechts und unten begrenzte Gebiet; liegt P,rechts unterhalb 
von @,, so ist. das Integral über das links und oben von den durch P, gehenden Achsenparallelen 
begrenzte Gebiet zu erstrecken; liegt P, zwischen €, und &,, so ist J (to|dole) =0, aber das 
{5 9,-Kriterium aussagenleer. Beim %,-Kriterium identifiziert man nun &, mit ? und die ein- 
zelnen _Integrationsgebiete superponieren sich zu einem einheitlichen Gebiet, zusammengesetzt 
aus dem unten und rechts von den P, durchsetzenden Achsenparallelen begrenzten Gebiet und 


17* 
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dem links und oben von den P, durchsetzenden Achsenparallelen beg 
Irrtumswahrscheinlichkeit des d,-Kriteriums ist also SE 


R ee 
BICHE | o(6) -F(60,]9)A0+J w(9) -U-F@@— Le 
© Disje) - 


ee © \ 
= Joa -saiaaraoy+) „ja srlaarae 


Wir gehen nun einen Schritt weiter und lassen auch die Grenze d, varlieren, d. h. wir 
wählen unter allen verschiedenen 9,-Aussagen, die durch Anwendung aller möglichen 9... 
Kriterien auf einen beobachteten i-Wert zu erreichen sind, jeweils die entschiedenste aus. 


ß 7, z = 


dı 


A 


“ 
0, 


Bild 1. 


Infolge der oben erwähnten Monotonieeigenschaften gibt es für jedes beobachtete T=t 
zwei lediglich durch die Form der Funktion F(7|®) bzw. f(T|®) bestimmte, dagegen von der 
Verteilung &(©) des Parameters © völlig unabhängige Grenzen, die sog. Vertrauens- oder 
Mutungsgrenzen: 


%,= Pille) <= Br (21), 

für welche die Gesamtwahrscheinlichkeit des beobachteten Prüfereignisses gerade e bzw. I—e 
beträgt: = 

F(t|0,) = 8, 1 — F(t|d,) a an a (22), 


und deren Bildpunkte inBild1 auf den Kurven @&, und @, liegen. d, ist dann der größte Wert, 
für den das ,-Kriterium in der Form Gl. (19a’) gerade noch wirksam ist, 9, der kleinste Wert, 
für den es in der Form Gl. (19a’’) gerade noch wirksam ist. Doppelte Anwendung des ®,-Krite- 
riums mit 9, — d, und 9, = d, liefert also für den unbekannten Parameterwert ® die Ungleichung 


Yil—)=i<d<dh- Pirna (23). 


Die Irrtumswahrscheinlichkeit derselben ergibt sich wieder durch Superponieren der allen 
möglichen ®,-Werten entsprechenden Integrale Gl, (20), d.h. anschaulich dadurch, daß man 
‚in Bild. 1 die beiden rechten Winkel entsprechend dem variirenden %, mit‘den Scheiteln P,, P, 
längs der Kurven @,,&, entlang gleiten läßt, zur Summe der beiden Integrale erstreckt über 
das gesamte links oben von der Kurve &, und über das gesamte rechts unten von €, gelegene 
Gebiet, welche gemäß der Definition Gl. (12) oder (21) jeweils genau gleich & sind: 


J()=/ ©(0) :F(6(0)9)|9) 40 +[ »() - [1 —F(8(9]1—2)|9)] 40 = 
Pa . (24). 
= [»(O)S(T|O)ATAo +] „ge TINATdO=e+e=2e 


Durch wiederholte Anwendung eines Zufallskriteriums gelangt man also auf diese Weise 
ohne einengende Voraussetzungen zuzweinurvonderbedingten Verteilung /(T|o) 
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abhängigen Grenzen 9,<%,, innerhalb deren man den wahren Parameterwert # 
vermuten darf. Die Wahrscheinlichkeit, mit einer solchen Entscheidung fehl- 
zugehen, bezogen auf die allen möglichen #-Werten und beobachteten i-Werten 
entsprechenden Entscheidungen, beträgthierbeigenau2e. Die einzige unentbehrliche 
Voraussetzung für dieses Resultat bildet die Monotonie von F(T |) in bezug auf ©, die jedoch 
zumeist erfülltist. Dagegen ist die Kontinuität der bedingten Verteilung f(7 |) nicht unbedingt 
erforderlich für Begriff und Existenz der Mutungsgrenzen. Ist nämlich die bedingte T-Verteilung 
diskret, also F(T|©) als Funktion von 7 unstetig, so existieren die Umkehrungen Gl. (12) und (13) 
von Gl. (11) nicht für alle, sondern nur für einzelne P-Werte; als Mutungsgrenzen gelten dann 
die Funktionen Y(t|1—e,), Y(t]e)) mit geeigneten oberhalb von 1-—-e bzw. unterhalb von & 
gelegenen Argumenten, so daß an die Stelle der exakten Gl. (22) die Ungleichungen 


Pula)<e. "Lore (22%) 


treten und die Irrtumswahrscheinlichkeit der Aussage Gl. (23) nach Gl. (24) nicht genau gleich 2e, 
sondern nur höchstens gleich 2e ist: 


Meder ee (240). 


Die Kontinuität der Anfangsverteilung &(®©) ist hingegen für den ganzen Gedankengang und 
sein Ergebnis völlig belanglos. 

Dasselbe Ziel erreicht R. A. Fisher [1a] auf verblüffend einfache Weise, indem er von 
der Tatsache, daß die über die Gebiete links oben von €, und rechts unten von’@, erstreckten 
Integrale von der Form der Anfangsverteilung (©) völlig unabhängig und genau gleich’ & sind, 
also von den direkt aus der Definition Gl. (12) der Umkehrfunktion ® folgenden Gleichungen 


} ®.9je) co Eon) 
a) [ [o(@)f(T|oydT d0= J w(9): F(6(9))|9) AO = e- | w(0)d0=e, 
ja ge Er 2 2 (24) 
b) I [ea strIHara0 = i-Feaal—e|M4E=T—(1-e]/w10)40=e 
—& O|1-e) —_o — © 


( 

ausgeht?). Hiernach ist die Wahrscheinlichkeit, bezogen auf alle ©- und 7-Werte, dafür, daß 
der dem wahren ©- und beobachteten T-Wert entsprechende Bildpunkt links oben von &, 
bzw. rechts unten von €, liege, daß also zwischen © und 7 die das betreffende Gebiet kenn- 
zeichnende Ungleichung 


a) F<6(0e) bzw... b) T>DlAl ev om (25) 
oder die damit gleichbedeutende 
aut bzw: bh)  SHZU Den 239 


gelte, genau e. 

Wiederum ist also die Irrtumswahrscheinlichkeit der Aussage Gl. (23) genau gleich 2e, 
wenn F(T}9) stetig in bezug auf T und monoton in bezug auf © ist. Für die Gültigkeit des 
Ergebnisses und die erforderlichen Voraussetzungen gelten wiederum die obigen Erwägungen. 

R.A. Fisher begnügt sich jedoch nicht mit der Definition der Mutungsgrenzen Gl. (21) 
und der Angabe ihrer. Irrtumswahrscheinlichkeit, sondern. geht einen wesentlichen Schritt 
weiter. Da die Irrtumswahrscheinlichkeit der beiden Ungleichungen (25*) jeweils genau gleich e, 
also der Sicherheitsgrad derselben 1—e und von ®(©) unabhängig ist, kann man das ursprüng- 
lich konstante P= s variieren lassen und die Mutungsgrenzen für verschiedene P-Werte be- 
stimmen. Durch die eindeutige Zuordnung der Mutungsgrenzen (25*) zu P=s entsteht dann 
formal eine ‚„‚Rückschlußverteilung‘‘ oder ‚‚Vertrauens-“ bzw. „Mutungswahrscheinlichkeits- 
verteilung‘“ von © mit der Wahrscheinlichkeitsdichte m (©]|T) und der Summenwahrscheinlichkeit 

@ 


Bone neDde.n..... . ee (26). 


—o 


Um diese zu bestimmen, müssen wir in Gl. (24) die Reihenfolge der Integrationen vertauschen: 


a) f fewm-sırioaoar =. 
Rz ne 
p) f Solo) fiTIY)daoaT=e 


2) Die hier gegebene ausführliche und. bewußt umständliche Darstellung bezweckte nur, innere Übereinstimmung 
bzw. Zusammenhang desneuen Fisherschen Gedankenganges mit dem älteren Prinzip der statistischen Zufallsprüfungen 
deutlich zu machen, 
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sowohl (24a*) als auch (24b*) führen dann zur Gleichung an 
BON=1—-P=e1-F(T0 0,0 


aus der sich durch partielle Integration 
LESE? 


ie Te 


EI en a TEEN u ee‘ 


ergibt. : 
Auf diese Weise erklärt Fisher formal eine kontinuierliche Mutungswahrscheinlichkeits- 
verteilung des Gesamtheitsparameters © mit der Wahrscheinlichkeitsdichte mw(©|7) und der 
Summenwahrscheinlichkeit ®(0©|T), die von der apriorischen ©-Verteilung o(9) völlig un- 
abhängig ist. So einfach sie rechnerisch-formal herzustellen ist, so schwierig ist es, ihre wahre 
Bedeutung zu begreifen. Natürlich ist W(O|t) nicht etwa die Wahrscheinlichkeit dafür, daß, 
wenn in der Stichprobe der Parameterwert T=t beobachtet worden ist, der Gesamtheitspara- 
meter © unterhalb ® liege;-denn diese-Wahrscheinlichkeit ist nichts anderes als das Intregral 
der Bayesschen Rückschlußwahrschinliechkeit, 


WERBEN VE 


[»(9) -f«|0) 40 
a - (28), i 
.[o(0)-fu]0) a0 


und hängt wesentlich von der Anfangsverteilung (©) ab. ®(9,T) ist vielmehr die Wahr- 
scheinlichkeit, bezogen auf alle ©- und 7T-Werte, dafür, daß, wenn die Gesamtheit irgendeinen 
Parameterwert © aufweist, der Parameterwert “7 einer Stichprobe aus derselben oberhalb von 
D(O1—W (dt), also speziell für O9=% oberhalb von t liege; oder umgekehrt die auf alle ©- 
und 7-Werte bezogene Wahrscheinlichkeit, mit der Aussage, bei beobachtetem Stichproben- 
parameter 7’ liege der wahre Parameterwert © unterhalb von Y(T11— (9]i)), also speziell 
für T7’=t unterhalb von d, eine richtige Entscheidung zu treffen. Folgerichtig bedeutet dem+ 
nach die Differenz 


80,9 8A) = Fu) — Fee)... >, 


die Wahrscheinlichkeit, bezogen auf alle ©- und T-Werte, mit welcher bei 
irgendeinem beobachteten 7 die Aussage, das wahre © liege in den Grenzen 


FTI—- Ro) <oO<Y(T1 — (0,1), 
also bei 7T=ti die Aussage 
9,<9 < 0 


richtig ist. Während sich die Bayessche Rückschlußverteilung Wz(#]i) auf die Betrachtung 
der Stichproben mit einem bestimmten: Stichprobenparameterwert T’=t beschränkt, bezieht 
sich also die Fishersche Verteilung ® (d]t) auf das Kollektiv aller möglichen £-Werte. 
Zweifellos ist die Fishersche Herleitung einer von der apriorischen ©-Verteilung »(®) 
unabhängigen Rückschlußverteilung Gl. (27) bzw. (27) außerordentlich. bestechend, scheint 
sie doch mit einem Schlage den uralten Streit um die Bayessche Rückschlußweise beizulegen 
und auf die einfachste Weise den alten Traum der Statistiker zu verwirklichen, den Bayes 
selbst wegen der bei seinem Vorgehen notwendigen allgemein nicht haltbaren Voraussetzung 
der Gleichverteilung (»(©) = const.) nicht erfüllen konnte. Trotzdem neige ich dazu, in dem 
Fisherschen Vorgehen doch eine mehr formal befriedigende Lösung des Problems zu erblicken, 
der Rückschlußverteilung Gl. (27°) eine greifbare reale Existenz abzusprechen und ihr vorzugs- 
weise formalen Wert zuzuerkennen®). Denn tatsächlich liefert eben ® (dt), wie oben aus- 
einandergesetzt wurde, keineswegs eine Wahrscheinlichkeit dafür, daß der dem beobachteten t 
zugrundeliegende wahre O-Wert unterhalb von d liege — eine solche müßte ja notwendig von 
der Anfangsverteilung (©) abhängen — sondern eine erheblich schwieriger zu erklärende 
Wahrscheinlichkeit, die zwar einerseits vom beobachteten T=t abhängt, aber andererseits 
gleichzeitig alle anderen 7’-Werte mit berücksichtigt. Die Mutungswahrscheinlichkeits- 
verteilung Gl. (27) ist praktisch gleichbedeutend mit der Angabe der Mutungs- 
grenzen Gl.(21) als Funktionen der willkürlich veränderlichen Zufallsgrenze 
P=e und setzt Stetigkeit von F(T|®) als Funktion von 7 und Monotonie von F(T|®) in 


») vgl, [2a] 8. 12, | 
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bezug auf © voraus. Während bei festem & die Konstruktion der Mutungsgrenzen auch im Falle 
eines als Funktion von 7 unstetigen F(7'|©) möglich ist, dann aber die Ungleichung (23) nicht 
mit der exakten Irrtumswahrscheinlichkeit 2 &, sondern nur mit einer höchstens 2 e betragenden 
Irrtumswahrscheinlichkeit erfüllt ist, kann bei der Herleitung der Rückschlußverteilung ® (@|7) 
auf die Kontinuität der bedingten Verteilung f(7|0)dT nicht allgemein verzichtet werden. 
Jedoch ist es W. Schäfer [7]für den Spezialfall des Rückschlusses von der-in einer n-gliedrigen 
Stichprobe beobachteten Merkmalshäufigkeit 7 auf die der Gesamtheit innewohnende Merk- 


 hsuinlglt © gelungen, trotz der Unstetigkeit der bedingten Summenwahrschein- 
ichkeit 


[RT] 


Fa Deu a En 29) 


ohne besondere Voraussetzungen eine näherungsweise geltende Rückschlußverteilung ®(©]|7T) 
zu gewinhen und so wenigstens in diesem — ältesten und vielleicht wichtigsten — Spezialfall 
in gewissem Sinne den Fisherschen Gedankengang auf diskrete Verteilungen f(7|©) aus- 
zudehnen‘®). 


Schließlich setzt aber der Begriff der Rückschlußverteilung ®(O|7T) — im Gegensatz 
zu dem der Mutungsgrenzen — ebenso wie die Bayessche Rückschlußverteilung Gl. (28) die 
Vorstellung des doppelten Bayesschen ‚Urnenschemas, d.h. das Vorhandensein einer aprio- 
rischen ©-Verteilung &(©) voraus, obwohl deren Form für das Resultat gleichgültig ist, und 
erfordert, daes an die Gl. (24) anknüpft, die Kontinuität von (©), denn für eine in eine diskrete 
Verteilung entartende, oder sogar auf einen einzigen ©-Wert zusammenschrumpfende Ver- 
teilung (©) ist die Vertauschung der Integrationen in Gl.(24) unsinnig. Begriff und Geltung 
der Mutungsgrenzen Gl. (23) bestehen unabhängig davon, ob man von der — manchmal allein 
berechtigten — Vorstellung einer isolierten Gesamtheit ausgeht, deren einziger, unbekannter 
O-Wert zu ermitteln sei, von der Vorstellung also, es existiere keine Anfangsverteilung »(©), 
oder aber von der Bayesschen Vorstellung einer höheren Gesamtheit von Gesamtheiten, deren 
O-Werte einer Anfangsverteilung (©) unterliegen. Der Begriff der Fisherschen Rückschluß- 
wahrscheinlichkeit ®(©|7) baut dagegen ebenso wie derjenige der Bayesschen Rückschluß- 
wahrscheinlichkeit W3(©|T) wesentlich auf dieser letzten Vorstellung auf und verliert andern- 
falls seinen Sinn. 
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als Axiom zugrunde legte (vgl. [2a]). 
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Zusammenfassender Bericht... an 


Über die Charakteristikenverfahren der Hydrodynamik. 5 ie : 


Teil II. | 
Von Kl. Oswatitsch in Graz. 


6. Anwendung des Verfahrens auf die achsensymmetrische adiabatischoe 
Überschallströmung. Be 
Führen wir anstatt der Geschwindigkeitskomponenten % und » als abhängige Veränder- 
liche den Geschwindigkeitsbetrag w und die Geschwindigkeitsrichtung d ein, so nehmen die Glei- 
chungen für die Kontinuität und für die Drehungsfreiheit bei der achsensymmetrischen adia- 
batischen stationären Überschallströmung folgende Gestalt an (x Achsenrichtung, y Richtung 
normal zur Achse): N 


W) e0sd “N sind: in dd 9-9,” sind 
(1%) eos 0, +15) sin w,— wsin tr weoös | Fe (28). 


sin d w, — cosd-w, + wcos9-d, + wsind-d, —=0 


Die Schallgeschwindigkeit c ist hierbei bekanntlich nur Funktion des Geschwindigkeits- 
betrages w. Die abhängigen Veränderlichen g und A aus Gl. (1) sind durch w und # ersetzt. 


rn mi; B=0. 
Y 


Für sehr große Werte von %, also für große Abstände von der Achse, erhalten wir ein homo- 
genes Gleichungssystem. Das Problem der achsensymmetrischen Strömung geht hier in jenes 
' der ebenen Strömung über. Die wesentlichen Eigenschaften dieser Strömung setzen wir im fol- 
genden als bekannt voraus. Die Charakteristiken des homogenen Gleichungssystems (der ebenen 
Strömung) & = konst., 7 = konst. erhalten wir durch Integration der Gl. (15). Wählen wir im 
Zustandsdiagramm als rechtwinklige kartesische Koordinaten nicht w und Ö sondern u und », 
so erhalten wir für die Charakteristiken des homogenen Systems die bekannten Epizykloiden. 
Denken wir uns also in Gl. (23) m durch «, n durch », g durch w und h durch # ersetzt, so müssen 
bekanntlich die Orthogonalitätsbedingungen erfüllt werden. Diese haben folgende Gestalt: 


2 2 2 
— "zwsin Dos 8: u0+ (1% Jwo+ (15 sin?) wv„—0; 


en we, \ 
— ut — "7 wsin D: cos 0+us+ [1 sin? 0) wns=0. 


Da die Lösung des ebenen Problems im allgemeinen als Felderverfahren bekannt ist (Prandtl- 


Busemann- Verfahren), wollen wir das Gitterpunktverfahren der ebenen adiabatischen 
Strömung: kurz erläutern. 


Bild 7. Bestimmung des Ortes 
von P, bei ebener adiabatischer 
Strömung. 


(—.—.— Verfahren 
erster Ordnung, 

Verfahren 
zweiter Ordnung.) 


X 


Wir denken uns die Strömungsrichtungen und die Geschwindigkeiten in zwei bekannten 
Punkten der Strömungsebene P, und P, gegeben (Bild 7). Die Richtung ist hierbei durch die 
strichpunktierte Linie gekennzeichnet. Wir können dann die Punkte P, und P, im Zustands- 
raum, als dessen Koordinaten wir u und v wählen (Hodograph), aufsuchen. Durch jeden Punkt 
gehen zwei Charakteristiken. In der Strömungsebene bezeichnet man diese als Machsche 
Linien; das sind jene Wellen, welche von dem beträchteten Punkt ausgehen würden, wenn dort 


= 
j 
4 


Faden A 
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eine kleine Störung, etwa in Form einer in Strömungsrichtung gestellten scharfen Schneide vor- 
handen sein würde. Die Machschen Linien in der Strömungsebene liegen symmetrisch zur 
Stromlinie und schließen mit dieser den Machschen Winkel & ein. Die Charakteristiken der 
Zustandsebene, die ja die Bilder der Machschen Linien sind, liegen symmetrisch zum Ge- 
schwindigkeitsvektor, der mit der zugehörigen Stromlinie stets gleiche Richtung hat. Die Zu- 
ordnung der Charakteristiken'ist nun so, daß’ die von links nach rechts laufenden Charakteristiken 
in der Strömungsebene und in der Zustandsebene einander entsprechen, wenn in beiden Fällen 
in Strömungsrichtung geblickt wird. Wir nennen sie die rechtsläufigen Charakteristiken und ent- 
sprechend dann die andere Schar die linksläufigen Charakteristiken. Um den Punkt P, im Hodo- 
graphen zu erhalten, müssen wir vom Punkte P, ausgehend längs einer rechtsläufigen und vom 
Punkte P, ausgehend längs einer linksläufigen Charakteristik fortschreiten. Der Schnittpunkt 
ist innerhalb der Zeichengenauigkeit exakt P,. Um nun noch den Ort von P,in der Strömungs- 
ebene kennenzulernen, nehmen wir beim Verfahren erster Ordnung im Punkte P, der Zustands- 
ebene die Richtung der Charakteristik, welche nicht nach P, führt. Die Normalrichtung dazu 
gibt uns die Richtung der Machschen Linie. Entsprechendes gilt für den Punkt P,. Das Ver-, 
fahren zweiter Ordnung unterscheidet sich nur durch eine etwas genauere Ortsbestimmung des 
Punktes P,, indem nämlich nicht die Richtung in den Punkten P, und P, des Zustandsdiagramms 
genommen wird, sondern in zwei Zwischenpunkten P,, von P, und P, und P,, von P, und P;. 
Bei diesem Gitterpunktverfahren der ebenen Strömung ist keineswegs angenommen, daß die 
Werte der Zustandsgrößen in den ganzen Feldern konstant sind. Man verschafft sich an irgend- 
einem Ort des Zustandsraumes die Werte am besten durch Interpolation. 

Betrachten wir nun das Charakteristikenverfahren für die achsensymmetrische Strömung, 
so ist also das entsprechende homogene Verfahren die eben besprochene Methode für ebene Strö- 
mung. Die Gl. (26a) und (26b) für die Änderung der Zustandsgrößen längs Charakteristiken in 
der Strömungsebene lauten, wenn wir ® und w als Zustandsgrößen wählen: 


Br a at 
D-r+ VE)” Asindsina=0 ..... (29a), 
1 w2 a: > 83 — 9 . . 
Re een Ve) Be 2 


(29a) gilt für die rechtsläufige, (29b) für die linksläufige Mach sche Linie. Ferner it „— 9% = 
”—ys=&a (Machscher Winkel) gesetzt. Das ebene Problem erhalten wir einfach, wenn die 
Strecken s,— s; und s;— s, verschwindend klein gegenüber % sind, also wenn die Abstände groß 
sind gegenüber den Abständen auf den Machschen Linien. Wir sehen außerdem, daß der Ein- 
fluß des letzten Gliedes klein wird, wenn die Strömung bei endlichem y in Achsenrichtung verläuft 
oder wenn der Machsche Winkel & sehr klein, die Machsche Zahl also sehr groß ist. 

Die GIn. (29a) und (29b) haben die Eigenschaft, daß die Koeffizienten von w;— w, und 
ws — w,lediglich von w abhängen, denn auch c ist Funktion von w allein. Die Gl. (29) können wir 
also so deuten, daß wir den Punkt P, genau so wie beim ebenen Problem finden, indem wir von 
den Punkten 

Sg — 


P;(w, d,+ y 


ausgehen. Wir erhalten P/ und P, aus P, und P, durch Verschieben längs Kurven w = konst., 
also längs Kreisen im Hodographen. Zur Berechnung der Korrekturen 


ge a 5 Sg — Sg. ee 
Isindsin «) und P- (1% Zune x) 
i Y 


83 —— 


Ss A sin Psina; Br ,y sin dsina 
brauchen wir außer den Längen s;— $), $s3— Ss; und y auch die von den Zustandswerten abhängige 
Größe sind sina. Wir können uns im Zustandsdiagramm Linien sin #sin x = konst. ein für 
allemal eingetragen denken, so daß wir sofort zu jedem Punkt, den wir in der Zustandsebene 
eintragen, das Produkt sind sina angeben können. Dies haben wir im folgenden nicht getan, 
um nicht zuviel Kurven im Zustandsdiagramm zu haben. Wir haben es vorgezogen, die Kreise 
w= konst. und die Geraden 9 = konst. zu zeichnen und am Rande außer den Werten von ® 
auch sin 9 und außer der bekannten Busemannschen „p-Zahl‘“ [2] auch sin & aufzutragen. 
So können wir also bei’jedem Punkt die Faktoren sin & und sin d leicht ablesen (Bild 8). 

Das Verfahren erster Ordnung ist nun außerordentlich einfach. Wir suchen die Punkte 
P, und P,im Hodographen auf, bilden dort die Charakteristikenrichtungen genau entsprechend 
Bild 7 und finden so P,in der Störungsebene. Dieser entnehmen wir nun die Abstände s,— S| 
und %,, ferner entnehmen wir dem Hodographen die Werte von sin & und sin d im. Punkte P, 
und finden nun durch Multiplikation und Division die Korrektur. Ist diese positiv, so wird P, 


= G = TE 
ee 
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ne “2 sin dsina in negativer d-Richtung vorzunehmen, um nach P} zu gelangen. Von P/und En 


ausgehend finden wir nun genau so wie beim ebenen Problem P, im Hodographen. 


7 
N 


Bild 8. 
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Eine Schilderung des Verfahrens zweiter Ordnüng erscheint uns hier nach unseren allgemeinen 
Ausführungen im vorhergehenden Abschnitt kaum mehr nötig. 

Eine kleine Schwierigkeit ergibt sich beim achsensymmetrischen Charakteristikenverfahren 
aus der Tatsache, daß die Korrektur auf der Achse unbestimmt wird, weil dort sowohl y als auch 
sin ö verschwindet. Wir sind, wenn der Punkt P, auf der Achse liegt, so vorgegangen, daß wir 
zunächst für diesen die Korrektur Null gesetzt und unter dieser Annahme Ort und Zustand des 
Punktes P, bestimmt haben. Dann haben wir die Korrektur im Halbierungspunkt der. Strecke 
‚P, P, bestimmt und mit dieser Korrektur weitergerechnet. Dieses Vorgehen führte stets zum Ziel. 
N ur müssen wir dabei beim Verfahren erster Ordnung einmal und beim Verfahren zweiter Ordnung 
zweimal iterieren, wobei wir allerdings die Werte von P, nach dem ersten Iterationsschritt un- 
geändert lassen können. Da auf der Achse ohnehin 
nur wenig Punkte liegen, kommt es auf einen Itera- 
tionsschritt mehr oder weniger nicht an. 

Wir können die Korrekturen $, und S, auch 
der Strömungsebene allein entnehmen. Diesempfiehlt 
sich vor allem dann, wenn uns eine Zeichenmaschine 
zur Verfügung steht, mit welcher man Winkel einfach 
übertragen kann (Kullmann - Zeichenmaschine) 
und bei der womöglich außerdem die Zentimeterein- 
Bild 9. Bestimmung der Korrektur bei adiabatischer teilung im Kreuzungspunkt der ‚inneren Begrenzun- 

achsensymmetrischer Strömung. der senkrecht aufeinander stehenden Lineale 
eginnt. 

Beim Verfahren erster Ordnung ist (s,— sı) sin & der, Normalabstand des Punktes P, von der 


Geraden, welche im Punkte P, die Strömungsrichtung hat (siehe Bild 9). ns ist jenes Stück 
dieser Geraden, welche zwischen ihrem Schnittpunkt mit der z-Ac re j 

Das Verhältnis dieser beiden Strecken gibt die Korektne en ee 
Die Strecke (s,;— 8) sina können wir bei einer entsprechend eingerichteten Zeichenmaschine 
sehr leicht finden. Wir stellen mit einem Lineal die Strömungsrichtung ein und legen dieses i i 
Punkte P, an. Nun lesen wir auf dem rechtwinklig dazu stehenden Lineal den gewünschten Ab. 
stand ab, indem wir dieses durch den Punkt P, gehen lassen. Wir können also bei einer geeionete Ä 
Zeichenmaschine die beiden Strecken ablesen, ohne irgendwelche Hilfslinien in die , x 
ebene einzutragen. Lediglich die Frage nach dem Vorzeichen der Korrektur ist noch offen. E & 


im Hodographen auf einem Kreis um den Ursprung in positiver ®-Richtung nach P/ verschoben. FR 
Entsprechendes machen wir für den Punkt P, Hier ist die Verschiebung bei positivem ; 
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und Y ‚ist stets positiv. (83— 5;) ist ebenfalls positiv, wenn wir in Strömungsrichtung rechnen, 
was jaim allgemeinen der Fallist. Das Vorzeichen der Korrektur 8, ist aleo unter diesen Umständen 
stets gleich dem Vorzeichen von 9. Bei 8, ist das Vorzeichen jenem von ® gerade entgegengesetzt. 


Beim Verfahren zweiter Ordnung können wir die Korrektur 8, und 8, in ähnlicher Weise - 
bestimmen. Nehmen wir an, wir hätten bereits eine erste Näherung: P/ von P, sowohl in der Zu- 
standsebene als auch in der Strömungsebene. Diese hätten wir mit dem eben besprochenen Ver- 
fahren erster Ordnung gewonnen. Nun bestimmen wir den Punkt P,im Ortsplan, indem wir die 
Charakteristikenrichtungen an Zwischenpunkten von P/P, und P/P, des Hodographen ent- 
nehmen und diese von P, und P, aus im Ortsplan auftragen. Um die Korrektur des Punktes P, 
noch einmal genauer zu erhalten, nehmen wir die Geschwindigkeitsrichtung im Zwischenpunkt 
von P, und P/ in der Zustandsebene. Ziehen wir von ?, mit dieser Richtung im Ortsplan eine 
Gerade, so ist der Normalabstand des Punktes P, von dieser (s3— s;) sin a in zweiter Näherung. 
Ziehen wir eine Gerade mit der gleichen Richtung vom Halbierungspunkt der Strecke PP, des 


Ortsplanes zur x-Achse, so ist die Länge dieses Geradenstückes En in zweiter Näherung. Bilden 


wir das Verhältnis dieser beiden Strecken, so erhalten wir also die Korrektur in zweiter Näherung. 
Entsprechendes gilt wieder für den Punkt P,. Ser 

Das anschließend wiedergegebene Beispiel zeigt, daß diese Korrektur in zweiter Näherung 
sich von jener der ersten Näherung nur ganz unwesentlich unterscheidet. Der Grund liegt darin, 
daß sich sin & auf der Strecke P, P, des Ortsplanes nur wenig ändert. Da ferner 9 auf der kleinen 
Strecke P, P, sich im wesentlichen linear ändert und das gleiche für y gilt, so kommt eine Änderung 
des Verhältnisses = r) bei der Bildung der zweiten Korrektur nur dadurch zustande, daß der 
Punkt von P/ nach P, verrückt, was aber im allgemeinen kaum ins Gewicht fällt. Für die achsen- 
symmetrische adiabatisehe stationäre Überfallströmung bringt also im allgemeinen das Ver- 
fahren zweiter ‘Ordnung beim Aufsuchen des Ortes von P, in der Strömungsebene wohl etwas; 
die Bestimmung des Zustandes in zweiter Ordnung kann man sich aber anscheinend stets ersparen. 
Dies zeigt auch das folgende Beispiel. 

Wir geben nun vier Berechnungen einer kleinen Aufgabe wieder. Es handelt sich um das 
Ausströmen eines Gases aus einem Rohr mit Parallelströmung in einen Raum mit Unterdruck. 
Im Rohr herrscht dabei eine Geschwindigkeit, welche etwa 1,17mal so groß ist als die Schall- 
geschwindigkeit. Die Gitterpunkte wurden ausgehend von der Rohrströmung auf folgende Arten 
ermittelt: 

1. Mit dem Verfahren erster Ordnung. 

2. Mit dem Verfahren erster Ordnung für die Zustände und dem Verfahren zweiter Ord- 
nung für den Ort. 

3. Mit dem Verfahren zweiter Ordnung für Zustand und Ort. 

Bei diesen ersten drei Methoden wurde die Korrektur stets auf geometrischem Wege ermittelt. 

4. Auf demselben Wege wie bei der zweiten Methode, nur wurde als erste Näherung das 
Verfahren zweiter Ordnung der ebenen Strömung gewählt. Die Korrekturen wurden 
ferner auf arithmetischem Wege bestimmt. 


Im Bild 10a, b und c werden die Resultate jeweils mit der genauesten Methode 3 verglichen. 
Zahlentafel 1 zeigt die auf den verschiedenen Wegen errechneten Werte in dem mit @ bezeich- 


Zahlentafell. 


Werte im Punkte Q 


Method Zeit 
ia P-Zahl w[c® v 

I Ordaunsstürrden Orb: 17-7 7 en 943,4 2,120 2399 3,6 min 
1. Ordnung für den Zustand: - : +» 

72..Ordnung für-den Orb... "we. 941,0 2,150 27207 4,0 min 
1. Ordnung für_den Zustand.  * 2 2 = e 

Bee Ordrung Türkden- Ort. a an 941,4 2,148 27° 6,0 min 
2. Ordnung für den Zustand: : 

BEBeSIchen Torte wer a ee ee ve 942,6 2,130 210710% 

Ta De ya eg ET EI re RR EE 940,9 2,158 280 20’ 


c* „kritische Schallgeschwindigkeit‘. 


Seite in einem Verfahren 
„Zeiten angegeben. Die Zeiten 
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Bild 10a. Vergleich des Ma chschen Netzes von Methode 3 ( 


Bild 10b. Vergleich des Machschen Netzes von Methode3 ( ) und Methode 2 (— — — —). 
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Bild10c. Vergleich des Machschen Netzes von Methode 3 ( 


) und Methode 4 (—— — —). 


ten gewonnen wurden. Die Rechnungen wurden dabei von zwei ausgebildeten Rechnerinnen 
unseres Instituts durchgeführt, wobei eine an der Zeichenmaschine, die andere am Rechnungs- 
protokoll arbeitete. Diese zweite Rechenkraft ist dabei allerdings nicht voll ausgenützt. Trotz- 
dem empfiehlt sich diese Zusammenarbeit, 
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Die Beispiele zeigen, daß die Ergebnisse durchaus zufriedenstellend sind. Da mit der Zeichen- 
maschine Winkel höchstens auf 10’ genau eingestellt werden können, liegt der Unterschied der 
Methode 2 und 3 gerade an der Grenze der Zeichengenauigkeit. Für den praktischen Gebrauch 
ist also die Methode 2 mit einfacher Bestimmung des Zustandes und doppelter Bestimmung des 
er am meisten zu empfehlen. Sie ist für die meisten Fälle genügend genau und erlaubt rasches 

echnen. 


% 


7. Anwendung des Verfahrens auf die instationäre adiabatische Stromfadenströmung. 


Die Gleichungen für die instationäre adiabatische Strömung haben mit den von uns zu- 
nächst verwendeten Zustandsgrößen der Dichte og und der Geschwindigkeit w nicht jene Symmetrie- 
eigenschaften wie etwa die Gleichungen für die ebene oder achsensymmetrische stationäre Über- - 
schallströmung, wenn wir diese für die Geschwindigkeitskomponenten «u und » schreiben. Dies 
ist auch gar nicht möglich, da o und w verschiedene Dimension haben. Es empfiehlt sich daher, 


anstatt der Dichte eine andere Zustandsgröße einzuführen, welche lediglich von der Dichte ab- 


hängt und die Dimension einer Geschwindigkeit hat. Diese Eigenschaft hat bei adiabatischer 
instationärer Strömung die Schallgeschwindigkeit c. In w und c geschrieben bekommen Euler- 
sche Gleichung und Kontinuitätsbedingung für ideale Gase folgende übersichtliche Gestalt (x das 
Verhältnis der spezifischen Wärmen): ; 


wma eg, (2 
5 5 a0) 
c w, u ar re AN) 
Tan 
Ist F# der Stromfadenquerschnitt, so hat f(x) folgende Bedeutung f (x) a 


Für die Richtungen der Charakteristiken in der x-i-Ebene erhalten wir folgende Werte: 
co yıgeWwirc:. un ee 
Wählen wir also die Zustandsgrößen w und c, so erhalten wir die Charakteristikenrichtungen 
einfach durch Bildung von Differenz und Summe der Zustandsgrößen. 
Für die Änderung der Zustandsgrößen längs den Charakteristiken erhalten wir nach (26a) 
und (26b) folgende Formeln: 


< > : (6; — &%) — (w —w,) = — (83 —Sı) Sin gs 


’ 192% 
(a6) + —) =— (88 — 5) sin y,weof(e) 

Beim Verfahren für das homogene Gleichungssystem, der einfachen Rohrströmung, ist f(x) = 0. 
Für diesen Fall wurden bereits von Schultz-Grunow [3] und von Sauer [4] ein- 
fache Verfahren ausgearbeitet. Schultz-Grunow verwendet w und c als Veränderliche. 
Man kann dabei im allgemeinen auf das Zustandsdiagramm verzichten, in welchem dann die 
Charakteristiken Gerade sind. Sauer geht zu neuen Zustandsgrößen w, q (w, c) über, wobei 
q so gewählt ist, daß die Orthogonalitätsbedingung (23) erfüllt wird. q ist dabei die Konstante 
des Energiesatzes der stationären Strömung. Natürlich kann man auch umgekehrt aus der For- 
derung nach der Orthogonalität diese Zustandsgrößen ableiten. Während bei Schultz- 
Grunow die Richtungen einfach durch Addition und Subtraktion der Zustandsgrößen gewonnen 
werden, werden bei Sauer die Richtungen der Charakteristiken genau in Analogie zum 
Prandtl-Busemann- Verfahren gebildet. 


Die Gleichungen (32) haben nun die angenehme Eigenschaft, daß die Koeffizienten der 
Klammerausdrücke (c;— C,), (wa— wı), (C3— 65) und (w;— w,) von den Zustandsgrößen nicht 
abhängen. Wir können beim Verfahren für die Stromfadenströmung den Punkt P, also so finden, 
daß wir die Punkte P, und P, im Zustandsdiagramm bei konstantem c von P, (w,, c,) nach PR 
Do, + (83— s)siny, wef(x),c] und von P,(w,c,) nach P, [w,— (s3- 8;) Sinyz wef(®), 65] 
verschieben. Die beiden Korrekturen sind dabei sehr einfach zu bilden. (s;—s,) sin y, und 
(s3;—8,) sin y, sind die Projektionen der Strecken PÄP,und P, P, in der Strömungsebene auf die 
t-Achse. Diese sind mit w und c und der Ortsfunktion f(x) zu multiplizieren. 


370 Oswatitsch, Über die Charakteristikenverfahren der Hydrodynamik pa, a5j27 Na NovılDen 1oar h 


Die Verschiebungen sind bei der Wahl von w und c als Zustandsgrößen parallel zur w-Achse 
vorzunehmen. Bei der Wahl der Sauerschen Zustandsgrößen w und q sind die Linien c=const. 
in das Charakteristikendiagramm des homogenen Systems einzutragen. Die Linien ce = const. 


w* e® 


sind Parabeln, welche aus der einhüllenden Parabel FIRE En a durch Parallelverschieben 


in Richtung der g-Achse entstehen. Län 
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gs dieser Parabeln sind also die Punkte um das Korrektur- 
stück auf der w-Achse zu verschieben (Bild 11). Der : 


Punkt P, wird dann genau so wie früher durch Fort- 
schreiten auf den Charakteristiken des homogenen 


Systems gefunden. Das Verfahren bleibt in seinen 
Grundzügen im übrigen dasselbe und wir können es. 


uns wohl ersparen, die näheren Einzelheiten noch weiter 
zu erläutern. 4 


8. Zusammenfassung. 


Für den Fall der ebenen adiabatischen statio- 
nären Überfallströmung und der instationären adia- 
batischen Rohrströmung gibt es besonders einfache 
Charakteristikenverfahren, die darauf beruhen, daß 
die Charakteristiken in der Zustandsebene unabhängig 
von den speziellen Anfangsbedingungen festliegen. 
Will man einen Punkt berechnen, der mit zwei be- 


BL 2 2 dl Be a a 


kannten Punkten je eine Charakteristik gemeinsam 
hat, so braucht man zur Ermittlung des Zustandes im 
gesuchten Punkt lediglich die bekannten Punkte in 
der Zustandsebene einzutragen. Der gesuchte Punkt 
ist nun leicht zu finden, weil er auch in der Zu- 
standsebene mit den bekannten Punkten je eine Cha- 
rakteristik gemeinsam hat. Zu diesem Zweck zeich- 
net man die Charakteristiken in der Zustandsebene 

ein für allemal ein. In dieser Arbeit werden auch 
für Fälle, bei denen die Charakteristiken im Zustandsdiagramm nicht von vornherein fest- 
liegen, einfache Methoden entwickelt, indem auf eins der einfachen Verfahren mit festliegen- 
den Charakteristiken im Zustandsdiagramm zurückgegriffen wird. Die Bestimmung des un- 
bekannten Punktes im Zustandsraum geschieht nun so, daß die beiden vorgegebenen Punkte 
zunächst um bestimmte Korrekturwerte verschoben werden und dann die Auffindung des ge- 
suchten Punktes in genau derselben Weise wie bei den einfachen Verfahren erfolgt. Die Be- 
stimmung des gesuchten Punktes im Ortsplan geschieht in genau derselben Weise wie bei den 
einfachen Verfahren. 


Die hier geschilderte Methode ist auf alle in der Hydrodynamik bisher bekannt gewordenen 
Charakteristikenverfahren anwendbar, auch für die Strömungen zwischen Verdichtungsstößen. 
Auf die Berechnung der Verdichtungsstöße selbst wurde nicht eingegangen, da sich hierbei gegen- 
über den bekannten Arbeiten keine neuen Gesichtspunkte ergeben. Eine Vereinfachung der 
Bereehnung von Verdichtungsstößen ist im übrigen nicht von jener Wichtigkeit, wie eine Ver- 
einfachung der Berechnung der gewöhnlichen Felder, da diese Rechnungen im allgemeinen be- 
deutend überwiegen. 


Nach unseren Erfahrungen ist die Arbeitsweise mit der hier geschilderten Methode sehr 
zufriedenstellend. Die Rechnungen und Konstruktionen können von guten Rechenhilfskräften 


ohne weiteres ausgeführt werden, wobei sich im allgemeinen eine Zusammenarbeit zweier Rechen- 
kräfte empfiehlt. 


(s3-5:)5in y We fir Go-SI5n 72 WC fer) 


Bild 11. Bestimmung des Zustandes in P, 
bei instationärer Stromfadenströmung. 
(—— — Linien ce = const) 
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KLEINE MITTEILUNGEN 


Eine einfache Konstruktion zur Auffindung des 
Trägheitsradius für eine beliebige Achse. 


Mit y und z als Hauptträgheitsachsen wird das Träg- 
heitsmoment einer Querschnittsfläche für eine unter 
dem Winkel gegen die y-Achse geneigte Gerade dar- 
gestellt durch 

g=Jy'cos®?o+ Js Sind... 1). 
Führt: man in (1) die Trägheitshalbmesser gemäß ihren 
Definitionsgleichungen J = F-i? ein, so wird 

ig =iycos?p + izsin?p en (2) 
Gl. (2) kann in verschiedener Weise zeichnerisch aus- 
gewertet werden!). Sind a und b die Halbachsen einer 
Ellipse mit der Mittelpunktgleichung 


so folgen bekanntlich die zur Ursprungsgeraden 
z=y'tgp 
parallelen Ellipsentangenten (siehe Bild) zu 


z=ytgp ti +tatlp ..... 


Aus (4) kommt für das Lotn vom Koordinatenursprung 
auf die Ellipsentangenten der Wert 
n = C08Yp' yB+ attg?p = y® cos?p + a?sin?p (B). 
Bringt man weiter die Ursprungsgerade 2 = y'tgp 
mit’der Ellipse im Punkt P zum Schnitt (siehe Bild 2), 
so kommt aus (3) mit Uf=r 
een > 

cos?p 62 cos?p + a?sin?p, 


ı) Man vgl. z.B.: A. Föppl: Vorlesungen über Techn. 
Mechanik, Bd.III. München und Berlin 1940, 8.107. 


» oder 


a? b2 


oder : 
b? cos?p + a?sin?p = are 


(0% 


Die Gl. (5) und (6) werden mit (2) identisch, wenn man 
d«=iz und b=iy i 


- macht.. Dieso erhaltene Ellipse wird bekanntlich Träg- 


heitsellipse genannt. Aus (6) kommt mithin 


»2. +2 
«2 3 ae Uydz 
= 2 ae 
I = 1ycos?p + izsing = Pr 
oder E Y 
Io = en a ee or 


In Bild 2ist gezeigt, wie man Gl. (7) zeichnerisch aus- 
wertenkann. Bringt man die Ursprungsgerade mit dem 
Kreis a über der kleinen Halbachse im Punkt Z zum 
Schnitt und zieht U # || A P, so wird 


OP:OE=04:00 


Eine zweite Möglichkeit folgt, wenn man die Ursprungs- 
gerade mit dem Kreis 5 über der großen Halbachse im 
Punkt F zum. Schnitt bringt und DFB zieht. 
Dann wird y 


OF:OP=0D:ÖB 
oder NEIL 
OP Tr 


Es kommt mithin 
ip = U [9 02% 
Man wird die erste Konstruktion vorziehen, wenn 
die Ursprungsgerade näher bei der z-Achse liegt, die 
zweite hingegen, wenn die Ursprungsgerade näher 
bei der y-Achse liegt. 
Nun sind aber die beiden zeichnerischen Verfahren 


‚gemäß Bild 1 und Bild 2 an das Zeichnen der Trägheits- 


ellipse gebunden und: wesentlich von deren Zeichen- 
genauigkeit abhängig. \ ; 


Bild 3. 


Im folgenden soll nun gezeigt werden (siehe Bild 3), 
wie man den Trägheitsradius für eine beliebige Achse 
ganz einfach zeichnerisch bestimmen kann. Bringt 
man die Ursprungsgerade mit dem Kreis (r, =iy) im 
Punkt A und mit dem Kreis (r, = iz) im Punkt B zum 
Schnitt, macht weiterhin 


00C=UB-sing=iz'sinp, 
UD=OA:c08p=iy:c08Q, 
so wird 
UD2 = iy cos? + 1zsin?p = u 
oder 
UÜD = ip. 


x - Zeichnen der 'Trä 
schnell und einfac ee 
Beatin 


> 


m als Parameter in = Form 3 | | re =. DE i 
: tn er Machen wir nun ; re 
Um nun die Gleichung der Umhüllungskurve zu erhal-  BOIAF, ZDIBE, 
‘ ten, hat man in bekannter Weise aus den beiden go wird Re 
Ben, 00:08-04:09, IR:OD=0B:0R 
anne mn (8a), und damit h ER “ 55: 3 
en) Et a =: 3 ab) OB- I iv ae 


“1 TR m 
ae me a er struktion der Scheitelkrüm 


Diese hstruktion hat: den 
wir nun igfür einen andere 


'kel@, so erhalten wir eine Schar von Geraden und wir 
fragen nun nach der Umhüllungskurve dieser C : 


schar,. Im Achsenkreuz y, z lautet die Abschnittsglei- PAST = a 
chung der Geraden 0 D E Be 
’ a) p 0 s 2 ” - b f h ’ 
63 iy C08 @ a izsinp 7 1 Tuer Ne En, (8). Der Volls 1 gl tw 
Setzen, wir werden, wie man auch « 


zeichnerisch finden 
so kommt aus (8) die Gleichung der Geradenschar mit er in He MEI, Bureokan: 


m zu eliminieren. Aus (8b) kommt 


welcher Wert in (8a) rückeingesetzt, den ge- 
suchten Zusammenhang zwischen y und z 
liefert. Es kommt zunächst 


2 als A £ 
02 en) yı & en 
iy iz T- .) E35 

iz 3 
oder nach einfacher Umformung 


Aral... 


Dies ist aber die Gleichung der Evolute einer 

Ellipse. Sind a und b die Halbachsen einer 

Ellipse, so lautet bekanntlich ihre Evolute 
an ®la bz %s 

(ee) -- )) Sal (9a). 


& 


NIX ZSEN 
NR IS 


A IND FON 
HAN LP 
7\ RI 


BR 
27 


Bl 2 


ZE 


“ 


In Bild 5 wird schließlich die Evolute aus (9) als Um- 
hüllungskurve der in Bild3 gezeichneten Geraden- 


Er 
2 


r HR 


} 


schar für vg gezeichnet. Bo 

Bild4. Karlsruhe. Johann Fadle. ER 
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